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RESUMEN 

En este articulo se presentan dos metodos de ultima generacion orientados a la 
evaluacion del dano sismico esperado tanto en edificios individuales como en grandes 
ciudades. El primer metodo, al que nos referimos como metodo del mdice de vulnerabilidad, 
define la accion sismica en terminos de intensidad macrosismica y caracteriza la vulnerabilidad 
de los edificios en terminos de un indice; las funciones que permiten estimar el dano para una 
intensidad y un indice de vulnerabilidad dados son funciones semiempiricas basadas en datos 
de dano observados en terremotos pasados. El segundo metodo es teorico y se fundamenta en 
el modelado y analisis mecanico del edificio que se caracteriza mediante su espectro de 
capacidad, a partir del que se obtienen las curvas de fragilidad; la accion sismica se define 
mediante el espectro de demanda que se obtiene a partir de la reduccion, por comportamiento 
inelastico del edificio, del espectro de respuesta 5% amortiguada del acelerograma que define 
el terremoto; el cruce entre los espectros de capacidad y de demanda define el punto de 
maxima respuesta del edificio cuando sufre el terremoto; este punto constituye el input en las 
curvas de fragilidad para estimar las probabilidades de ocurrencia de los diferentes estados de 
dano; nos referimos a este segundo metodo como metodo del espectro de capacidad. El 
metodo del mdice de vulnerabilidad es muy eficiente en la caracterizacion de la vulnerabilidad 
permitiendo incorporar, de forma sencilla, informacion complementaria de edificios del mismo 
tipo, pero la definition de la accion sismica en terminos de intensidad es pobre. La 
caracterizacion de la accion sismica es mas fina en el metodo del espectro de capacidad pero 
la definicion de la fragilidad del edificio es computacionalmente costosa ya que requiere 
sofisticados analisis estructu rales. Asi pues, en una primera parte se establece el marco teorico 
de los estudios de riesgo sismico detallando la definicion del terremoto, del edificio y de las 
funciones de dano para ambos metodos. En una segunda parte se efectua la aplicacion a la 
ciudad de Barcelona en Espana, donde se ha podido construir una extensa y depurada base de 
datos que contiene las caracteristicas de situacion, geometricas, estructurales y constructivas 
de los edificios de viviendas. Con ambos metodos se analizan dos escenarios sismicos 
definidos en terminos de intensidad macrosismica y en terminos de espectros de respuesta. El 
primer escenario corresponde a un terremoto historico ocurrido a 25 km al norte de la ciudad. 
El segundo corresponde al terremoto con una probabilidad de ocurrencia del 10% en 50 afios, 
es decir, aquel con un periodo de retorno de 475 anos. Nos referimos al primer escenario como 
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escenario determinista y al segundo como escenario probabilista. A pesar de las diferencias 
metodologicas, los escenarios de dano obtenidos son consistentes y ponen de manifiesto la 
validez y robustez de ambos metodos. Las ciudades que, como Barcelona, estan situadas en 
regiones de peligro sismico entre bajo y moderado no suelen tener conciencia del peligro 
sismico por lo que la ausencia total de precauciones de proteccion sfsmica conduce a una 
elevada vulnerabilidad de sus edificios y, consecuentemente, a un elevado riesgo. 



ABSTRACT 



Two last generation seismic risk evaluation methods are discussed in the paper, useful 
both for individual buildings and cities. The first method, known as vulnerability index method, 
defines the seismic action in terms of the macroseismic intensity and characterizes the 
vulnerability of the buildings in terms of an index. The functions which allow evaluating damage 
for a given intensity and a vulnerability index are semi-empiric, based on observed damages 
during past earthquakes. The second method, namely the capacity spectrum method, is a 
theoretical one and it is based on structural models and on mechanical analyses of buildings, 
which are characterised by their capacity spectra, from where fragility curves are developed. 
The seismic action is defined by means of a demand spectrum obtained by reducing, in function 
of the inelastic behaviour of the building, of the elastic response spectrum corresponding to a 
damping ratio of 5%. Crossing the capacity spectrum with the demand spectrum, the building’s 
maximum seismic response point is obtained (performance point). This point is the input in the 
fragility curves in order to estimate the occurrence probabilities of the different damage states. 
The vulnerability index method is very efficient when characterising vulnerability, as it allows 
incorporating quite easily additional information on similar buildings; nevertheless, it uses a 
definition of the seismic action in terms of intensity which is quite poor. On the contrary, this 
seismic action definition is quite accurate in the capacity spectrum index, although the definition 
of the fragility of the buildings is computationally expensive, as it requires sophisticated 
structural analyses. In the first part of the article, the theoretical frame of the seismic risk 
evaluation is established, giving a detailed definition on the earthquake, building and damage 
functions for both methods. In the second part, the evaluation methods are applied to the city of 
Barcelona (Spain), where an extensive and depurated data base was built, which contains 
characteristics about the location, geometry, structure and construction of dwelling. Two seismic 
scenarios, defined in terms of macroseismic intensity and response spectra, are analysed using 
the two methods. The former scenario corresponds to a historic earthquake that took place 25 
km north of the city. The latter is an earthquake with an occurrence probability of 1 0% in 50 
years, that is, an earthquake with a return period of 475 years. The former scenario will be 
referred as the deterministic scenario, while the latter will be called probabilistic scenario. In 
spite of the differences in methodologies, the damage scenarios obtained are consistent and 
show the validity and robustness of both methods. Cities such as Barcelona, which is located in 
a low to moderate seismic risk region, are not usually aware of the seismic risk. That is why the 
lack of adequate seismic protection measures leads to higher vulnerability of buildings and, 
thus, to higher risk. 



1. INTRODUCCION 



Los estudios de riesgo sismico se han convertido en una herramienta excelente para 
predecir y gestionar los efectos catastroficos de los terremotos en ambientes urbanos. Sin duda 
alguna, el siglo XX, con sus grandes luces y sombras, ha contemplado el crecimiento de la 
poblacion mundial y la transformacion de la sociedad agrfcola y rural en industrial, tecnologica y 
urbana. Las grandes aglomeraciones humanas y el estilo de vida moderno aumentan el riesgo 
de la humanidad. En el mundo, durante este siglo XX han ocurrido mas de 1100 terremotos 
fuertes que han causado mas de millon y medio de vfctimas. Ejemplos de terremotos que 
produjeron grandes perdidas, solo en el ultimo cuarto de siglo, son los de Rumania (1977), 
Chile (1985), Mexico (1985), Armenia (1988), Estados Unidos (1987, 1989, 1994), Colombia 
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(1983, 1994, 1999), Japon (1995), Turquia y Grecia (1999), India (2001), Iran (2004) y el 
terremoto de Indonesia (2004), con mas de 150.000 victimas producidas por el tsunami 
generado. 

La sociedad actual vive en la cultura del riesgo y conoce relativamente bien como 
predecir, cuantificar y gestionar los riesgos de periodos de recurrencia cortos, e incluso, las 
sociedades desarrolladas que conviven con los terremotos, como por ejemplo Japon y el oeste 
de America del Norte, han aprendido a minimizar los efectos daninos de los terremotos sobre 
los edificios y la poblacion. Cada vez mas, la ocurrencia de terremotos pone de manifiesto la 
diferencia clara del patron de danos entre paises desarrollados y los paises en vias de 
desarrollo. En la primera mitad del siglo XX aproximadamente el 75 % de las victimas por 
causa sismica se producian en paises en vias de desarrollo mientras que en la segunda mitad 
menos del 10% se producen en paises desarrollados. Asi pues las naciones en vias de 
desarrollo y otras, cuya sismicidad es moderada o con largos periodos de recurrencia, siguen 
siendo victimas de las catastrofes naturales. En el primer caso, como por ejemplo El Salvador, 
su ritmo de crecimiento es inferior a las intensas, repetidas y enormes perdidas economicas y 
humanas producidas por el impacto de fenomenos como huracanes, ciclones y terremotos, lo 
que, unido a sus escasos recursos, les impide salir de la rueda fatidica del subdesarrollo. En el 
segundo caso, la baja o moderada sismicidad o los largos periodos de retorno de grandes 
terremotos producen una nula preocupacion por el riesgo sismico, una ausencia total de 
memoria historica y un descuido de cualquier precaution sismica. En ambos casos, ello induce 
un importante aumento de la vulnerabilidad del entorno edificado y del riesgo sismico. 
Ciertamente la principal causa de perdidas humanas y economicas que producen los 
terremotos es el comportamiento sismico inadecuado de las estructuras. Concretamente, en el 
siglo XX el colapso de edificios produjo el 75% de las perdidas humanas por terremotos. El 
progreso en el diseno de estructuras se aplica a las estructuras nuevas y en menor medida a la 
rehabilitation de estructuras existentes; pero existen muchas mas estructuras antiguas que 
edificios de nueva en construction. La unica forma de reducir el riesgo es reducir la 
vulnerabilidad de los edificios. Sin duda el analisis del riesgo debe abordarse desde una 
perspectiva multidisciplinar. Los aspectos geologicos, sismologicos e ingenieriles permiten 
estimar el peligro y el dano fisico directo esperado. Pero este no es completo si no se incluyen, 
ademas otros aspectos como el numero de victimas y otros factores sociales, organizacionales 
e institucionales, relacionados con el desarrollo y la capacidad de respuesta de las 
comunidades. Pero el primer paso es evaluar el dano fisico directo. Este trabajo se refiere a la 
estimation del dano sismico directo en edificios. En el se describen y aplican dos modelos 
avanzados de evaluation de escenarios de dano entendidos como la prediction de la escena 
creada por la ocurrencia de un determinado terremoto en una poblacion. El primer metodo 
describe la action mediante la intensidad macrosismica y caracteriza el edificio mediante un 
indice de vulnerabilidad; las funciones y matrices que permiten estimar el dano esperado ante 
una action sismica se basan en la escala de intensidad EMS’98 (Grunthal, 1998). El segundo 
metodo describe la action sismica en terminos de espectros de respuesta y el edificio en 
terminos de espectros de capacidad. La evaluation del dano esperado se efectua mediante 
curvas de fragilidad. 



2. EL RIESGO SISMICO. MARCO TEORICO 



McGuire (2004) define el riesgo sismico como la probabilidad de que una comunidad 
experimente una perdida o de que su entorno edificado sea danado por causa de terremotos. 
Normalmente estas probabilidades representan el grado de perdida o dano que es igualado o 
excedido en un periodo de tiempo. Esta perdida o dano debe ser cuantificable y cuantificada y 
puede ser un coste economico absoluto, el numero de victimas en una ciudad o region o el 
coste de reparation de una instalacion o servicio. El riesgo sismico es calculado mediante un 
conjunto de terremotos, el coste o dano asociado y la correspondiente probabilidad de 
ocurrencia o de excedencia. Por consiguiente un estudio de riesgo sismico requiere la 
consideration de los siguientes aspectos: 1) el terremoto; 2) el edificio o bien expuesto, 3) las 
funciones de dano o funciones de perdida; 4) analisis de riesgo; 5) toma de decisiones. Este 
estudio se refiere a los tres primeros: el terremoto, el edificio y el dano esperado. Los dos 
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restantes requieren estudios multidisciplinares en los que participen los responsables de la 
proteccion civil y de la planificacion y gestion de emergencias. 

2.1. El terremoto: peligrosidad sismica. 

El terremoto define la accion sismica mediante una variable o un conjunto de variables 
C, las cuales pueden ser definidas desde una optica probabilista o determinista. En el caso 
probabilista se define la frecuencia anual de excedencia de un determinado nivel del 
movimiento del suelo. Tambien suele darse en terminos de probabilidad de excedencia durante 
un perfodo de tiempo o mediante su perfodo de retorno. Por ejemplo las prescripciones de las 
normas sfsmicas suelen usar tfpicamente la accion esperada con una probabilidad de 
excedencia del 10 % en 50 anos; que corresponde a un perfodo de recurrencia de 475 anos. 
Los 50 anos corresponden a la vida media de los edificios y el 10% suele ser un nivel de riesgo 
aceptable. Ambos parametros pueden variar en funcion de las caracterfsticas del analisis de 
riesgo. En el caso determinista se establece un terremoto caracterfstico en base al cual se 
efectua el analisis. Habitualmente, se utiliza el terremoto mas fuerte ocurrido en la region, 
tambien llamado maximo terremoto historico. Los escenarios sismicos definen los niveles de 
movimiento esperado del suelo en la zona de estudio, es decir, definen la distribucion espacial 
del parametro o parametros que definen la accion sismica. Es obligatorio que los escenarios 
sismicos incluyan los efectos locales, por ejemplo de suelo, debidos a las caracterfsticas to- 
pograficas o geologicas de la ciudad. Las incertidumbres en la definicion de la accion son altas. 
Particularmente en zonas de peligro sfsmico moderado debido a los enormes perfodos de recu- 
rrencia de los terremotos y, en consecuencia, a la pobreza y crudeza de los datos disponibles. 

Otro aspecto fundamental y controvertido es la caracterizacion del movimiento del 
suelo. Definir la accion mediante un solo parametro, llamese intensidad macrosfsmica o 
maxima aceleracion del suelo, es util y comodo pero insuficiente. La ocurrencia de los 
terremotos y el dano observado han puesto de manifiesto la complejidad de la aceleracion 
sismica y de la respuesta de los edificios. La tendencia actual es definir la accion mediante 
espectros de respuesta compatibles con acelerogramas reales. Con todo, especialmente en 
aquellas regiones donde solo existen datos macrosfsmicos, el amplio perfodo historico cubierto 
por la intensidad macrosfsmica, puede aconsejar basarse en este parametro para estudios de 
riesgo que permiten aproximaciones razonables a la prediction del dano esperado. En zonas 
con registros de aceleracion o con conocimientos de fuente y camino es preferible usar 
espectros de respuesta observados o sinteticos. McGuire (2004) resume los principales 
metodos para la estimation probabilista del peligro sfsmico. Ejemplos de metodos de analisis 
de dano y riesgo que utilizan la intensidad son: el metodo italiano (Angeletti et al., 1988; 
Benedetti y Petrini, 1984; Bendetti et al., 1988; Barbat et al. 1996, 1998; Barbat 1998) y las 
instrucciones ATC-13 (1985) y ATC-25 (1996). Hazus’99 (2002) define la accion en terminos 
espectrales y en el proyecto Risk-UE, que ha proporcionado herramientas para el analisis de 
riesgo sfsmico en ciudades europeas, se proponen metodos basados en la intensidad 
macrosfsmica y en valores espectrales (Milutinovic y Trendafiloski, 2003; Giovinazzi y 
Lagomarsino, 2002; Giovinazzi, 2005; Lagomarsino y Penna, 2003) 

En este trabajo se definen 4 escenarios sismicos para la ciudad de Barcelona. Dos de 
ellos se determinan mediante metodos deterministas y dos mediante tecnicas probabilistas. 
Nos referiremos a estos 4 escenarios como determinista y probabilista de nivel 1 y determinista 
y probabilista de nivel 2. Los de nivel 1 definen la accion mediante intensidades macrosfsmicas 
EMS’98 (Grunthal 1998) los de nivel 2 la definen en terminos espectrales compatibles con el 
eurocodigo (EC-8, 2004). 

2.2. El edificio: matrices tipologicas, vulnerabilidad y capacidad. 

El concepto de vulnerabilidad se relaciona con la capacidad resistente. Dos edificios 
tienen la misma vulnerabilidad cuando sometidos a la misma accion sismica sufren el mismo 
dano. Un edificio es mas vulnerable cuando el dano esperado es mayor. Por lo tanto los 
conceptos de vulnerabilidad y dano sfsmico no son independientes y, frecuentemente esta 
relation puede inducir a una cierta confusion conceptual. De hecho las iniciativas americanas 
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como ATC’s y Hazus obvian el concepto de vulnerabilidad y, para cada tipo de edificio, definen 
directamente funciones de dano del tipo: 

D=f(c\e) (1) 

donde D indica el grado de dano esperado por un edificio tipo e cuando se halla sometido a 
una action sismica cuantificada mediante el parametro c. 

Este tipo de funciones, que por ejemplo el ATC-13 y ATC-25 definen para un gran 
numero de edificios e instalaciones, se utiliza para predecir el dano esperado para un escenario 
sismico concreto. Algunas iniciativas europeas, como por ejemplo el metodo italiano, 
introducen la vulnerabilidad mediante una variable explicita. Cada edificio, estructura, 
infraestructura o servicio puede calificarse mediante una clase de vulnerabilidad (Grunthal, 
1998) o puede cuantificarse mediante un indice que, en general, vale cero para el caso nada 
vulnerable, es decir no se espera dano, y vale 1 para el caso mas vulnerable, es decir para 
intensidades sismicas relativamente bajas se espera un dano importante o incluso el colapso. 
En estos casos, en la determination de funciones de dano, es decir en las funciones del tipo de 
la ecuacion (1), la variable e, que define el edificio, se sustituye por un indice de vulnerabilidad 
Iv. Sandi (1983), Wilches-Chaux (1989), Cardona y Barbat (2000) y Barbat et at. (2006) 
discuten aspectos conceptuales relacionados con la tipificacion de los edificios y su 
vulnerabilidad. Definimos a continuation las variables y conceptos utilizados en este trabajo, en 
el que solo nos referimos a edificios convencionales residenciales. 

BTM (Building Typology Matrix) o matriz tipo logics 

Dada una ciudad o entorno edificado, una BTM o matriz tipologica es una clasificacion 
de los edificios segun clases o tipos. Probablemente el ejemplo mas sencillo de matriz 
tipologica es la de la escala MSK, donde los edificios se agrupan de acuerdo a solo tres clases: 
A, B y C. En el metodo de nivel II los edificios se clasifican mediante una BTM. Para cada tipo 
de edificio se define su capacidad y fragilidad. Los conceptos de capacidad y fragilidad se 
definen mas abajo. 

Clases e Indices de Vulnerabilidad 

Constituyen una forma alternativa de clasificar los edificios. Por ejemplo la escala 
macrosismica europea EMS’98 establece 6 clases de vulnerabilidad: A, B, C, D, E y F. A 
diferencia de la clasificacion de acuerdo a una BTM de tipos constructivos, el concepto, clases 
e indices de vulnerabilidad hacen referencia, generalmente implicita, al dano esperado y, por 
consiguiente, diferentes tipos constructivos pueden estar bien representados por una misma 
clase o indice de vulnerabilidad. En el metodo de nivel I, los edificios se caracterizan mediante 
un indice de vulnerabilidad. La estimation del dano se efectua a partir de la definition de 
funciones de dano del tipo definido en la ecuacion (1). 

Curvas y Espectros de Capacidad 

Las curvas y espectros de capacidad permiten caracterizar el comportamiento 
resistente de los edificios de forma analitica. Aunque un edificio de planta rectangular, cuando 
es sometido a un movimiento no estacionario tridimensional, tiene diferentes caracteristicas 
segun las dos direcciones ortogonales en general, para estimar el dano esperado, el edificio 
puede ser bien representado mediante un modelo sencillo de un grado de libertad. Los 
metodos analiticos de estimation de dano tienen en cuenta el comportamiento no lineal de una 
estructura durante movimientos sismicos fuertes. Una curva fuerza-deformacion cuantifica, de 
forma simplificada, la no linealidad y es caracteristica de la estructura que representa. La curva 
de capacidad muestra la fuerza lateral que se requiere para causar un determinado 
desplazamiento. Una forma conveniente y equivalente de representar la cun/a de capacidad es 
en el formato (Sd, Sa), siendo respectivamente Sd y Sa, las aceleraciones y desplazamientos 
de respuesta o espectrales. 

Este tipo de diagramas muestran la aceleracion espectral asociada a un rango de 
desplazamientos espectrales y se llaman espectros de capacidad. Es frecuente mostrar el 
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espectro de capacidad de forma bilineal simplificada. En esta forma el espectro de capacidad 
se aproxima mediante dos tramos de recta definidos por los puntos de cedencia (Sd y , Sa y ) y de 
capacidad ultima (Sd u , Sa u ). La Figura 1 muestra un ejemplo de espectro de capacidad. Su 
forma bilineal se muestra en trazo discontinue y sobre ella se indican los puntos de cedencia y 
de capacidad ultima. El primer tramo recto que pasa por el origen define el perfodo propio del 
edificio y representa su comportamiento lineal para desplazamientos y aceleraciones 
espectrales pequenos. El tramo entre los puntos de cedencia y de capacidad ultima modela el 
comportamiento no lineal. Cuanto mayores sean las aceleraciones de cedencia y de capacidad 
ultima, mayor es la resistencia de la estructura. Como veremos mas adelante, los espectros de 
capacidad permiten estimar el desplazamiento espectral esperado cuando el edificio sufre un 
terremoto caracterizado por su espectro de respuesta, pero ademas, tambien permite estimar 
curvas de fragilidad y, por consiguiente, el dano esperado. 




Figura 1 .- Espectro (trazo contfnuo) y espectro bilineal (trazo discontinue) de capacidad. Se 
muestran tambien los puntos de cedencia (Sd y , Sa y ) y de capacidad ultima ( Sd u , 
Sa u ). 



2.3. El dano 

Desde mucho antes del desarrollo de la sismologia y de la ingenierfa sfsmica el dano 
causado por los terremotos ha sido el observable mas usado para describir y cuantificar los 
terremotos. Jacobo Castaldi y Robert Mallet se consideran precursores de las escalas de 
intensidad y, en sentido amplio, podemos decir que son tambien precursores de los escenarios 
de riesgo sismico. El primero recogio, en un mapa, los efectos del terremoto de Niza ocurrido el 
20 de julio de 1561, el segundo se desplazo con un equipo de expertos para estudiar los 
efectos del terremoto de Napoles ocurrido en diciembre de 1857. En esta epoca ya era 
conocido que la distribution de los efectos macrosfsmicos de un terremoto pueden ser 
representados mediante Ifneas que delimitan zonas de igual intensidad o zonas en las que los 
efectos causados por los terremotos son similares. Con todo, la description y cuantificacion del 
dano no ha dejado aun de tener grandes dosis de subjetividad que dificultan su tratamiento 
cientifico y preciso. Existen pues diversas formas de describir y cuantificar el dano. Aquf 
haremos referencia a la forma adoptada en este estudio. En general las escalas de intensidad y 
los metodos empfricos de analisis de riesgo basados en el concepto de vulnerabilidad y en 
observaciones macrosfsmicas definen el dano mediante los llamados estados de dano. Los 
metodos analiticos basados en espectros de capacidad y curvas de fragilidad heredan la 
consideration del dano mediante grados o niveles de dano. Con todo, en ambos casos, el dano 
suele tambien cuantificarse mediante diferentes parametros o indices de dano cuyo significado 
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es tambien diverso. Las funciones de dano (ecuacion (1)) deben permitir predecir, para cada 
edificio y para cada escenario sfsmico, el dano esperado. Describimos a continuation los 
conceptos y formalismos empleados en nuestro caso de aplicacion. Distinguiremos el metodo 
de nivel I y el metodo de nivel II. El primero es representative de los metodos que usan indices 
de vulnerabilidad y funciones de dano empfricas basadas en escalas macrosfsmicas. El 
segundo es representative de los que usan espectros de capacidad y curvas de fragilidad. 



2.3.1. Estados y distribuciones de dano 

La escala de intensidad MSK-64 (Sponheuer y Karnik, 1964) establece cinco clases de 
dano no nulo. La escala EMS'98, que debe considerarse como una actualization de la MSK- 
64, confirma estos cinco grados o estados de dano. Asf pues, incluyendo el grado de dano nulo 
o sin dano, cualitativamente y cuantitativamente se postulan los 6 grados de dano siguientes: 
0.- sin dano, 1 .- Leve, 2.- Moderado, 3.- Severo, 4.- Extensivo y 5.- Completo. Por otra parte, la 
escala macrosfsmica EMS’98, establece que el dano causado por los terremotos, se distribuye 
siguiendo una distribution binomial. Es decir: 



p (N, d) = — 1)! — d k (1 - df-'- k 
k (N- \-k)\k\ 



k = 0-(N-\) (2 ) 



donde N es el numero de los estados de dano y P k (N,d) define la probabilidad del estado de 
dano k. En nuestro caso N vale 6 y k va de 0 a 5. De esta forma la ecuacion (2) se transforma 
en: 



P,(6,d) = — — — d k (l-d) 5 ~ k k = 0---5 (3) 

* (5 -k)\k\ K} 

d es un parametro que define completamente la distribution y toma valores entre 0 y 1 . Para 
d=0 la probabilidad del estado de dano nulo es la unidad, mientras que para d=1 la probabilidad 
del estado de dano completo es la unidad. Para valores intermedios de d las probabilidades de 
los diferentes estados de dano toman valores entre cero y la unidad. Otro parametro 
frecuentemente utilizado es el llamado grado de dano medio o estado de dano esperado, que 
se define como d*=d (N-1). A partir de la distribution de probabilidad de dano d* se puede 
obtener mediante la siguiente ecuacion: 



d* = ^kP k (6,d) 

k = 0 



(4) 



d* toma valores entre 0 y 5 y representa el estado de dano mas probable. Asf, d=d*= 0 
representa una probabilidad 1 del estado de dano Nulo (0) y 0 para los otros estados. Un valor 
de d=1 corresponde a d*=5, a una probabilidad igual a la unidad para el estado de dano de 
Colapso (5) y una probabilidad nula para los otros estados de dano. Valores intermedios entre 
0 y 1 para d y entre 0 y 5 para d* indican una distribution de dano. En este caso, d* es una 
especie de grado de dano medio o grado de dano mas probable. En cualquier caso lo 
importante es que d*determina de forma unfvoca la distribution de las probabilidades de dano. 
Esta sencilla hipotesis sobre la distribution del dano permite construir distribuciones de 
probabilidad de dano a partir de descripciones parciales o incompletas en las escalas de 
intensidad. Asf por ejemplo la escala MSK-64 describe de la siguiente manera los danos 
sufridos por los edificios de tipo A cuando ocurre un terremoto de intensidad VII. “Muchas 
construcciones del tipo A sufren dahos graves (clase 3) y algunas incluso destruccion (clase 
4)”. Entendiendo por muchas el 50% y por algunas el 5%, podemos predecir la distribution 
binomial probable de forma que la suma de las probabilidades de los estados de dano grave y 
destruccion sumen 55%. La Tabla 1 y la Figura 2 muestran un ejemplo de este tipo de 
construction de la distribution binomial que corresponde a d=0.5553 y d*=2.78. 
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Una distribucion de probabilidad mas versatil es la distribution Beta, cuya funcion de 
densidad viene dada por: 



, . r(f) (x-a) q 1 (b-x)‘ 9 1 



para a<x<b 



(5) 



r(f) es la funcion gamma, ay t) estan relacionados con el intervalo de definicion de la 

distribucion y, por lo tanto, con el numero de estados de dano. t y q permiten modular la forma 
de la distribucion. t esta relacionado con la dispersion. La probabilidad de que se de un nivel 
de dano x, viene dada por la integral entre a y xde la funcion p p (x), es decir: 



i p( x )= \pA z ) d& ( 6 ) 



Tabla 1 .- Completando una distribucion de probablilidad de dano. Ver explicacion en el texto. 



Edificios Tipo A 


Grado de 
dano 


Intensidad 


VII 


VII 


Nulo (0) 




0.02 


Insignificante 

(1) 




0.11 


Leve (2) 




0.27 


Moderado (3) 


0.50 


0.34 


Severo (4) 


0.05 


0.21 


Colapso (5) 




0.05 


d=0.5553, d*=2.78 




Nulo Leve Moderado Severo Destruction Colapso 

Grado de dano 



Figura 2.- Distribucion binomial de probabilidad 
correspondiente a la Tabla 1 . 



En el caso de estados de dano discretos definidos por los niveles de dano k, (k=0,1 ,5) las 

probabilidades de cada estado de dano se definen mediante: 



= Pp(k + l)-PAk) 



Pk = r p 



(7) 



y la probabilidad de que se iguale o exceda un determinado estado o nivel de dano se define 
mediante: 



P(D > D k )= 1 - P p (k) (8) 

Tomando a=0, b=6, t=8 se consigue una distribucion binomial equivalente que se 
ajusta imponiendo que el factor de dano medio o nivel de dano esperado sea el mismo. La 
siguiente ecuacion establece la equivalencia entre las distribuciones binomial de la ecuacion (3) 
y la distribucion Beta definida en la ecuacion (5): 
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q = ?|o.007<f* 3 - 0.052 d* 2 + 0.2875^* j (9) 

donde d es el factor de dano tal como ha sido definido en la ecuacion (4). De esta forma el 
dano esperado y la distribution de probabilidad de los diferentes estados puede definirse 
mediante un solo parametro. 

2.3.2. Curvas de fragilidad 

Una forma equivalente de definir el dano es mediante curvas de fragilidad (Hazus’99, 
2002). Para un tipo de edificio y un estado de dano dado, estas curvas definen la probabilidad 
de que se iguale o exceda el estado de dano considerado en funcion de la intensidad sfsmica. 
Es decir: 



F D k A l ) = P[GD > D k | i] (10) 

T indica la tipologfa considerada, D k es el grado de dano de la curva de fragilidad 
correspondiente e /' es la variable que define la intensidad sfsmica. La curva de fragilidad del 
estado de dano nulo es trivialmente la unidad. Se supone que para los estados de dano no 
nulo, las curvas de fragilidad siguen una distribution de probabilidad lognormal: 



Co 

II 

e 


—In 


( s 


i “ j 


Pels 


e 

\^d,ds ) _ 



( 11 ) 



donde ahora S d es el desplazamiento espectral o parametro de la action sfsmica, S d ds es el 

valor medio del desplazamiento en el que el edificio alcanza un umbral del estado de dano ds, 
I3 ds es la desviacion estandar del logaritmo natural del desplazamiento espectral del estado de 
dano ds y ® es la funcion de distribution acumulativa normal estandar. Por lo tanto, estas 
curvas quedan definidas mediante solo dos parametros: el valor medio S d ds y la desviacion 

tfpica fids- El valor medio define el punto en el que la probabilidad de igualar o exceder el estado 
de dano es igual al 50%, la desviacion tfpica nos da una idea de la dispersion. La siguiente 
ecuacion muestra la funcion de densidad de probabilidad correspondiente a las curvas de 
fragilidad. 



Fr « ASd) =j^j. sr 1 * 



2 



' 1 , Sd^ 2 

— ln^ 

P Sd 



( 12 ) 



donde Sd es el desplazamiento espectral, Sd es el valor medio de la distribution de 
probabilidad, y b es la desviacion tfpica del In Sd . La curva de fragilidad, es decir, la 

probabilidad de que un edificio iguale o exceda el estado de dano considerado, viene dada por 
la integral entre 0 y Sd de la funcion de densidad de probabilidad dada en la ecuacion (12), es 
decir: 



Sd 

P (Sd) = P T Gdi [ GD > Gdi \ Sd] = | F T Gdi (Sd) d (Sd) ( 1 3) 

o 



La Figura 3 muestra un ejemplo de curvas de fragilidad. 
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Desplazamiento espectral 



Figura 3.- Ejemplo de curvas de fragilidad. Se consideran 4 estados de 
dano diferentes del grado de dano nulo o sin dano. 



Finalmente, conocidas las curvas de fragilidad, es posible estimar la probabilidad de 
ocurrencia de cada estado de dano mediante la siguiente ecuacion: 



P[GD=D k \i,T] = P(GD>D k \i,T)-P(GD>D (k+l) \i,T) k = 0-(N-l) (14) 



Para el estado de dano de colapso o k=N, la correspondiente curva de fragilidad 
define su probabilidad de ocurrencia, ya que no hay grados mayores. 



2.4. Funciones y matrices de probabilidad de dano 

El punto clave para una prediccion fiable del dano esperado esta en el ajuste de 
funciones de dano. Estas funciones definen, para cada edificio o servicio, el dano previsible 
cuando el edificio sufre un determinado terremoto. Pueden referirse a cualquier aspecto 
relacionado con el dano, como por ejemplo el dano ffsico directo, los danos inducidos, el dano 
socio-economico, el numero de vfctimas, el volumen de escombros, el numero de familias sin 
hogar o los tiempos de recuperation de un servicio, entre otros. Aquf nos referimos a dano 
ffsico directo esperado en edificios destinados a viviendas. Existen varios metodos para calibrar 
estas funciones de dano. ATC-13 (1985), por ejemplo, usa la opinion de expertos. Aquf 
describimos brevemente el tratamiento de este problema en los dos metodos escogidos para 
este estudio. 



2.4.1. Metodo del fndice de vulnerabilidad 

En este metodo la action sfsmica se caracteriza mediante la intensidad macrosfsmica 
EMS’98 y el edificio mediante un fndice de vulnerabilidad. La information de danos observados 
en terremotos pasados se ha utilizado para obtener la siguiente funcion semi-empfrica 
propuesta por Giovinazzi y Lagomarsino (2002): 



d* = 2.5 



1.0 + tank 



1 + 6.25 Vj - 13.1 
23 



(15) 



La Figura 4 muestra las funciones de dano para diferentes indices de vulnerabilidad. 
Observese como, dada una intensidad sfsmica, el grado de dano aumenta con el fndice de 
vulnerabilidad. Por ejemplo, para una intensidad X el grado de dano mas probable es 0.52 para 
un fndice de vulnerabilidad de 0.1 mientras que es 4.5 para un fndice de vulnerabilidad de 0.9, 
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indicando que, para esta intensidad, los edificios tipificados mediante el primer indice de 
vulnerabilidad apenas sufriran danos mientras que los danos seran muy severos en los 
edificios caracterizados mediante el segundo indice. 



Figura 4.- 

Ejemplo de funciones de 
dano para diferentes 
indices de vulnerabilidad. 
d* se obtiene a partir de la 
ecuacion (15). 




Este grado de dano d* es el grado de dano medio definido en la ecuacion (4), que 
permite recuperar la probabilidad de cada estado de dano, o matriz de probabilidad de dano, a 
partir de la ecuacion (3) o a partir de la distribution Beta definida en la ecuacion (5). De esta 
forma la estimation del grado de dano mas probable y de la distribution de probabilidad de los 
diferentes estados de dano se obtiene de forma directa a partir del indice de vulnerabilidad y de 
la intensidad. Por lo tanto es crucial el afinar bien en la asignacion del indice de vulnerabilidad 
a los edificios. 



2.4.2. Metodo del espectro de capacidad 

El metodo del espectro de capacidad considera solo 5 estados de dano, incluyendo el 
estado de dano nulo. Se considera en este trabajo, que el grado de dano completo, incluye los 
estados de dano extensivo y completo, por entender que cuando el dano es extensivo es muy 
probable que haya que demoler el edificio dado que deja de ser reparable. El dano se 
cuantifica mediante las curvas de fragilidad definidas mas arriba. Se usa una forma simplificada 
de obtencion de estas curvas a partir de los espectros bilineales de capacidad. El espectro 
bilineal de capacidad se define mediante el punto de cedencia (Sd y , Sa y ) y el punto de ultima 
capacidad. (Sd u , Sa u ). La definition de la distribution lognormal que define la curva de 

fragilidad requiere estimar el valor medio o umbral del estado de dano (Sd,) y la desviacion 
tfpica ($). En este trabajo se sigue un procedimiento simplificado propuesto en el proyecto 

Risk-UE. Los umbrales de los estados de dano se definen mediante las relaciones de la Tabla 
2 (vease tambien Figura 5). Para estimar las desviaciones tfpicas se supone que las 
probabilidades de los estados de dano en los umbrales, se distribuyen de acuerdo a una 
distribution Binomial o Beta equivalente. De esta forma, fijando la probabilidad de excedencia 
del umbral considerado sea del 50% se obtienen las probabilidades de los otros estados 
mediante la ecuacion (6). 
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Tabla 2.- 

Estados de dano y umbrales a partir 
de los puntos de cedencia y de 
capacidad ultima de un espectro de 
capacidad. 



Umbral del estado de 
dano 


Estado 
de dano 


~Sd~ l = 0.7 Sdy 


Leve 


— 


Moderad 


Sd 2 = Sdy 


0 


Sd 3 = Sdy + 0.25 (Sdu - Sdy) 


Severo 


Sd 4 -Sdu 


Colapso 




Figura 5.- Umbrales de los estados de dano. 



La Tabla 3 muestra los valores obtenidos. Un ajuste por rmnimos cuadrados de la 
ecuacion 13 permite estimar las desviaciones tfpicas. La Figura 6 ilustra este proceso. Los 
puntos corresponden a los valores de la Tabla 3, las Ifneas continuas corresponden a los 
ajustes obtenidos. 

De esta forma para cada edificio definido por su espectro de capacidad se obtienen las 
curvas de fragilidad. El problema de la estimacion de la matriz de probabilidad de dano 
correspondiente a un escenario sfsmico o terremoto concreto se reduce ahora a obtener el 
desplazamiento espectral esperado, que permite entrar en las curvas de fragilidad y obtener las 
probabilidades de cada estado de dano mediante las ecuaciones (13) y (14). El metodo del 
espectro de capacidad permite realizar esta estimacion mediante la obtencion del llamado 
punto de desempeno (“performance point’) que proporciona el nivel de desplazamiento 
espectral producido en el edificio por el terremoto del escenario considerado (Fajfar, 2002; 
Fajfar y Gaspersic 1996). El ATC-40 (1996) proporciona tres metodos para determinar este 
punto, Milutinovic y Trendafiloski (2003) explican, paso a paso, la aplicacion de estos metodos. 
Todos ellos tienen en cuenta el comportamiento no lineal de los edificios cuando se hallan 
sometidos a fuerzas severas que producen la incursion en su comportamiento no elastico y, por 
lo tanto, requieren la reduccion del espectro de respuesta lineal 5% amortiguado. 



Tabla 3.- 

Distribucion de probabilidades 
de excedencia fijando la 
probabilidad del umbral en un 
50%. 



Pp(1) 


0.50 


0.90 


0.99 


1.00 


Pp(2) 


0.12 


0.50 


0.87 


0.99 


Pp(3) 


0.01 


0.14 


0.50 


0.88 


Pp(4) 


0.00 


0.01 


0.10 


0.50 


d* 


0.91 


1.92 


3.08 


4.09 




Sd 1 Sd 2 Sd 3 Sd 4 

Desplazamiento espectral 



Figura 6.- Curvas de fragilidad ajustadas a los puntos de 

la Tabla 3. 
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La Figura 7 muestra un ejemplo del procedimiento. El edificio se representa por el 
espectro de capacidad mientras que la action sfsmica se representa mediante el espectro de 
respuesta lineal 5% amortiguada. El cruce entre ambos espectros permite efectuar la reduction 
del espectro de respuesta para obtener el espectro de respuesta no lineal, llamado tambien 
espectro de demanda, y el punto de desempeno. El desplazamiento espectral que define el 
punto de desempeno representa el desplazamiento producido en el edificio por el terremoto. 
Este desplazamiento espectral permite entrar en la curva de fragilidad (Figura 3 y Figura 6) 
para, usando las ecuaciones (13) y (14), determinar las probabilidades de ocurrencia de todos 
y cada uno de los estados de dano o matrices de probabilidad de dano (Whitman et al., 1973). 
A partir de estas matrices puede estimarse el grado de dano medio o grado de dano mas 
probable (ecuacion (4)). 




Figura 7.-Ejemplo de estimation del punto de desempeno que define la 
respuesta maxima del edificio para la action sfsmica mediante 
el espectro de respuesta. Tambien se obtiene el espectro de 
demanda o espectro de respuesta inelastica. 

En resumen, en el metodo del indice de vulnerabilidad, la action se cuantifica mediante 
la intensidad macrosfsmica y el edificio mediante un fndice de vulnerabilidad; el ajuste de 
funciones semiempfricas permiten establecer la relation que predice el grado de dano medio 
esperado en el edificio cuando ocurre el terremoto. La hipotesis de una distribution estadfstica 
de las probabilidades de los estados de dano permite estimar las matrices de probabilidad de 
dano. En el metodo del espectro de capacidad, el edificio se representa por su espectro bilineal 
de capacidad y el terremoto mediante su espectro de respuesta lineal 5% amortiguada. El 
espectro de capacidad permite obtener las curvas de fragilidad. El cruce de los espectros de 
capacidad y de respuesta produce el desplazamiento maximo esperado en el edificio sometido 
a la action sfsmica. Entrando en las curvas de fragilidad se obtienen las matrices de 
probabilidad de dano. El metodo del fndice de vulnerabilidad pertenece al conjunto de metodos 
de analisis de dano sfsmico llamados empfricos o semiempfricos ya que se basan en danos 
observados y en consideraciones analfticas sobre los edificios. El metodo del espectro de 
capacidad pertenece a otro grupo de metodos denominados analfticos o mecanicos. Estos 
metodos se fundamentan en el analisis estructural y en la estimation del dano a partir de 
consideraciones mecanicas (Giovinazzi, 2005). Ambos metodos se aplican, a continuation, al 
parque de edificios de viviendas de la ciudad de Barcelona. 



3. APLICACION A LA CIUDAD DE BARCELONA 
3.1. Un poco de historia 

Existen indicios de que hubo habitantes en Barcelona desde el ano 2.500 a JC, con un 
primer poblado datado entre los anos 2.000 y 1.500 a JC. Entre los siglos III y II a JC existfan 
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dos poblados Ibericos. Uno, llamado Barcilo, en un monticulo, conocido como monte Taber, 
donde se halla el nucleo antiguo de la actual ciudad. Otro, llamado Laye, en MontjuTc. Los 
romanos llegaron a Barcelona 218 anos a JC prosperando la colonia del Mons Taber, ya en el 
siglo I a. J.C. Tras ser destruida, la ciudad fue reconstruida, probablemente en el ano 263, y 
fortificada antes de finalizar el siglo IV. Su muralla la convertian en una ciudad segura. La 
ciudad paso a manos musulmanas probablemente el ano 717 6 718. La ciudad mantuvo su 
aspecto durante la epoca Carolingia desde los inicios del siglo IX y hasta principios de siglo XI. 
En el ano 1 .000, Barcelona conserva su nucleo fortificado pero nuevos barrios surgen ya en su 
exterior. En la primera mitad del siglo XIII Barcelona empieza a construir una nueva muralla 
que contiene estas aldeas. En 1.249 Jaime I creo el municipio y las primeras instituciones 
ciudadanas de gobierno de la ciudad. En el ano 1 .333 la peste se cebo en la ciudad causando 
mas de 10.000 muertes. Sucesivas pestes y epidemias azotaron la ciudad de forma 
intermitente hasta los siglos XV y XVI. En 1.359 se construye una tercera muralla. Con todo, la 
zona izquierda de las Ramblas, conocida como el arrabal (el Raval), siguio siendo una zona 
rural que no se edifico y poblo hasta la expansion industrial de 1.780. El siglo XVIII fue un mal 
siglo para Barcelona que tuvo que afrontar importantes modificaciones polfticas y urbanfsticas. 
A este perfodo corresponde la creacion de La Barceloneta. A partir de 1774 se empieza a 
derribar la muralla interior que separaba la ciudad del Raval y se poblaba y urbanizaba este 
sector a la izquierda de la Rambla. En la primera mitad del siglo XIX se intensifica la 
urbanizacion del interior de la ciudad que, progresivamente, se dota del aspecto e 
infraestructuras de una ciudad moderna. Pero lo que da paso al nacimiento de la Barcelona 
actual es el ensanche por el llano. Entre 1858 y 1868 se derriban las murallas de la ciudad y se 
inicia el gran proyecto urbanistico para la construction de la Barcelona moderna, una nueva 
ciudad abierta e industrial. El Ensanche (“I’Eixample’) de Barcelona, fue proyectado por el 
ingeniero lldefonso Cerda (Cerda, 1968). Este proyecto de expansion de la ciudad permitio unir 
el nucleo de la ciudad con los pueblos del llano que daran nombre a los distritos y barrios de la 
ciudad actual. El censo de 1 .900 establece en 553.000 el numero de habitantes de esta ciudad 
ampliada. Entre 1.909 y 1.929 Barcelona experimenta una colosal expansion. De los 587.411 
habitantes de 1.910 se pasa a 1.005.565 en 1930. Entre 1940 y 1960 se produce un 
considerable crecimiento sostenido y la ciudad gana medio millon de habitantes. Este 
incremento de la poblacion es debido principalmente al aumento de la natalidad y al gran 
numero de inmigrantes procedentes del Sur de Espana. A pesar de que la Carta municipal de 
1961 otorga mayor autonomia a la ciudad, durante la alcaldfa de Porcioles (1957-1973) se 
produjo un crecimiento especulativo y desordenado con grandes deficiencias estructurales. En 
1975 fallece el general dictador Franco y, en 1978, se restablece la monarqufa y la democracia 
en Espana iniciandose un mayor control de la sociedad sobre las instituciones y el gobierno de 
la ciudad. En la actualidad Barcelona esta encajada entre Montjuic, Collserola y el rfo Besds, 
que limitan claramente su crecimiento en superficie, dejando como unica alternativa la 
reconversion de zonas deprimidas, en desuso o industriales. Los ultimos ejemplos de 
transformation urbanfstica de la ciudad son la Vila Ollmpica de 1992 y la apertura al mar de la 
Avinguda Diagonal con el Forum de las culturas en el 2004. El distrito del Eixample destaca 
como la principal area de Barcelona donde se acumula una importante poblacion una notable 
actividad economica y gran patrimonio cultural (Garcia Espuche, 1990; Paricio, 2001). 
Actualmente, el municipio de Barcelona se extiende en una superficie de 9.907 Ha. con una 
poblacion de 1.5 millones de habitantes segun el censo de 2001. 

3.2. Los datos 

La informacion utilizada en este estudio ha sido facilitada, en su mayor parte, por el 
Institut Municipal d’lnformatica (IMI) y por el Departament dEstadlstica del Ayuntamiento de 
Barcelona (DEAB). Estas dos instituciones Se complementan en la gestion y analisis de la 
informacion referente a la ciudad. La unidad basica de estudio es la parcela catastral edificada 
o edificio. Para su caracterizacion se han usado tres bases de datos: el catastro, las edades y 
las caracteristicas constructivas. La base de datos del catastro esta bien documentada en el 
informe Infocca (1999) y permite determinar las caracteristicas geometricas, en planta y en 
altura de los edificios de cada parcela. A partir de la correlation de la informacion contenida en 
estas bases de datos se dispone de informacion completa de 70.905 edificios que representan 
el 93,38 % de los edificios del censo del 2001 . De estos edificios se conoce su configuration en 
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planta y en altura, su situation dentro de la manzana, formando parte de un agregado o no, asf 
como la information referente a los edificios colaterales; la tipologfa constructiva de la 
estructura y de los forjados, el ano de construction y el estado de conservation. La calidad, 
representatividad y resolution de nuestra base de datos global, permite efectuar los analisis 
edificio por edificio, pero por claridad, simplicidad, utilidad y, sobre todo, para no perder de vista 
el sentido probabilista de los analisis y conclusiones, los resultados se van a presentar de 
acuerdo a las tres grandes divisiones administrativas definidas y utilizadas por el ayuntamiento 
de la ciudad: Distritos, Barrios y Zonas de Estudio (ZRP). Barcelona consta de 10 distritos, 38 
barrios y 238 ZRP. Toda la documentation recopilada se usara para caracterizar la 
vulnerabilidad, capacidad y fragilidad de los edificios de la ciudad. 

3.3. Los terremotos: escenarios sismicos 

Para definir el movimiento sfsmico esperado de los escenarios considerados se ha 
tenido en cuenta la zonificacion sfsmica de los suelos. Cid (1998), clasifica los suelos de la 
ciudad en 4 zonas. La zona R constituida por afloramientos rocosos paleozoicos y terciarios 
(ver tambien Cid et al. 2001). La zona I corresponde a suelos blandos caracterizados por un 
basamento paleozoico a una profundidad inferior a 350 m sobre el que se halla una gruesa 
capa de materiales terciarios y una capa de materiales deltaicos del holoceno con un espesor 
de entre 25 y 70 m. La calidad de los suelos aumenta en las zonas II y III, que son de transition 
entre los suelos blandos de los deltas y los afloramientos rocosos. 





Figura 8.-Zonificacion sfsmica de los Figura 9.- Funciones de amplification 

suelos de Barcelona (Cid 1 998). espectral de los suelos de 

Barcelona (Cid 1998). 

Los efectos sobre la action sfsmica en terminos de intensidad se han tenido en cuenta 
incrementando la intensidad esperada en el basamento rocoso en media unidad en las zonas 
de suelos tipo II y tipo III y en una unidad en la zona I. En terminos espectrales se han tenido 
en cuenta las funciones de transferencia propuestas por Cid (1998) para los suelos de las 
zonas I, II y III. La Figura 8 muestra el mapa de zonificacion sfsmica de los suelos de la ciudad, 
la Figura 9 muestra las funciones de amplification espectral. 

El escenario determinista en terminos macrosfsmicos, corresponde a un terremoto 
historico ocurrido en 1448 en Cardedeu, a una distancia de 25 km del centra de Barcelona y a 
una profundidad de 7 km (Olivera et al. 2006). Este es el sismo de mayor intensidad epicentral 
ocurrido a menor distancia de la ciudad. La Figura 10 muestra la distribution de intensidades 
esperadas en la ciudad. Se observa el efecto de la atenuacion sfsmica hacia el Sur y los 
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efectos de suelo que amplifican la serial hacia el Este, en el frente de mar. Las intensidades 
varian entre VI en el afloramiento rocoso de la montana de Montjui'c y VIII el el Noreste de la 
ciudad donde coinciden la zona mas cercana al epicentro y los suelos mas blandos. 




Figura 10.- Escenario determinista en 
terminos macrosismicos. 




Periodo, s 

Figura 11.- Escenario determinista en 
terminos espectrales. 





Figura 12.- Escenario probabilista en terminos Figura 13.- Escenario probabilista en 
macrosismicos. terminos espectrales. 



En terminos espectrales el terremoto de Cardedeu se ha complementado con otro 
terremoto ocurrido en 1 .428 en la provincia de Girona con una intensidad epicentral IX (MSK) 
que, pese a producirse a mayor distancia, fue percibido en la ciudad (Olivera et al. 2006). La 
idea es unir las caracteristicas espectrales de un sismo cercano, el de Cardedeu, con las de un 
sismo lejano, como el de 1 .428. Para cada periodo del espectro de respuesta, se ha tornado el 
mayor valor de ambos terremotos para asi definir, de forma conservadora, todo el rango de 
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periodos, con especial cuidado en los periodos largos. La Figura 11 muestra los espectros de 
respuesta zonificados. 

Para el escenario probabilista se ha utilizado el terremoto con una probabilidad del 10% 
en 50 anos que corresponde a una intensidad en el basamento rocoso entre VI y VII (Secanell 
et al. 2004). La Figura 12 muestra el escenario porobabilista en terminos de intensidad 
considerando los efectos de suelo. Las intensidades varian entre VI-VII en los afloramiento 
rocosos y VII-VIII en los terrenos blandos de los deltas y frente de mar. Irizarry (2004) obtiene 
los escenarios en terminos espectrales. La Figura 13 muestra los espectros de respuesta 
zonificados corrrespondientes al escenario probabilista. Tanto los espectros de la Figura 
llcomo los de la Figura 13 son espectros elasticos de respuesta 5% amortiguada, que 
constituyen el punto de partida para estimar los espectros de demanda y los puntos de 
desempeno. 

3.4. Los edificios 

Aproximadamente, el 92% de los edificios de viviendas de Barcelona pueden agruparse 
en dos grandes clases: edificios de mampostena u obra de fabrica de ladrillo y edificios de 
hormigon armado. Las bases de datos disponibles de los edificios de la ciudad han permitido 
caracterizarlos tanto mediante un indice de vulnerabilidad (Lantada, 2006) como mediante 
espectros de capacidad (Bonett, 2003; Bonett et al., 2004; Moreno et al. 2004a, b, Moreno, 
2006). Resumimos a continuation los principales resultados obtenidos por estos autores. 

3.4.1. Indices de vulnerabilidad 

El indice de vulnerabilidad V, de cada edificio se obtiene como suma de un indice 
basico V* , un modificador regional A M R y un modificador de comportamiento A M c . Es decir: 

Vj = Vj + AM r + AM C (16) 

El indice basico es diferente para las diferentes tipologias constructivas; los 
modificadores regionales tienen en cuenta las peculiaridades de la region o periodo de 
construction mientras que los modificadores de comportamiento incorporan otros aspectos, 
como por ejemplo caracteristicas geometricas o de position del edifico que influyen en su 
comportamiento sismico. Lantada (2006) describe, detalladamente, el procedimiento de 
asignacion del indice de vulnerabilidad. El modificador regional se fundamenta en la existencia 
de normativas de diseno sismico en la fecha de construction del edificio, penalizando los 
edificios antiguos construidos sin ninguna conciencia colectiva del peligro sismico y bonificando 
los edificios mas modernos que, a pesar de que frecuentemente no existia ninguna obligation 
de diseno sismorresistente, se beneficiaban de mejoras en las practicas constructivas. Se 
consideran dos modificadores de comportamiento: modificador de edificio y modificador de 
position. El modificador de edificio hace referenda al edificio aislado y cuantifica y puntua 
propiedades como el numero de plantas, la irregularidad en planta y en altura, la longitud de la 
fachada y el estado de conservation. El modificador de position tiene en cuenta otras 
caracteristicas del edificio en relation con el agregado al que pertenece. Asi, por ejemplo, se 
tiene en cuenta la diferencia de altura entre edificios colindantes y la position en referencia al 
conjunto, agregado o manzana a la que pertenece. Concretamente, este modificador penaliza 
los edificios esquina y terminales, aumentado su indice de vulnerabilidad respectivamente en 4 
y en 6 centesimas y bonifica los edificios que se hallan en medio de otros dos edificios 
disminuyendo su indice de vulnerabilidad en 4 centesimas. 

Los indices de vulnerabilidad de los edificios de mampostena toman valores entre 0.7 y 
casi la unidad, con un valor medio de 0.87. Los indices de los edificios de hormigon armado 
son menores, en el rango entre 0.4 y 0.85 y con un valor medio de 0.65. La Figura 14 
muestra, por zonas de estudio (ZRP), los indices medios de los edificios de mamposteria. Se 
observa un patron radial, con los mayores indices de vulnerabilidad en el centra historico de la 
ciudad donde existe una importante concentration de edificios antiguos de una calidad sismica 
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deficiente. La Figura 15 muestra la distribution de los indices de vulnerabilidad medios de las 
zonas de estudio para los edificios de hormigon armado. Se observa una ligera diminution de 
la vulnerabilidad y se pierde el patron radial. Este analisis de vulnerabilidad pone de manifiesto 
la baja calidad sfsmica de los edificios de viviendas de la ciudad, hecho que es tfpico de las 
ciudades situadas en zonas de escasa peligrosidad sfsmica y que, en consecuencia, son 
ajenas a cualquier preocupacion por la protection sismorresistente. 




Vulnerabilidad 
_ 0 . 6-0 7 
□ 0 . 7 - 0.8 
■ 0.8 - 0.9 
0 . 9 - 1.1 
Noediftcios- 
Qtras tipologias 



Figura 14.- Indices de vulnerabilidad 
medios para los edificios de 
mamposteria. 




indice Vulnerabilidad 
I I 0.4 - 0.6 
0 . 6 - 0.7 
0 . 7 - 0.8 
0.8 - 1.0 
No edificios- 
Otras tipologias 



Figura 15.- indices de vulnerabilidad 
medios para los edificios de 
hormigon armado. 



3.4.2. Espectros de capacidad y curvas de fragilidad 

Para aplicar el metodo del espectro de capacidad los edificios se clasifican en altos, 
medianos y bajos y se han usado pianos estructurales y constructivos para su modelado y 
estudio de su comportamiento dinamico mediante programas de analisis estructural no lineal. 
Bonett (2003) obtiene los espectros de capacidad para las tres categorfas de edificios de 
mamposteria usando el programa TreMuri (Galasco et al. 2002). La Figura 16 muestra los 
espectros bilineales de capacidad para los edificios medianos y altos especificando sus 
caracterfsticas de altura. 




Figura 16.- Espectros de capacidad para los 
edificios de mamposteria 
medianos y altos. 




Figura 17.- Espectros de capacidad para los 
edificios de hormigon armado 
bajos, medianos y altos. 
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Los edificios bajos son casas unifamiliares de una o dos plantas con una altura inferior 
a los 6 metros y su espectro bilineal no se muestra en la Figura 16 debido a que su elevada 
resistencia, comparada con la de los edificios medianos y altos, complica su representation a 
la misma escala. Para este edificio, de pequena altura, los puntos de cedencia y de capacidad 
ultima respectivamente son: 

( Sdy,Say ) = (0.27 cm, 0. 65 g) y ( Sdu,Sau ) = (1 36cm, 0.56g) 

Se observa como la capacidad de los edificios de mamposteria disminuye con la altura. 
La Figura 17 muestra las caracteristicas en altura y los espectros de capacidad de los edificios 
de hormigon armado (Moreno et al. 2004a; Moreno, 2006). De nuevo se observa como la 
capacidad disminuye con la altura. A partir de los espectros de capacidad se estiman las 
curvas de fragilidad mediante el procedimiento descrito en el apartado 2.4.2 (ver tambien Tabla 
2 y Tabla 3). La Tabla 4 muestra los parametros de las cun/as de fragilidad para los diferentes 
tipos de edificios considerados. 

Tabla 4.- Parametros que definen las curvas de fragilidad para cada tipo de edificio 



Tipo de edificio 




Parametros d 


le las curvas de fra 


gilidad 




Sd x 


A 


Sd 2 


A 


Sd 3 


A 


Sd 4 


A 


Mamposteri 

a 


Bajo 


0.19 


0.28 


0.27 


0.37 


0.54 


0.54 


1.36 


0.72 


Mediano 


0.44 


0.40 


0.63 


0.50 


1.20 


0.75 


2.91 


0.70 


Alto 


0.46 


0.30 


0.68 


0.65 


1.68 


0.65 


2.61 


0.65 


Hormigon 


Bajo 


0.49 


0.28 


0.70 


0.37 


1.84 


0.82 


5.24 


0.83 


mediano 


0.99 


0.28 


1.42 


0.36 


2.34 


0.50 


5.11 


0.61 


Alto 


1.33 


0.28 


1.89 


0.29 


2.59 


0.34 


4.68 


0.45 



3.5. El dano fi'sico directo: escenarios 

Una vez determinados los indices de vulnerabilidad y los espectros de capacidad los 
edificios han quedado caracterizados para poder aplicar los metodos descritos en el apartado 
2.4 que permiten evaluar las matrices de probabilidad de dano correspondientes a los 
escenarios considerados, tanto en terminos de intensidad macrosismica como en terminos 
espectrales. Estas matrices de probabilidad de dano son las que nos van a permitir desarrollar 
escenarios de dano fisico directo. 

3.5.1. Metodo del indice de vulnerabilidad 

Para cada edificio se conoce su indice de vulnerabilidad y la intensidad sismica 
esperada en el emplazamiento. La ecuacion (15) permite estimar el estado de dano medio o 
mas probable d* que, a su vez, permite estimar las matrices de probabilidad de dano mediante 
las ecuaciones (9), (5), (6) y (7). El indice de vulnerabilidad minimo de 0.4 se obtiene para un 
edificio de hormigon, mientras que el maximo de 1.24 se obtiene para un edificio de 
mamposteria. Las intensidades de los escenarios determinista y probabilista varian entre VI y 
VIII. La Tabla 5 muestra las matrices de probabilidad de dano para un edificio con un indice de 
vulnerabilidad 0.4. La Tabla 6 es analoga a la anterior para un indice de vulnerabilidad de 1.24. 
La Figura 18 muestra el factor de dano medio para las zonas de estudio (ZRP) de la ciudad 
para el escenario determinista. La Figura 19 es analoga a la anterior para el escenario 
probabilista. Para facilitar la comparacion con los escenarios correspondientes al espectro de 
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capacidad el factor central de dano se ha normalizado a 4 en lugar de a 5. 

3.5.2. Metodo del espectro de capacidad 

Cruzando los espectros de respuesta que definen los escenarios y los espectros de 
capacidad que definen los edificios se obtienen los puntos de desempeno. Entrando en las 
curvas de fragilidad con el correspondiente desplazamiento espectral, se estiman las 
probabilidades de cada estado de dano. Bonett (2003) obtuvo los puntos de desempeno y las 
matrices de probabilidad de dano para los edificios de mampostena. La Tabla 7 muestra la 
matriz de probabilidad de dano correspondiente a los edificios de mampostena. 



Tabla 5.- Matriz de probabilidad de dano para un fndice de vulnerabilidad de 0.4. 
D* es el factor central de dano normalizado a 4. 



Intensidad 


d* 


D* 


Probabilidad del estado de dano | 


Nulo 

(0) 


Leve 

(D 


Moderado 

(2) 


Severo 

(3) 


Extensivo 

(4) 


Completo 

(5) 


VI 


0.0899 


0.0719 


0.9680 


0.0282 


0.0035 


0.0003 


0.0000 


0.0000 


VI-VII 


0.1376 


0.1101 


0.9459 


0.0473 


0.0063 


0.0006 


0.0000 


0.0000 


VII 


0.2093 


0.1674 


0.9063 


0.0803 


0.0121 


0.0012 


0.0001 


0.0000 


VII-VIII 


0.3162 


0.2530 


0.8365 


0.1360 


0.0245 


0.0029 


0.0001 


0.0000 


VIII 


0.4721 


0.3777 


0.7199 


0.2212 


0.0510 


0.0074 


0.0005 


0.0000 



Tabla 6.- Matriz de probabilidad de dano para un fndice de vulnerabilidad de 1 .24. 
D* es el factor central de dano normalizado a 4. 



Intensidad 


d* 


D* 


Probabilidad del estado de dano I 


Nulo 

(0) 


Leve 

(D 


Moderado 

(2) 


Severo 

(3) 


Extensivo 

(4) 


Completo 

(5) 


VI 


3.1883 


2.5506 


0.0016 


0.0380 


0.1685 


0.3372 


0.3468 


0.1078 


VI-VII 


3.6553 


2.9242 


0.0002 


0.0104 


0.0747 


0.2366 


0.4115 


0.2666 


VII 


4.0382 


3.2306 


0.0000 


0.0023 


0.0255 


0.1232 


0.3511 


0.4979 


VII-VIII 


4.3320 


3.4656 


0.0000 


0.0004 


0.0072 


0.0497 


0.2179 


0.7248 


VIII 


4.5462 


3.6370 


0.0000 


0.0001 


0.0018 


0.0162 


0.1004 


0.8816 



Se observa como el dano esperado para un terremoto relativamente pequeno es 
relativamente alto. Por ejemplo, para el escenario determinista, aproximadamente un 45% de 
los edificios altos y medianos de mampostena no reforzada localizados en la Zona II sufrirfan 
un estado de dano entre moderado y severo. Para el escenario probabilista los edificios 
danados estarfan alrededor del 75%. Este hecho confirma la alta vulnerabilidad de este tipo de 
edificios. Moreno (2006) efectua un estudio analogo para los edificios de hormigon armado, 
obteniendo las matrices de probabilidad de dano correspondientes. A partir de estas matrices 
de probabilidad de dano, se han construido los escenarios de dano ffsico directo. La Figura 20 
muestra el escenario de dano determinista. La Figura 21 muestra el caso probabilista. 

El detalle y calidad de la base de datos de los edificios de la ciudad nos ha permitido 
realizar el analisis de vulnerabilidad y de dano edificio por edificio. Por lo tanto, el sistema de 
informacion geografica utilizado permite obtener mapas detallados de cualquier distrito, barrio o 
division administrativa a nivel de edificio. Una forma tambien util de representar el dano es 
mediante mapas de probabilidad de ocurrencia de un determinado estado de dano. La Figura 
22 y la Figura 23 muestran respectivamente las probabilidades del estado de dano moderado y 
severo para el escenario determinista y para el distrito mas antiguo de la ciudad ( Ciutat Vella). 
La Figura 24 y la Figura 25 muestran las probabilidades de los mismos estados de dano para 
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el distrito mas emblematico de Barcelona ( L’Eixample ) donde se acumula una importante 
actividad economica y una elevada densidad de poblacion. 
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Dano medio esperado 
0-0.5 No dano 
0.5 -0.7 Leve 
0.7 -0.9 
0.9 -1.1 

1.1 -1.3 
1 3 - 1.5 

1.5 - 1.7 Moderado 

1. 7- 1.9 
1.9 -2.1 

2.1 -2.3 

2.3 - 2.5 

2.5 - 2.7 Severn 

2.7- 29 
2.9 -3.1 

3.1 -3.3 

3.3- 35 

3.5-4 Completo 
^7ZZ\ Sin edificios- Otras tipologfas 




Figura 18.- Distribucion del grado de dano medio esperado por ZRP 
para el escenario determinista. Metodo del indice de 
vulnerabilidad. 
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Dano medio esperado 



0 - 0.5 No dano 

□ 0.5 -0.7 Leve 

I 1 0.7 -0.9 

□ 0.9 - 1.1 
Z3 n-13 

□ 1.3 - 1.5 

□ 1.5 - 1.7 Moderado 
^ 17 - 1.9 

| 1 1.9-2. 1 

I I 2.1 -2.3 

|H] 2 3-2.5 

□ 2.5 -2.7 Severe 
|~1 2.7-29 

_ 2 9-3.1 
H 3.1 -3.3 



3.3-35 



im 3.5 - 4 Completo 

Sin edificios- Otras tipologfas 



Figura 19.- Distribucion del grado de dano medio esperado por ZRP para el 
escenario probabilista. Metodo del indice de vulnerabilidad. 
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Dano medio esperado 

0 - 0.5 No dano 
0.5 - 0.7 Leve 
0.7 -0.9 
0.9-1. 1 

1.1 -1.3 

1 .3- 1.5 

1.5 - 1.7 Moderado 
1.7-1. 9 
1. 9-2.1 

2.1 -2.3 

2.3 - 2.5 
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Figura 20.- Distribution del grado de dano medio esperado por ZRP 
para el escenario determinista. Metodo del espectro de 
capacidad. 
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Dano medio esperado 



0 - 0.5 No dano 
0.5 - 07 Leve 
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I I 1.1 -1-3 

| 1 1. 3-1.5 

□ 1.5 - 1.7 Moderado 
1 17-1.9 



ZD 



1 .9 - 2.1 

2.1 -2.3 

2.3 - 2.5 
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27-2.9 
2.9 -3.1 

3.1 -3.3 

3.3 - 3.5 

3.5 -4.0 Completo 



Figura 21.- Distribution del grado de dano medio esperado por ZRP para el 
escenario probabilista. Metodo del espectro de capacidad. 
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Tabla 7.- 

Matrices de probabilidad de dano para los edificios de mamposterta de 
Barcelona obtenidas mediante el metodo del espectro de capacidad. d*e s el 
factor central de dano considerando 4 estados de dano diferentes del estado de 

dano nulo. 





Escenario determinista 


Escenario probabilista | 






Probabilidades de los estados de dano 


D 


Probi 

estad 


ibilidades de los 
os de dahos 


d* 




Zona 


0 


1 


2 


3 


4 




1 






D 


Edificios 

bajos 


D 


0,950 


0,037 


0,011 


0,002 


0,000 
















D 


0,737 


0,189 


0,063 


0,009 


0,001 




0.287 


0.416 


0.249 


0.042 


0.006 




m 


0,917 


0,061 


0,018 


0,003 


0,001 
















R 


1,000 


0,000 


0,000 


0,000 


0,000 
















Edificios 

medianos 


D 


0,003 


0,166 


0,399 


0,353 


0,079 




0.000 




0.322 








D 


0,121 


0,384 


0,289 


0,189 


0,017 




0.007 


0.160 










m 


0,273 


0,364 


0,215 


0,139 


0,009 








0.395 








R 


0,623 


0,193 


0,105 


0,076 


0,003 














1.600 


Edificios 

altos 


D 


0,003 


0,145 


0,389 


0,371 


0,092 




0.000 


0.019 


0.248 








D 


0,135 


0,388 


0,281 


0,178 


0,018 




0.002 


0.133 


0.385 


0.381 


0.099 


2.441 


m 


0,307 


0,369 


0,195 


0,120 


0,009 




0.014 


0.23 


0.386 








R 


0,647 


0,205 


0,086 


0,059 


0,003 




0.632 


0.257 


0.095 


0.014 


0.002 




Zonas correspondientes a la Microzonificacion stsmica de 1 


a ciud 


ad: 



R - Afloramiento rocoso. I - Suelos blandos. II y III - Suelos intermedios. 
Estados de dano: O-Nulo, 1-Leve, 2-moderado, 3-Severo, 4-Completo. 
d* es el estado de dano medio o factor central de dano. 



4. RESUMEN, DISCUSION Y CONCLUSION 

En este trabajo se han expuesto dos metodologfas de generacion y analisis de 
escenarios de dano sismico entendidos como la prediccion de la situacion esperable en caso 
de ocurrir un terremoto. Ambos metodos se han aplicado a una ciudad situada en una region 
donde el peligro sismico es bajo. 

El metodo del tndice de vulnerabilidad considera la accion en terminos de intensidad 
macrostsmica y caracteriza la vulnerabilidad del edificio mediante un tndice. La prediccion del 
dano esperado se efectua a partir del ajuste de funciones semiemptricas que permiten estimar, 
de forma completa, tanto el dano esperado o estado de dano mas probable, como las 
probabilidades de ocurrencia de cada uno de los estados de dano. El metodo es 
extremadamente versatil para caracterizar el comportamiento sismico de los edificios dado que 
permite modular el tndice de vulnerabilidad representative de la clase, mediante modificadores 
regionales y de edificio. De esta forma un mismo tipo de edificio puede recorrer un amplio 
rango de indices y, por lo tanto, de acuerdo a la escala macrostsmica EMS’98, un amplio rango 
de vulnerabilidades. La aplicacion a grandes nucleos urbanos es rapida y sencilla y permite 
obtener escenarios robustos cuando se interpretan desde la optica probabilista. 

El metodo del espectro de capacidad se basa en el analisis del comportamiento 
mecanico de los edificios. Efectivamente puede entenderse como un metodo avanzado que 
han incorporado nuevas metodologtas de analisis de dano y de riesgo sismico en grandes 
aglomeraciones urbanas. Con todo requiere una mayor cantidad y calidad de informacion, tanto 
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en lo que se refiere a la action sfsmica como en lo que se refiere a la caracterizacion de la 
fragilidad de los edificios. Para lo primero se requieren datos de aceleracion sfsmica; para lo 
segundo se requiere analizar modelos mecanicos de los edificios mediante programas de 
analisis no lineal de estructuras, para lo que se requieren tambien importantes simplificaciones. 
Se trata pues de un metodo potente pero costoso en tiempo de computo y en cantidad y 
calidad de datos. 




No edificio- 
otras tipo logfas 



Figura 22.- Probabilidad del estado de 
dano Moderado en el distrito de 
Ciutat Vella para el escenario 
determinista y el metodo del 
espectro de capacidad. 



Figura 23.- Probabilidad del estado de dano 
Severn en el distrito de Ciutat 
Vella para el escenario 
determinista y el metodo del 
espectro de capacidad. 




| 10-20 
| 20-30 
| 30-40 
140-50 
No edificio- 
otras tipologfa 













3040 
4050 

No edificio- 
otras tipologfas 



Figura 24.- Probabilidad del estado de 
dano moderado en el distrito de 
I’Eixample para el escenario 
probabilista y el metodo del 
espectro de capacidad. 
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40-50 

No edificio- 
otras tipologfas 




Figura 25.- Probabilidad del estado de dano 
Severn en el distrito de I’Eixample 
para el escenario probabilista y el 
metodo del espectro de 
capacidad.. 



La aplicacion de ambos metodos a la ciudad de Barcelona ha proporcionado resultados 
consistentes. El analisis de un escenario determinista y de un escenario probabilista ha 
permitido predecir, con un importante nivel de detalle, la situation que la ocurrencia de estos 





26 



L. G. Pujades, A. H. Barbat y N. Lantada 



terremotos producirfa en la ciudad. A pesar de que el metodo del espectro de capacidad puede 
considerarse mas puntero e ingenieril que el metodo del fndice de vulnerabilidad, su aplicacion 
a la ciudad se basa solo en el modelado y caracterizacion de 6 tipos de edificios. Todos los 
edificios de la ciudad se suponen bien representados por alguna de estas 6 clases 
estructurales. El metodo del fndice de vulnerabilidad permite afinar mas y, de forma mas 
sencilla, permite representar una amplia gama de tipos de edificios asf como discriminar, dentro 
de cada tipologfa, las caracterfsticas geometricas y de posicion dentro de agregados de los 
edificios individuales. Asf pues, consideramos que el metodo del fndice de vulnerabilidad en el 
caso de Barcelona proporciona una imagen mas detallada y realista de los escenarios de dano. 
Los principales motivos de esta afirmacion residen por una parte en la definicion de la accion 
sfsmica y por otra en la versatilidad y facilidad en la caracterizacion de los edificios. Barcelona 
no dispone de ningun dato de aceleracion y sf de intensidad. La caracterizacion de los edificios 
es mas fina en el metodo del mdice de vulnerabilidad que en el metodo del espectro de 
capacidad. Con todo el metodo del espectro de capacidad tambien permite obtener excelentes 
resultados, afina mas en la definicion de la accion sfsmica y en el comportamiento dinamico de 
los edificios. 

A pesar de la dificultad de comparar de forma directa los resultados obtenidos mediante 
ambos metodos, comparacion que se complica al considerar diferentes estados de dano, 
ambos metodos han dado resultados consistentes, poniendo de manifiesto la alta 
vulnerabilidad y fragilidad de los edificios de Barcelona donde, el relativamente bajo peligro 
sfsmico y el rapido crecimiento de la poblacion han producido un entorno urbano altamente 
vulnerable a la accion de los terremotos haciendo predecibles danos importantes para 
terremotos relativamente pequenos. 

El dano ffsico directo es el punto de partida para la consideracion del riesgo que incluye 
danos a la poblacion, a las instituciones y servicios asf como al funcionamiento del sistema 
economico y social que sobrepasa el ambito definido por la division territorial de una ciudad. 
Existen estudios (ver por ejemplo ATC-13 1985; Coburn y Spence, 1992; 2002) que 
proporcionan herramientas para incorporar al analisis del dano otros aspectos como por 
ejemplo el numero de heridos y muertos, el numero de personas sin hogar, el volumen de 
escombros y el costo economico entre otros. Estas herramientas, en general, suelen ser 
formulas empfricas compatibles con datos de terremotos ocurridos en el pasado y suelen 
basarse en el conocimiento de las probabilidades de ocurrencia de los estados de dano ffsico. 
Por lo tanto, tanto los metodos expuestos, como el caso de aplicacion aquf presentado 
constituyen una herramienta importante previa a otros analisis mas sofisticados de riesgo 
sfsmico. 
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RESUMEN 

Se presentan los resultados de un estudio realizado a diferentes disenos de un edificio 
regular de 15 pisos con base en marcos ductiles de concreto reforzado. Dicho estudio tiene 
como proposito detectar las diferencias entre los disenos que contemplan e ignoran la 
disposicion de las secciones 7.3 y 7.3.5. 1 de la propuesta de las Normas Tecnicas 
Complementarias para Diseno y Construccion de Estructuras de Concreto (NTCC-01) del 
Reglamento del Distrito Federal de Mexico relacionada con el incremento en 70 por ciento de la 
fraccion de la carga axial debida a sismo para la obtencion de la carga axial de diseno de 
columnas de marcos ductiles disenados con un factor de comportamiento sismico Q = 4, asi 
como comentar las implicaciones de la aplicacion de esta disposicion. 



SUMMARY 

This paper presents the results of a research study devoted to identify similarities and 
differences between the designs of special moment-resisting frames when sections 7.3 and 
7.3.5. 1 of the draft version for NTCC-OI’s of the Mexico’s Federal District Code is incorporated 
or omitted. These sections request to increment in 70 percent the axial force’s portion related to 
the lateral loads caused by earthquakes for the design of columns in special moment resisting 
frames designed with a response modification factor Q=4. Regular 15-story buildings structured 
with special moment resisting reinforced concrete frames were design for this purpose. 
Nonlinear dynamic analyses of such designs were conducted in order to assess the differences 
and discuss the implications that the application of these sections has on the design process 
and on the expected nonlinear behavior. 
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1 INTRODUCCION 

Las estructuras se disenan para resistir sin colapso, aunque con dano estructural y no 
estructural, a lo que se denomina como sismo de diseno. Esto implica aceptar un 
comportamiento no lineal de importancia y ser racional en el diseno, reduciendo las fuerzas 
laterales con un factor de reduccion de fuerzas por comportamiento inelastico (Q’ en los 
reglamentos de Mexico, en la literatura internacional) para suministrar a la estructura una 
resistencia varias veces menor a la requerida para mantenerla en el intervalo elastico (Figura 
1 )- 



Una de las razones que establecen los reglamentos de diseno sismico de Mexico para 
poder dividir entre Q’ es la reduccion de fuerzas laterales que se justifica por el comportamiento 
no lineal global de la estructura (supuesto de manera simplista como elastoplastico, Figura 1), 
el cual es afectado por la no simultaneidad y no uniformidad en la formacion de las 
articulaciones plasticas. 

Fuerza Sfsmica 




Figura 1 Conceptualization del: a) factor de comportamiento sismico Q y, b) factor de reduccion de las 
fuerzas sismicas con fines de diseno Q', de los reglamentos de diseno sismico de Mexico, 
considerando que, por simplicidad, se supone que el sistema inelastico equivalente tiene un 
comportamiento elastoplastico perfecto. 

En el diseno de estructuras de marcos esbeltos o flexibles existe generalmente un 
aumento de los momentos de volteo y consecuentemente de las cargas axiales en las 
columnas, incremento que en ocasiones no se toma en cuenta en el analisis elastico realizado 
con fuerzas laterales reducidas, y que es particularmente importante en las columnas de los 
primeros entrepisos de edificios altos y esbeltos. 

Este fenomeno ha motivado la propuesta de algunos enfoques para el diseno de 
columnas de marcos ductiles de concreto con los que se pretende cubrir la posibilidad que la 
fuerza sfsmica horizontal exceda a la obtenida dividiendo el coeficiente sismico entre el factor 
de reduccion de las fuerzas sismicas con fines de diseno, Q’, a pesar que estos marcos no 
formen parte de estructuras flexibles o esbeltas. 

Estas propuestas, que como se comenta en Luna (2000) y Luna y Tena (2002), estan 
basadas en estudios metodologicamente inapropiados, han sido impuestas desde 1987 en las 
Normas Tecnicas Complementarias para Diseno y Construction de Estructuras de Concreto 
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(NTCC) del Reglamento para Construcciones del Distrito Federal de Mexico, sin evaluar su 
impacto en los disenos que se obtienen al aplicarlas. Las NTCC son ampliamente adoptadas 
sin mayores variantes por reglamentos estatales y municipales en Mexico. 

De esta manera, en la revision de las NTCC publicada en 1996 (NTCC-96), se 
establecio que al calcular la carga axial de diseno de columnas de marcos ductiles para la cual 
se value el momento resistente M c de una columna, la fraccion de dicha carga debida al sismo 
deberia tomarse igual a 1 .7 veces la calculada a partir del analisis cuando esto conduzca a un 
momento resistente menor. 

Esta disposicion ha tenido un impacto muy negativo, pues complica enormemente el 
proceso de diseno de marcos ductiles. Los disenadores mexicanos han tenido que realizar 
predisenos practicamente de forma manual por no contar con el software adecuado, ya que al 
no existir esta disposicion en ningun otro reglamento vanguardista en el ambito internacional y 
siendo los programas de analisis y diseno mas utilizados en los despachos de ingenieria en 
Mexico aun de procedencia extranjera, no incorporan estas disposiciones. Todo esto, ente otras 
razones, desincentivaba el que se disenara con marcos ductiles de concreto reforzado en 
Mexico (Rioboo, 1995). 

Por esta razon, en la propuesta de NTCC-01, se establecio que el incremento del 
setenta por ciento de la carga axial debida al sismo aplicaria solo para el diseno de estructuras 
con Q=4, lo que libera de esta disposicion sin sustento a los marcos ductiles disenados con 
Q=3. Sin embargo, en opinion de los autores, esta modificacion es claramente insuficiente e 
insatisfactoria a la luz de lo expuesto anteriormente. 

Por lo anterior, se realizo un estudio en diferentes disenos hechos para un edificio 
regular de 15 pisos con base en marcos ductiles de concreto reforzado, disenados con factor de 
comportamiento sismico Q = 4. Dicho estudio tiene como proposito detectar las diferencias 
entre los disenos que contemplan e ignoran la disposicion de las secciones 7.3 y 7.3.5. 1 de las 
NTCC-01 , relacionada con el incremento en 70 por ciento de la fraccion de la carga axial debida 
a sismo para la obtencion de la carga axial de diseno empleada en el dimensionamiento de las 
columnas de marcos ductiles disenados con Q = 4, asi como comentar las implicaciones de la 
aplicacion de esta disposicion. 

Para tal efecto, se disenaron las columnas del edificio en estudio considerando y 
omitiendo esta disposicion y apegandose totalmente al resto de los criterios de la propuesta de 
las NTCC-01 y de las Normas Tecnicas Complementarias para Diseno por Sismo (NTCS). Se 
considero que los modelos de edificio se encuentran ubicados en la subzona lllb del Distrito 
Federal, pues a esta le corresponde el espectro de diseno con mayor coeficiente sismico, y en 
esta zona esta ubicada la Secretaria de Comunicaciones y Transportes (SCT), lugar donde se 
obtuvo lo que es hasta la fecha el acelerograma mas intenso registrado en terreno blando en 
Mexico y en el mundo, del cual se tomo la componente E-W del acelerograma registrado 
durante el sismo del 19 de septiembre de 1985, con el fin de tener un registro representative y 
de caracteristicas similares al sismo maximo contemplado por el Reglamento del Distrito 
Federal, del cual se pretenden evaluar sus criterios. 

Para efectuar los analisis no lineales se utilizo el programa DRAIN-2DX (Prakash et al. 
1992), utilizando modelos que consideran la accion del diafragma rigido y por tanto distribuyen 
las fuerzas sismicas entre los marcos que resisten en una misma direccion. Se consideraron los 
efectos P-A y no se tomo en cuenta la interaccion suelo-estructura. El estudio se hizo 
considerando resistencias nominales y fuentes adicionales de resistencia (sobrerresistencia). 
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Cabe senalar que, recientemente, los miembros del Comite de las Normas de Concreto 
decidieron modificar esta disposition y establecer ahora en las NTCC-04 vigentes que, en 
marcos ductiles disenados con Q=4 con relation altura-base mayor que 2.0, y que se 
encuentren en la zona III, al dimensionar por flexocompresion, se incremente 50 por ciento la 
fuerza axial y el momento flexionante debidos al sismo. Aunque con este cambio no se tienen 
tantas complicaciones en el diseno como con el criterio anterior, los autores la consideran una 
medida arbitraria, dado que la misma no esta justificada en estudios especificos previos a su 
aprobacion. 



2 EDIFICIO EN ESTUDIO Y METODOLOGIA DE DISENO 
2.1 Description de la estructura y modelos de analisis 

Se trata de un edificio de oficinas regular en planta y en elevacion (figura 2). La 
estructu radon del edificio es totalmente simetrica, con cuatro crujias en la direccion x y tres en 
la direccion y con claros de 8 y 7 m respectivamente. En la direccion x se dispusieron ademas 
vigas secundarias. El edificio consta de 15 niveles con alturas tipicas de entrepiso de 3.4 m, lo 
que da una altura total de 51 m (figura 2). 

La aplicacion u omision del incremento en 70 por ciento de la fraction de la carga axial 
debida a sismo para la determination de la carga axial de diseno, empleada en el 
dimensionamiento de las columnas establece la diferencia entre los modelos de estudio. El 
modelo denominado DQ4CP1.7 es el que incorpora el 70% de incremento en la carga axial, 
mientras que el modelo DQ4CP ignora este incremento. Q4 y CP denotan que el factor de 
comportamiento sismico es cuatro y que el analisis sismico se hizo conforme a la propuesta de 
diseno tradicional que establece el cuerpo principal de la propuesta de NTCS-01. 
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a) Planta y elevacion tipo 

Figura 2 Edificio en estudio. 




b) Modelo de analisis 



3 METODOLOGIA DE DISENO 

El calculo estructural se realizo con el programa de analisis y diseno de edificios de 
concreto ECOgcW version 2.06 (Corona, 2003), el cual aplica las disposiciones del Reglamento 
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para las Construcciones del Distrito Federal (RCDF) y sus Normas Tecnicas Complementarias 
para Estructuras de Concreto (NTCC) para el diseno de miembros a flexion y flexocompresion 
biaxial. Con este software se determinaron las fuerzas sismicas laterales mediante un analisis 
dinamico modal espectral y se incluyeron los efectos P-A. En la figura 2 se presenta el modelo 
tridimensional del edificio en estudio en el programa ECOgcW. 



3.1 Consideraciones para el analisis estructural 

Con la finalidad de obtener estimaciones razonables de los periodos fundamentales de 
cada modelo en su condicion inicial (elastica) y ultima (inelastica) y consecuentemente tener 
una buena estimacion de las demandas sismicas incorporando tanto un sismo frecuente (de 
servicio) como uno intenso (de diseno), se consideraron dos modelos de analisis para cada uno 
de los edificios tipo estudiados: (1) uno en el cual se modelo a todos los miembros estructurales 
con propiedades de secciones brutas, tomando en cuenta que todo edificio nuevo que aun no 
experimenta sismo alguno se encuentra en tales condiciones antes de afrontar a su primer 
sismo (sea intenso o no) y, (2) otro en el cual se modelo de manera muy cruda y aproximada el 
agrietamiento de las vigas, tomando el momento de inercia de las mismas igual a la mitad del 
momento de inercia de la seccion bruta, tal y como se establece en el articulo 1.3.1 de la 
propuesta de NTCC-01. Tal modelacion considera ya un degradacion importante de la rigidez 
de las vigas, que para una estructura disenada con reglamentos modernos solo puede estar 
asociada a un estado ultimo de deformacion, por lo que con tal modelacion se obtiene una 
rigidez secante equivalente que lidera a obtener el periodo esperado al estado ultimo de 
deformacion esperado por el Reglamento (8 U , figura 1). Esta modelacion propuesta por las 
NTCC y otros reglamentos internacionales como el ACI es imperfecta (aunque endiosada por 
algunos colegas), pues entre otras cosas, es bien conocido que la degradacion de rigidez de las 
vigas de los marcos en elevacion y en planta no es uniforme, como se ilustrara secciones mas 
adelante (seccion 5.2, fig 27). Cabe senalar que una perdida del 50% de la rigidez de las vigas 
esta asociada a rotaciones plasticas y agrietamientos importantes, observables a simple vista. 

En este trabajo, en los modelos de analisis estructural se consideraron que las 
secciones de las vigas son rectangulares; es decir, no se tomo en cuenta la participation de la 
losa, lo cual es representative de lo que se hace en muchos despachos de calculo, segun se ha 
obtenido en varias encuestas realizadas en despachos de diseno estructural de Mexico 
(Fuentes 2000, Correa 2005). 

Por fines de absoluta congruencia con las hipotesis hechas en las NTCS, se reviso el 
cumplimiento de las distorsiones de entrepiso permisibles segun lo establecen las NTCS 
(multiplicar por Q las distorsiones obtenidas a partir de un analisis elastico con las fuerzas 
sismicas ya reducidas por Q') con los modelos que consideran secciones brutas de las vigas. 
Tambien se determinaron para estos modelos con secciones brutas las demandas preliminares 
de refuerzo para los miembros estructurales, realizando una revision final empleando para cada 
uno el modelo que tomaba en cuenta el agrietamiento de las vigas. Finalmente, con los 
modelos que consideran el agrietamiento de las vigas (modelo secante equivalente) tambien se 
reviso que las distorsiones obtenidas, sin multiplicar por Q, no sobrepasaran los limites de 
distorsion permisible. Los suministros de refuerzo quedaron asi definidos de tal suerte que se 
satisficieran simultaneamente las demandas de ambos modelos, considerando el agrietamiento 
de las vigas y sin considerar el agrietamiento de las vigas. 

Es importante mencionar que otro de los efectos del agrietamiento de las estructuras es 
el aumento del amortiguamiento y los espectros para diseno sfsmico ofrecidos por las NTCS 
estan definidos para un valor fijo de 5% del amortiguamiento critico, por lo que no se puede 
modelar bien este efecto. 
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Con estas consideraciones el periodo fundamental de la estructura considerando las 
secciones brutas fue T e =1.61s mientras que al considerar el agrietamiento aumento a 
T agr =2.15s (direccion x). En la direccion y se tuvieron los periodos T e =1.17s y T agr =1.48s 
considerando secciones brutas y agrietadas respectivamente. 



3.2 Consideraciones para el diseno estructural 

Se supuso que los elementos no estructurales estan perfectamente desligados de la 
estructura, por lo que se adopto como distorsion de entrepiso h'mite el valor: 

yQ = 0.012 (1) 

donde y es la distorsion de entrepiso producida por las acciones sfsmicas calculadas con las 
ordenadas espectrales reducidas por Q' y Q es el factor de comportamiento sismico (figura 1 ). 

El criterio para determinar el tamano de las secciones transversales de columnas y 
vigas fue satisfacer casi al h'mite, pero sin salirse de lo practico y entrar en lo inconstruible, las 
distorsiones de entrepiso permisibles en el cuerpo principal de las NTCS-01, a la vez que las 
secciones obtenidas admitian el refuerzo necesario dentro de los h'mites establecidos en la 
propuesta de NTCC-01 para tener la resistencia requerida. Esto con la finalidad que los 
resultados del estudio correspondan a una situacion h'mite o cri'tica. De esta manera, la seccion 
transversal de las columnas se vario cada cinco niveles. 

En un principio, como el edificio tiene forma rectangular, se penso en utilizar columnas 
rectangulares orientadas de tal manera que su momento de inercia mayor coincidiera con el 
lado corto del edificio para dar mayor rigidez global al edificio en esta direccion; sin embargo, 
como se noto posteriormente que las vigas en la direccion corta demandaban de grandes 
secciones, fue posible utilizar columnas de seccion transversal cuadrada, ya que al ser un 
edificio estructurado con marcos momento-resistentes, la rigidez lateral de la estructura en la 
direccion corta aumento notablemente con el aumento de las secciones transversales de las 
trabes en esta direccion. 




Distorsion en direccion y 

Lfmites del cuerpo principal 

de las NTCS-01 

Figura 3 Distorsiones de entrepiso de diseno (modelo con secciones brutas) 



El hecho que las vigas paralelas a la direccion y resultaran con una seccion importante, 
le confirio al edificio una rigidez lateral mayor en este sentido que en el largo, lo cual se 
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evidencia con las distorsiones de entrepiso que sufre la estructura con las secciones elegidas 
(figura 3). 

Las trabes paralelas al eje y de todos los modelos tuvieron problemas con resistencia a 
fuerza cortante, ya que al recibir la descarga de las vigas secundarias, estan sometidas a 
cortantes elevados. Este problema fue mas frecuente en los primeros entrepisos, lugar donde el 
cortante debido a sismo es mayor. 

Las dimensiones de las vigas paralelas el eje x fueron establecidas para dar rigidez 
lateral en ese sentido a la estructura. Para su dimensionamiento no rigieron los criterios de 
resistencia, sino que se busco tener secciones tales que no fuera necesario incrementar las 
dimensiones de las columnas para cumplir con las distorsiones de entrepiso permisibles por el 
RCDF. Como se menciono arriba, aunque la direccion x corresponde al lado largo del edificio, 
resulto ser finalmente la direccion mas flexible, por lo que si se hubieran adoptado secciones 
menores de trabes en esta direccion, aunque hubieran admitido el refuerzo longitudinal y 
transversal dentro de los limites permisibles teniendo la resistencia requerida, hubiera sido 
necesario incrementar la seccion de las columnas en este sentido para cumplir con las 
distorsiones de entrepiso permisibles. 

Para el diseno y detallado de los elementos estructurales de los modelos se aplicaron 
las disposiciones para marcos ductiles de la propuesta de NTCC-01 . En el diseno de trabes por 
fuerza cortante, asi como para el diseno de las columnas por flexocompresion y cortante, se 
empleo siempre el metodo riguroso de la propuesta de NTCC-01 que se basa en los principios 
de diseno por capacidad. Para la revision de la resistencia a flexion minima de las columnas de 
la seccion 7.3.2 de la propuesta de NTCC-01 se tomo la carga axial correspondiente a la 
combinacion que rigio para el diseno. Es importante mencionar que la propuesta de NTCC-01 
no indica que carga axial debe adoptarse para realizar esta revision. 



4 DISENO ESTRUCTURAL 

4.1 Descripcion de las secciones de diseno 

Las columnas de los modelos DQ4CP y DQ4CP1.7 presentan tres secciones 
transversales diferentes; 100 x 100 cm (columnas Cl) de los entrepisos uno a cinco, 80 x 80 cm 
(columnas C2) de los entrepisos seis a diez y 65 x 65 cm (columnas C3) de los entrepisos once 
a quince (tablas 2 y 4). El refuerzo longitudinal se encuentra uniformemente distribuido en todos 
los casos (figuras 4 y 5). 

Como para la adopcion de las secciones de las columnas de los modelos DQ4CP y 
DQ4CP1.7 rigio el cumplimiento de la distorsion de entrepiso permisible, no fue necesario 
incrementar la seccion de las columnas del modelo DQ4CP1.7 con respecto a las del DQ4CP, 
limitandose asi las diferencias entre las columnas de uno y otro modelo a la cuantia de refuerzo 
longitudinal. El factor 1.7 interviene solo en el diseno de las columnas, razon por la cual las 
vigas de los modelos DQ4CP y DQ4CP1.7 resultaron identicas. 

Las diferencias entre los disenos DQ4CP y DQ4CP1.7 se presentan fundamentalmente 
en las columnas ubicadas en los primeros entrepisos y en la periferia del edificio. Las columnas 
del modelo DQ4CP1 .7 situadas en los lugares comentados demandan de mayores cantidades 
de acero. En las tablas 1 y 3 se han sombreado las celdas correspondientes a las columnas 
que resultaron con armados diferentes. En las tablas 2 y 4 se observa la diferencia cuantitativa 
en cuanto al refuerzo entre estas columnas. 
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Si bien algunas otras columnas de los primeros entrepisos de los modelos DQ4CP y 
DQ4CP1.7 no sombreadas en las tablas de tipificacion de columnas presentaron diferencias en 
las demandas de refuerzo, estas eran tan pequenas que se desvanecfan al disenarlas 
sujetandose a los diametros de varilla comerciales. 

Las ultimas columnas del primer grupo (columnas Cl) de ambos modelos, que se 
encuentran del entrepiso tres al cinco, estan reforzadas por flexocompresion con la cuantfa 
minima, por lo cual no se observan diferencias entre el modelo DQ4CP y DQ4CP1.7. Las 
columnas centrales permanecen sin cambio en estos modelos, aun con cuantfas de refuerzo 
longitudinal superiores a la minima. 

Las columnas de los niveles superiores (columnas C2 y C3) de los modelos DQ4CP y 
DQ4CP7.3 no presentan diferencias en las demandas de refuerzo longitudinal, a pesar de estar 
reforzadas por flexocompresion con cuantfas superiores a la minima. 



Tabla 1 Tipificacion de las columnas del modelo DQ4CP 



Nivel 


Columnas 


Al y D1 


A2y D2 


A3 y D3 


A4y D4 


A5y D5 


B1 y Cl 


B2yC2 


B3y C3 


B4y C4 


B5yC5 


1 


C1-R1 


C1-R2 


C1-R2 


C1-R2 


C1-R1 


C1-R3 


C1-R2 


C1-R2 


C1-R2 


C1-R3 


2 


C1-R1 


C1-R4 


C1-R4 


C1-R4 


C1-R1 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


3 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


5 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


6 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


7 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


8 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


9 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


10 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


11 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R3 


C3-R2 


C3-R2 


C3-R2 


C3-R3 


12 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R3 


C3-R2 


C3-R2 


C3-R2 


C3-R3 


13 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R3 


C3-R2 


C3-R2 


C3-R2 


C3-R3 


14 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


15 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 



Tabla 2 Refuerzo de las columnas del modelo DQ4CP 



Seccion 


Refuerzo para flexocompresion 


Refuerzo para fuerza cortante 
(en ramas de estribo) 


Tipo 


Dimensiones 

(cm) 


Barras 


Cuantfa (p) 


en los extremos 


en el centra 


C1-R1 


100 x 100 


48 #8 


0.0243 


8 #4 @ 10 cm 


8 # 4 @ 30 cm 


C1-R2 


100 x 100 


34 #8 


0.0172 


8 #4 @ 10 cm 


8 # 4 @ 30 cm 


C1-R3 


100 x 100 


36 #8 


0.0183 


8 #4 @ 10 cm 


8 # 4 @ 30 cm 


C1-R4 


100 x 100 


26 #8 


0.0132 


6 # 4 @ 1 0 cm 


6 # 4 @ 20 cm 


C1-R5 


100 x 100 


32 #8 


0.0162 


6 #4 @ 10 cm 


6 # 4 @ 20 cm 


C2-R1 


80x80 


16 # 8 


0.0127 


5 #4 @ 10 cm 


5 # 4 @ 25 cm 


C2-R2 


80x80 


20 #8 


0.0158 


6 # 4 @ 1 0 cm 


6 # 4 @ 30 cm 


C2-R3 


80x80 


32 #8 


0.0254 


6 # 4 @ 1 0 cm 


6 # 4 @ 30 cm 


C3-R1 


65x65 


12 # 8 


0.0144 


4 # 4 @ 10 cm 


4 # 4 @ 25 cm 


C3-R2 


65x65 


16 # 8 


0.0192 


5 #4 @ 10 cm 


5 # 4 @ 30 cm 


C3-R3 


65x65 


24 #8 


0.0288 


4 # 4 @ 1 0 cm 


4 # 4 @ 25 cm 



Tabla 3 Tipificacion de las columnas del modelo DQ4CP1.7 



Nivel 


Columnas 


Al y D1 | A2 y D2 | A3 y D3 | A4 y D4 | A5 y D5 | B1 y Cl | B2 y C2 | B3 y C3 | B4 y C4 | B5 y C5 
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C1-R1 


C1-R2 


C1-R2 


C1-R2 


C1-R1 


C1-R1 


CI-R3 


CI-R3 


CI-R3 


C1-R1 




C1-R1 


C1-R4 


C1-R4 


C1-R4 


C1-R1 


C1-R1 


C1-R4 


C1-R4 


C1-R4 


C1-R1 




C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 




C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 




C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R4 


C1-R5 


C1-R4 


C1-R4 


C1-R4 


C1-R5 




C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 




C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 




C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 




C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


0 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R1 


C2-R3 


C2-R2 


C2-R2 


C2-R2 


C2-R3 


1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R3 


C3-R2 


C3-R2 


C3-R2 


C3-R3 


2 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R3 


C3-R2 


C3-R2 


C3-R2 


C3-R3 


3 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R3 


C3-R2 


C3-R2 


C3-R2 


C3-R3 


4 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


5 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 


C3-R1 



Tabla 4 Refuerzo de las columnas del modelo DQ4CP1.7 



Seccion 


Refuerzo para flexocompresion 


Refuerzo para fuerza cortante 
(en ramas de estribo) 


Tipo 


Dimensiones 

(cm) 


Barras 


cuantia (p) 


en los extremos 


en el centra 


C1-R1 


100 x 100 


48 #8 


0.0243 


8#4@ 10 cm 


8 # 4 @ 30 cm 


C1-R2 


100 x 100 


36 #8 


0.0183 


8#4@ 10 cm 


8 # 4 @ 30 cm 


C1-R3 


100 x 100 


34 #8 


0.0172 


8#4@ 10 cm 


8 # 4 @ 30 cm 


C1-R4 


100 x 100 


26 #8 


0.0132 


6 # 4 @ 1 0 cm 


6 # 4 @ 20 cm 


C1-R5 


100 x 100 


32 #8 


0.0162 


6 # 4 @ 10 cm 


6 # 4 @ 20 cm 


C2-R1 


80x80 


16 # 8 


0.0127 


5 #4 @ 10 cm 


5 # 4 @ 25 cm 


C2-R2 


80x80 


20 #8 


0.0158 


6 # 4 @ 1 0 cm 


6 # 4 @ 30 cm 


C2-R3 


80x80 


32 #8 


0.0254 


6 # 4 @ 1 0 cm 


6 # 4 @ 30 cm 


C3-R1 


65x65 


12 # 8 


0.0144 


4 # 4 @ 10 cm 


4 # 4 @ 25 cm 


C3-R2 


65x65 


16 # 8 


0.0192 


5 #4 @ 10 cm 


5 # 4 @ 30 cm 


C3-R3 


65x65 


24 #8 


0.0288 


4 # 4 @ 10 cm 


4 # 4 @ 25 cm 



Algunas columnas que bajo carga gravitacional estan sometidas a cargas axiales 
pequenas por estar en la periferia del edificio y corresponderles un area tributaria reducida, por 
sismo podrian estar sometidas a cargas axiales altas, inclusive de tension. Mediante las 
combinaciones de carga podria compensarse y llevar a tener una carga axial de diseno baja. 
Esto, aunado a que muchas de estas columnas ubicadas en los entrepisos superiores resultan 
tener el refuerzo minimo, provoca que tengan una pobre resistencia a flexion. Fue asi que al 
realizar la revision de la resistencia minima a flexion de las columnas conforme al articulo 
7.3.2. 1 de la propuesta de NTCC-01, se detecto que las columnas B1, Cl, B5 y C5 de los 
niveles seis en adelante, no tenian la resistencia minima a flexion requerida, por lo que fue 
necesario incrementar la cuantia de refuerzo longitudinal. 

Aunque las columnas B2, C2, B4 y C4 no se encuentran en la periferia del edificio y, por 
ende, tienen areas tributarias grandes, se conectan a ellas vigas por todas direcciones, lo cual 
les demanda una resistencia a flexion elevada para el cumplimiento del requisito de resistencia 
a flexion minima para las columnas, razon por la cual tambien fue necesario incrementar la 
cuantia de refuerzo longitudinal de estas columnas para cumplir con este requisito. 
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Cl-Rl C1-R2 C1-R3 C1-R4 C1-R5 





C2-R1 C2-R2 C2-R3 C3-R1 C3-R2 C3-R3 

Figura 4 Secciones y refuerzo de las columnas del modelo DQ4CP 



Cl-Rl 



C1-R2 



C1-R3 



C1-R4 



C1-R5 



C2-R1 C2-R2 C2-R3 C3-R1 C3-R2 C3-R3 

Figura 5 Secciones y refuerzo de las columnas del modelo DQ4CP1 .7 



Como se menciono anteriormente, las trabes orientadas en la direccion y resultaron con 
una seccion grande y, por ende, una resistencia importante a flexion, siendo precisamente 
estas trabes las que demandaron de una mayor resistencia a flexion de todas las columnas y 
que provocaron en un principio el incumplimiento de las columnas arriba mencionadas. 



Aunque las columnas B1, Cl, B5 y C5 comprendidas entre los niveles tres y cinco 
habian resultado con la cuantia minima y cumplian con la resistencia minima a flexion 
estipulada en la propuesta de NTCC-01, se aumento su cuantia de tal manera que tuvieran por 
lo menos el mismo numero de varillas de las columnas que se encuentran encima de ellas, las 
cuales sufrieron un incremento en el refuerzo longitudinal por no cumplir con la resistencia a 
flexion minima estipulada en el Reglamento. 



Las columnas de los primeros entrepisos cumplieron holgadamente con el requisito de 
resistencia a flexion minima debido a que tienden a ser mas robustas para cumplir con las 
distorsiones de entrepiso permisibles, ademas que presentan porcentajes de acero superiores 
al minimo por los esfuerzos a que estan sometidas por su ubicacion en elevacion, todo lo cual 
redunda en un momento flexionante resistente mayor pues ademas soportan cargas axiales de 
consideracion. 






Evaluacion del criterio que incrementa en 70% la carga axial por sismo en columnas de... 
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Los requisitos para el confinamiento del nucleo del concreto rigieron la separacion del 
refuerzo transversal en los extremos de las columnas. De esta manera no hubo diferencias 
entre los modelos DQ4CP y DQ4CP1 .7. 

Las cuantia de refuerzo longitudinal para momento negativo de las trabes en direccion x 
(tabla 5) fluctua entre p=0.0048 y p=0.0124 (trabes T1, tabla 6). Como se menciono 
anteriormente, para la definicion de las secciones transversales de estas trabes rigieron los 
requisitos de deformacion, lo cual se refleja en el hecho que unicamente los criterios de 
resistencia a fuerza cortante rigieron en la separacion de los estribos en los extremos de las 
trabes ubicadas entre los niveles uno al diez, mientras que los requisitos de confinamiento del 
nucleo del concreto rigieron la separacion de estribos en los extremos a partir del nivel once. No 
se vario la dimension de la seccion transversal de estas trabes en elevacion por lo que todas 
presentan una seccion de 40 x 80 cm. 



Tabla 5 Tipificacion de las vigas paralelas al eje xde los modelos DQ4CP y DQ4CP1.7 



Nivel 


Marcos A y D, cruji'a: 


Marcos B y C cruji'a: 


1-2 


2-3 


3-4 


4-5 


1-2 


2-3 


3-4 


4-5 


1 


T1-R1 


T1-R1 


T1-R1 


T1-R1 


T1-R2 


T1-R2 


T1-R2 


T1-R2 


2 


T1-R1 


T1-R1 


T1-R1 


T1-R1 


T1-R2 


T1-R2 


T1-R2 


T1-R2 


3 


T1-R1 


T1-R1 


T1-R1 


T1-R1 


T1-R2 


T1-R2 


T1-R2 


T1-R2 


4 


T1-R1 


T1-R1 


T1-R1 


T1-R1 


T1-R2 


T1-R2 


T1-R2 


T1-R2 


5 


T1-R1 


T1-R1 


T1-R1 


T1-R1 


T1-R2 


T1-R2 


T1-R2 


T1-R2 


6 


T1-R1 


T1-R1 


T1-R1 


T1-R1 


T1-R2 


T1-R2 


T1-R2 


T1-R2 


7 


T1-R3 


T1-R3 


T1-R3 


T1-R3 


T1-R4 


T1-R4 


T1-R4 


T1-R4 


8 


T1-R3 


T1-R3 


T1-R3 


T1-R3 


T1-R4 


T1-R4 


T1-R4 


T1-R4 


9 


T1-R5 


T1-R5 


T1-R5 


T1-R5 


T1-R6 


T1-R6 


T1-R6 


T1-R6 


10 


T1-R5 


T1-R5 


T1-R5 


T1-R5 


T1-R6 


T1-R6 


T1-R6 


T1-R6 


11 


T1-R7 


T1-R7 


T1-R7 


T1-R7 


T1-R8 


T1-R8 


T1-R8 


T1-R8 


12 


T1-R7 


T1-R7 


T1-R7 


T1-R7 


T1-R8 


T1-R8 


T1-R8 


T1-R8 


13 


T1-R9 


T1-R9 


T1-R9 


T1-R9 


T1-R10 


T1-R10 


T1-R10 


T1-R10 


14 


T1-R9 


T1-R9 


T1-R9 


T1-R9 


T1-R10 


T1-R10 


T1-R10 


T1-R10 


15 


T1-R9 


T1-R9 


T1-R9 


T1-R9 


T1-R10 


T1-R10 


T1-R10 


T1-R10 



La seccion transversal de las vigas paralelas al eje y (tabla 7) se vario en elevacion con 
el fin de tener las dimensiones que fueran absolutamente necesarias para resistencia a corte y 
una sobrerresistencia a flexion lo mas baja posible en las mismas. 
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Tabla 6 Refuerzo de las vigas paralelas al eje xde los modelos DQ4CP y DQ4CP1.7 



Seccion 


Refuerzo para flexion 


Refuerzo para fuerza cortante 
(en ramas de estribo) 


Tipo 


Dimensiones 

(cm) 


Lecho sup. 


Lecho inf. 


en los extremos 


en el centra 


T1-R1 


40x80 


5 # 10 


3#10 + 2#8 


4 # 2.5 @ 10 cm 


4 # 2.5 @ 12.5 cm 


T1-R2 


40x80 


5 # 10 


6 # 8 


4 # 2.5 @ 10 cm 


4 # 2.5 @ 12.5 cm 


T1-R3 


40x80 


2#10 + 4#8 


1 # 10 + 4 # 8 


4 # 2.5 @ 10 cm 


4 # 2.5 @ 15 cm 


T1-R4 


40x80 


2#10 + 4#8 


5 # 8 


4 # 2.5 @ 10 cm 


4 # 2.5 @ 12.5 cm 


T1-R5 


40x80 


6 # 8 


4 # 8 + 1 #6 


4 # 2.5 @ 15 cm 


4 # 2.5 @ 15 cm 


T1-R6 


40x80 


6 # 8 


4 # 8 


4 # 2.5 @ 15 cm 


4 # 2.5 @ 15 cm 


T1-R7 


40x80 


4 # 8 + 1 #6 


6 # 6 


4 # 2.5 @ 15 cm 


4 #2.5 @20 cm 


T1-R8 


40x80 


4 # 8 + 1 #6 


4 # 6 + 1 #5 


4 # 2.5 @ 12.5 cm 


4 #2.5 @20 cm 


T1-R9 


40x80 


4 # 6 + 2 # 5 


2 # 6 + 2 # 5 


4 # 2.5 @ 12.5 cm 


4 #2.5 @35 cm 


T1-R10 


40x80 


4 # 6 + 2 # 5 


4 # 5 


4 # 2.5 @ 12.5 cm 


4 #2.5 @35 cm 



Tabla 7 Tipificacion de las vigas paralelas al eje yde los modelos DQ4CP y DQ4C P1 .7 



Nivel 


Marcos 1 y 5, crujia: 


j Marcos 2 y 4, crujia: j 


Marco 3, crujia: | 


A-B 


B-C 


C-D 


A-B 


B-C 


C-D 


A-B 


B-C 


C-D 


1 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


2 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


3 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


4 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


5 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


6 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


7 


T1-R1 


T1-R2 


T1-R1 


T1-R3 


T1-R4 


T1-R3 


T1-R5 


T1-R6 


T1-R5 


8 


T2-R1 


T2-R2 


T2-R1 


T2-R3 


T2-R4 


T2-R3 


T2-R5 


T2-R6 


T2-R5 


9 


T2-R1 


T2-R2 


T2-R1 


T2-R3 


T2-R4 


T2-R3 


T2-R5 


T2-R6 


T2-R5 


10 


T2-R1 


T2-R2 


T2-R1 


T2-R3 


T2-R4 


T2-R3 


T2-R5 


T2-R6 


T2-R5 


11 


T2-R7 


T2-R8 


T2-R7 


T2-R9 


T2-R10 


T2-R9 


T2-R1 1 


T2-R12 


T2-R1 1 


12 


T2-R7 


T2-R8 


T2-R7 


T2-R9 


T2-R10 


T2-R9 


T2-R1 1 


T2-R12 


T2-R1 1 


13 


T3-R1 


T3-R2 


T3-R1 


T3-R3 


T3-R4 


T3-R3 


T3-R5 


T3-R6 


T3-R5 


14 


T3-R1 


T3-R2 


T3-R1 


T3-R3 


T3-R4 


T3-R3 


T3-R5 


T3-R6 


T3-R5 


15 


T3-R1 


T3-R2 


T3-R1 


T3-R3 


T3-R4 


T3-R3 


T3-R5 


T3-R6 


T3-R5 



Se tuvieron tres secciones distintas, 65 x 100 cm (trabes T1) de los niveles uno al siete, 
60 x 90 cm (trabes T2) de los niveles ocho a doce y 45 x 80 cm (trabes T3) de los niveles trece 
a quince. La separacion de los estribos en los extremos es baja (tabla 8), rigiendo para la 
misma los requisitos de resistencia a fuerza cortante de los niveles uno a once, mientras que de 
los niveles doce a quince rigio el requisito de confinamiento del nucleo. Las cuantias de 
refuerzo longitudinal para momento negativo fluctuan entre p=0.0057 y p=0.0072 (trabes T1), 
p=0.0044 y p=0.0075 (trabes T2) y p=0.0052 y p=0.0060 (trabes T3). 







Evaluation del criterio que incrementa en 70% la carga axial por sismo en columnas de... 
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T abla 8 Refuerzo de las vigas paralelas al eje y de los modelos DQ4CP y DQ4CP1 7 



Section 


Refuerzo para flexion 


Refuerzo para fuerza cortante 
(en ramas de estribo) 


Tipo 


Dimensiones 

(cm) 


Lecho sup. 


Lecho inf. 


en los extremos 


en el centra 


T1-R1 


65 x 100 


4#10 + 3#8 


1 # 10 + 6 # 8 


4 # 3 @ 1 0 cm 


4 # 3 @ 10 cm 


T1-R2 


65 x 100 


4#10 + 3#8 


1 # 10 + 6 # 8 


4#3@ 10 cm 


4 # 3 @ 35 cm 


T1-R3 


65 x 100 


3#10 + 4#8 


6 # 8 


4 # 3 @ 8 cm 


4 # 3 @ 20 cm 


T1-R4 


65 x 100 


2#10 + 5#8 


6 # 8 


4 # 3 @ 10 cm 


4 # 3 @ 30 cm 


T1-R5 


65 x 100 


4 # 10 + 1 #8 


5 # 8 


4 # 3 @ 9 cm 


4 # 3 @ 20 cm 


T1-R6 


65 x 100 


4 # 10 + 1 #8 


5 # 8 


4#3@ 10 cm 


4 # 3 @ 40 cm 


T2-R1 


60x90 


8 # 8 


6 # 8 


4#3@ 10 cm 


4 # 3 @ 12.5 cm 


T2-R2 


60x90 


8 # 8 


6 # 8 + 1 #6 


4 # 3 @ 10 cm 


4 # 3 @ 25 cm 


T2-R3 


60x90 


5# 10 


5 # 8 


4 # 3 @ 9 cm 


4 # 3 @ 20 cm 


T2-R4 


60x90 


4 # 10 + 1 #8 


5 # 8 


4#3@ 10 cm 


4 # 3 @ 30 cm 


T2-R5 


60x90 


2#10 + 4#8 


3 # 8 + 2 # 6 


4 # 3 @ 9 cm 


4 # 3 @ 20 cm 


T2-R6 


60x90 


2#10 + 3#8 


3 # 8 + 2 # 6 


4 # 3 @ 10 cm 


4 # 3 @ 40 cm 


T2-R7 


60x90 


6 # 8 


3 # 8 + 2 # 6 


4#3@ 12.5 cm 


4 # 3 @ 40 cm 


T2-R8 


60x90 


4 # 8 + 1 #10 


3 # 8 + 2 # 6 


4 # 3 @ 15 cm 


4 # 3 @ 40 cm 


T2-R9 


60x90 


2#10 + 3#8 


2 # 8 + 2 # 6 


4 # 3 @ 1 0 cm 


4 # 3 @ 25 cm 


T2-R10 


60x90 


1 # 10 + 4 # 8 


2 # 8 + 2 # 6 


4 # 3 @ 10 cm 


4 # 3 @ 40 cm 


T2-R1 1 


60x90 


1 # 10 + 4 # 8 


4 # 6 + 1 #5 


4 # 3 @ 1 0 cm 


4 # 3 @ 30 cm 


T2-R12 


60x90 


4 # 8 + 1 #6 


4 # 6 + 1 #5 


4#3@ 12.5 cm 


4 # 3 @ 40 cm 


T3-R1 


45x80 


4 # 8 


3 # 6 + 2 # 5 


4 # 3 @ 12.5 cm 


4 # 3 @ 30 cm 


T3-R2 


45x80 


4 # 8 


4 # 6 + 1 #5 


4#3@ 12.5 cm 


4 # 3 @ 35 cm 


T3-R3 


45x80 


4 # 8 


2 # 6 + 2 # 5 


4 # 3 @ 10 cm 


4 # 3 @ 15 cm 


T3-R4 


45x80 


2 # 8 + 4 # 6 


2 # 6 + 2 # 5 


4#3@ 12.5 cm 


4 # 3 @ 20 cm 


T3-R5 


45x80 


2 # 8 + 4 # 6 


5 # 5 


4#3@ 10 cm 


4 # 3 @ 1 5 cm 


T3-R6 


45x80 


2 # 8 + 3 # 6 


5 # 5 


4#3@ 12.5 cm 


4 # 3 @ 20 cm 
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5 ANALISIS DINAMICOS NO LINEALES 

5.1 Consideraciones para los analisis no lineales 

Para efectuar los analisis no lineales se empleo el programa DRAIN-2DX (Prakash etal. 
1992), utilizando modelos que incluyen la distribucion de fuerzas laterales entre los marcos por 
accion del diafragma rfgido (figuras 6 y 7). La accion de la losa como diafragma rigido se 
modelo a partir de elementos elasticos infinitamente rfgidos, teniendo como unica funcion 
transmitir las fuerzas cortantes entre cada uno de los marcos, sin posibilidad de disipar energfa 
por comportamiento inelastico ni por amortiguamiento. Se consideraron los efectos P-A y se 
ignoro la interaccion suelo-estructura. 



0 '0 0 0 0 0 0 0 0 0 
- 0 - 0 - 0 ' 0 0 - 1 - 0 0 0 0 




Figura 6 Modelo con distribucion de fuerzas laterales entre marcos en direccion x 

0000 0000 
0 0 0 0 0 0 0 0 " 0 0 0 0 




Figura 7 Modelo con distribucion de fuerzas laterales entre marcos en direccion y 

Para los analisis no lineales se contemplo tanto las resistencias nominales como 
fuentes de sobrerresistencia y se utilizo la componente E-W del acelerograma registrado 
durante el sismo del 19 de septiembre de 1985 en la Secretarfa de Comunicaciones y 
Transportes (SCT), mostrado en la figura 8, con el fin de tener un acelerograma representative 
y de caracteristicas similares al sismo maximo contemplado por el Reglamento del cual se 
pretenden evaluar sus criterios. 




Evaluacion del criterio que incrementa en 70% la carga axial por sismo en columnas de... 
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bi 



h-1 

W 

U 



0.20 

0.10 

0.00 

- 0.10 

- 0.20 




Tiempo (s) Periodo (s) 

Figura 8 Componente E-W del acelerograma registrado en SCT durante el sismo de 1985 



5. 1. 1 Determination de las resistencias nominales 



Se determinaron las resistencias a flexion de las trabes y flexocompresion y 
flexotension de las columnas de cada modelo considerando: 

• El area real provista de refuerzo, incluyendo la aportacion del refuerzo a 
compresion de las vigas. 

• El nuevo bloque equivalente de esfuerzos del concreto estipulado en la 
propuesta de las NTCC, que es mas parecido al del reglamento ACI-318. 

• Un comportamiento elastoplastico del acero de refuerzo, como se especifica en 
la propuesta de NTCC. 

Aunque se han designado como resistencias nominales a aquellas obtenidas aplicando 
las tres consideraciones arriba mencionadas, es importante mencionar que incluyen cierto 
grado de sobrerresistencia; aquellas inherentes al uso de la propuesta de RCDF. Por lo 
anterior, la respuesta de la estructura con consideracion de resistencias nominales presentara 
efectos de sobrerresistencia. 

Los factores de sobrerresistencia incorporados a este nivel serian: 

• Factor debido al uso de factores de combinacion de carga vertical y lateral (FC 
= 1 . 1 ). 

• Factor debido al uso de de factores de reduction de resistencia asociados a un 
modo de falla ductil por flexion para vigas 1/0.9 = 1.11 y para columnas 1/0.8 = 
1.25). 

• Factor de area real de acero. 

• Factor asociado al uso de secciones y elementos tipo. 

• Contribution del acero a compresion. 

• Efectos tridimensionales. 

• Cargas vivas normativas mayores que las reales. 
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Los factores de sobrerresistencia asociados a la posibilidad de dominio de condiciones 
de carga gravitacional en el diseno, requerimientos de capacidades resistentes mi'nimas, 
requerimientos de rigideces por criterios de deformacion y requerimientos de cuantias de acero 
mi'nimas se encuentran englobados en los factores asociados al uso de secciones y elemento 
tipo y al area real de acero. 



5. 1.2 Determination de las resistencias incorporando fuentes de sobrerresistencia 

Para tener una estimacion de la sobrerresistencia maxima, se determino la resistencia 
de los miembros estructurales considerando las siguientes fuentes adicionales de 
sobrerresistencia: 

• El factor de comportamiento real del acero empleando el modelo de Rodriguez 
y Botero (1994) y, 

• El factor de comportamiento real del concreto utilizando el modelo de Kent y 
Park modificado (Kenty Park, 1971). 

La aportacion de la losa a la resistencia a flexion de las vigas, en particular cuando la 
losa se encuentra en tension, es otra de las fuentes de sobrerresistencia reconocidas; sin 
embargo, existe gran incertidumbre en cuanto aporta en efecto la losa, debido al hecho que el 
ancho y, consecuentemente, el refuerzo de la losa comprendido en este ancho, es variable y 
esta en funcion de la deformacion que el sismo impone a la estructura, lo que dificulta una 
determinacion razonable de la cuantia efectiva de refuerzo de la losa que pudiera incrementar 
la resistencia a flexion de las trabes. Entre mayores sean las rotaciones en las articulaciones 
plasticas adyacentes a las caras de las columnas, mas barras de la losa alejadas de la columna 
contribuiran a la resistencia a flexion. Ademas, tambien interviene en esto la presencia de vigas 
transversales. Por lo anterior, no se tomo en cuenta la participation de la losa en la resistencia 
a flexion de las trabes en este estudio. 

Tampoco se consideraron los factores de sobrerresistencia relacionados con: 

• Resistencias de los materiales mayores que las nominales. 

• Incremento de la resistencias de los materiales ante cargas rapidas. 



5.2 Resultados de los analisis no lineales 

Se estudiaron los siguientes parametros dinamicos para todos los modelos no lineales: 
(a) curvas de histeresis de entrepiso, (b) mapeo de la fluencia de elementos, en tiempos 
especificos, en envolventes, y su monitoreo en el dominio del tiempo, incluyendo la 
determinacion de la magnitud de las rotaciones inelasticas maximas, (c) envolventes de 
distorsiones de entrepiso maximas (A), (d) envolventes de cortantes de entrepiso maximos 
(V/W T ), (e) demandas maximas de ductilidad de entrepiso (p), (f) distorsion de entrepiso 
asociada a la primera fluencia de cualquier elemento estructural que forma parte del entrepiso 
(Ytiuencia o A fe ), (g) distorsion de entrepiso maxima asociada al medio ciclo de histeresis de carga 
maxima (y m axima o A mi ), (h) cociente del numero de medios ciclos donde se presento respuesta 
inelastica (MC ine iasticos) entre el numero total de ciclos a que fue sujeta la estructura ante un 
acelerograma dado (Total MC), (i) cociente entre la rigidez de entrepiso secante minima 
inelastica o “pico a pico” (k m[nim a 6 k imin ) asociada a los ciclos histereticos de amplitud maxima y 
la rigidez elastica del entrepiso (k e i 4s t ica 6 k e i) y, (j) cociente entre el promedio de las rigideces de 
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entrepiso secantes asociadas a los medios ciclos histereticos inelasticos (k prome dio) y la rigidez 
elastica del entrepiso (k e ia S tica 6 k e i). Algunos de estos conceptos se ilustran esquematicamente 
en la figura 9. 





Figura 9 Definition de parametros a partir de ciclos y medio ciclos de histeresis, generalmente 

asimetricos 



Para ambos modelos (DQ4CP y DQ4CP1.7), en direccion xel comportamiento no lineal 
fue mucho mas pronunciado que el correspondiente a la direccion y, cuando se consideran 
resistencias nominales. Ademas, se observaron las siguientes tendencias en el comportamiento 
(figuras 10 a 13): 

• Los marcos interiores tienden a fluir a niveles de distorsion menores a los de los 
marcos exteriores. 

• Las distorsiones de entrepiso maximas correspondientes a los marcos 
interiores son ligeramente mayores a las de los marcos exteriores. 

• Los marcos exteriores tienden a tomar mas cortante que los interiores. 

• Las demandas de comportamiento inelastico son mas altas para los marcos 
interiores. 

• Las mayores demandas de ductilidad estan asociadas a los marcos interiores. 

• Los marcos interiores presentan mayor degradation de rigidez de entrepiso que 
los marcos exteriores. 
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Figura 10 Envolventes de respuesta maxima, marcos A y B del modelo DQ4CP (nominal) 
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Figura 1 1 Envolventes de respuesta maxima, marcos 1 , 2 y 3 del modelo DQ4CP (nominal) 
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Figura 12 Envolventes de respuesta maxima, marcos A y B del modelo DQ4CP1 .7 (nominal) 
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Figura 13 Envolventes de respuesta maxima, marcos 1 , 2 y 3 del modelo DQ4CP1 .7 (nominal) 
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De estas observaciones se concluye que los marcos interiores, en cuyo diseno rigen 
generalmente las combinaciones por cargas verticales, presentan cierto grado de subdiseno por 
sismo cuando estos disenos se basan en los resultados de analisis con modelos elasticos 
tridimensionales, como lo permiten las normas. 

Al incorporar los efectos de las fuentes adicionales de sobrerresistencia en los analisis 
no lineales, se observo que los modelos sufren niveles de distorsion dinamica maxima muy 
similares a cuando se tomo en cuenta las resistencias nominales. La distorsion de entrepiso 
asociada a la primera fluencia aumenta, las distorsiones maximas se abaten notablemente, la 
estructura aumenta su capacidad para tomar mas cortante y la respuesta inelastica disminuye, 
por lo que disminuyen las demandas de ductilidad y la degradacion de rigidez por trabajo 
inelastico. El numero de fluencias tambien disminuye, asi como las demandas de rotacion 
plastica acumulada (figuras 14 a 17). 

Las distorsiones de entrepiso dinamicas maximas de los marcos A, B y el edificio en 
direccion x del modelo DQ4CP rebasan ligeramente a las distorsiones de diseno en los 
entrepisos tres, cuatro y cinco, superando apenas en los entrepisos cuatro y cinco al h'mite 
y=0.012 establecido en cuerpo principal de las NTCS-01 (figuras 18, 19 y 23), estando muy por 
debajo del Ifmite y=0.030 estipulado en el Apendice Normativo A de las NTCS-01 para el 
cumplimiento del estado Ifmite de seguridad contra colapso. Aunque las distorsiones de 
entrepiso de los modelos DQ4CP y DQ4CP1.7 son muy parecidas, las correspondientes al 
modelo DQ4CP1.7 son ligeramente menores a las del modelo DQ4CP en los primeros 
entrepisos, lo que ayuda a que en el entrepiso tres no se rebasen las distorsiones de diseno 
como en el caso del modelo DQ4CP (figuras 18, 19 y 23). En direccion y las distorsiones de 
entrepiso dinamicas no rebasan a las distorsiones de diseno (figuras 20, 21, 22 y 24). En los 
primeros seis entrepisos, las distorsiones de entrepiso dinamicas maximas del modelo 
DQ4CP1.7 son ligeramente mayores a las correspondientes del modelo DQ4CP. 
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Figura 14 Envolventes de respuesta maxima, modelo DQ4CP direccion x (sobrerresistencia) 
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Figura 17 Envolventes de respuesta maxima, modelo DQ4CP1.7 direccion y (sobrerresistencia) 



Tanto en direccion x como en direccion y, las diferencias entre las distorsiones de 
entrepiso dinamicas maximas de los modelos DQ4CP y DQ4CP1.7 son practicamente 
imperceptibles al considerar las fuentes adicionales de sobrerresistencia (figuras 25 y 26). La 
distorsion de fluencia de los modelos DQ4CP y DQ4CP1.7 en direccion xes muy parecida; sin 
embargo, se aprecia que en el primer entrepiso se retrasa ligeramente la fluencia del modelo 
DQ4CP1.7 con respecto a la del modelo DQ4CP, ya que esta se presenta para un valor de 
distorsion mayor (figuras 18, 19 y 23). En direccion y no se aprecian diferencias entre los 
modelos DQ4CP y DQ4CP1 .7 (figuras 20, 21 , 22 y 24). 



Cuando se consideran las fuentes adicionales de sobrerresistencia, la distorsion de 
fluencia del modelo DQ4CP1.7 es mayor que la correspondiente al modelo DQ4CP en el primer 
entrepiso en direccion x (figura 25), sin notarse diferencia alguna en la direccion y (figura 26). 



Las distorsiones maximas en direccion xde ambos modelos rebasan a las distorsiones 
de diseno de los entrepisos uno a nueve, alcanzando valores maximos de y=0.020 en el 
entrepiso cinco. Al comparar los valores de distorsion maxima de los modelos DQ4CP y 
DQ4CP1.7, se observa que en el primer entrepiso, el modelo DQ4CP1.7 presenta un valor 
menor al del modelo DQ4CP (figuras 18, 19 y 23). En la direccion y las distorsiones maximas 
rebasan a las de diseno solo de los entrepisos dos a cinco, alcanzando valores maximos de 
7=0.008 en el entrepiso tres (figuras 20, 21, 22 y 24). Del entrepiso uno al seis, el modelo 
DQ4CP1.7 en direccion y presenta mayores valores de distorsion de entrepiso maxima. 

Considerando las fuentes adicionales de sobrerresistencia, no se aprecian diferencias 
entre los valores de distorsion maxima de los modelos DQ4CP y DQ4CP1.7 en ambas 
direcciones (figuras 25 y 26). 
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Figura 18 Envolventes de respuesta maxima del marco A de los modelos DQ4CP y DQ4CP1 .7 
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Figura 19 Envolventes de respuesta maxima del marco B de los modelos DQ4CP y DQ4CP1.7 
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Figura 20 Envolventes de respuesta maxima del marco 1 de los modelos DQ4CP y DQ4CP1.7 
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Figura 21 Envolventes de respuesta maxima del marco 2 de los modelos DQ4CP y DQ4CP1 .7 
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Figura 23 Envolventes globales en direccion x, modelos DQ4CP y DQ4CP1.7, resist, nominal 
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Figura 25 Envolventes globales en direccion x, modelos DQ4CP y DQ4CP1.7, sobrerresistencia 
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El marco A del modelo DQ4CP1.7 toma en el primer entrepiso mas cortante que el 
marco A del modelo DQ4CP; sin embargo, el marco B del modelo DQ4CP1.7 toma menos 
cortante que su homologo del modelo DQ4CP en el mismo entrepiso (figuras 18 y 19). Al 
considerar el comportamiento global en direccion x, se aprecia que los dos modelos toman 
valores similares de cortante, notandose que el modelo DQ4CP1.7 toma un poco mas de 
cortante en el primer entrepiso que el modelo DQ4CP, mientras que en el quinto entrepiso 
sucede lo contrario (figura 23). 



En direccion y, al comparar los valores de cortante de los modelos DQ4CP y DQ4CP1 .7 
se observan diferencias importantes unicamente en el primer entrepiso, siendo que el cortante 
que toma el marco exterior uno del modelo DQ4CP es mayor al correspondiente al modelo 
DQ4CP1.7 y ocurriendo lo contrario en los marcos interiores dos y tres (figuras 20 a 22). Al 
considerar el modelo global en direccion y, los valores de cortante de entrepiso son muy 
similares (figura 24). 



Al considerar las fuentes adicionales de sobrerresistencia en direccion x, el modelo 
DQ4CP1.7 toma en el nivel uno mas cortante que el modelo DQ4CP, sin que se aprecie en la 
direccion y diferencia alguna (figuras 25 y 26). 



En cuanto a las demandas maximas de accion inelastica en los entrepisos 
(MCineiasticos/Total MC), en ninguna direccion se aprecian diferencias entre los modelos DQ4CP y 
DQ4CP1.7 cuando se considera el comportamiento global de la estructura (figuras 23 y 24). 



Al considerar las fuentes adicionales de sobrerresistencia en direccion x, se observa 
para el primer entrepiso que la actividad inelastica del modelo DQ4CP1 .7 es inferior a la del 
modelo DQ4CP pues el modelo DQ4CP1.7 tiene aproximadamente 2.5% de demanda de 
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medios ciclos inelasticos (MC ine iasticos/Total MC), mientras que el modelo DQ4CP presenta cerca 
del 10% (figura 25). En la direccion y se tienen porcentajes de demandas de medios ciclos 
inelasticos muy bajos, alcanzando algunos entrepisos valores maximos de aproximadamente 
2.5% (figura 26), sin observar diferencias de importancia entre los modelos DQ4CP y 
DQ4CP1.7. 

Las demandas de ductilidad de entrepiso en direccion x se encuentran por debajo del 
factor de comportamiento sismico Q=4, alcanzando valores maximos cercanos a tres en el 
sexto entrepiso. El modelo DQ4CP1.7 presenta en el primer entrepiso una demanda de 
ductilidad de p=1.8, mientras que el modelo DQ4CP tiene en ese mismo entrepiso una 
demanda de ductilidad del orden de p=2.4, sin apreciar diferencias significativas en el resto de 
los entrepisos (figura 23). Se observa una clara correspondencia entre los valores mas altos de 
demanda de ductilidad y la mayor degradacion de rigidez por trabajo inelastico, siendo que a 
las maximas demandas de ductilidad corresponden valores de rigidez inferiores a la mitad de la 
rigidez elastica si se toma la rigidez minima, 6 65% de la rigidez elastica si se utiliza la rigidez 
promedio (figura 23). En direccion y las demandas de ductilidad alcanzan valores maximos de 
tres en el cuarto entrepiso, correspondiendo a estas demandas de ductilidad rigideces minimas 
del orden del 35% de la rigidez elastica y rigideces promedio de aproximadamente 50% de la 
rigidez elastica (figura 24). En esta direccion, el modelo DQ4CP1.7 presenta en el primer 
entrepiso una demanda de ductilidad ligeramente menor a la del modelo DQ4CP. 

Cuando se consideran las fuentes adicionales de sobrerresistencia, se tiene en 
direccion x valores maximos del orden de p=1.5 en el quinto nivel asociados a rigideces 
minimas de aproximadamente 65% de la rigidez elastica y rigideces promedio de 80% de la 
rigidez elastica. La variacion de las demandas de ductilidad de entrepiso con la altura es muy 
similar en los modelos DQ4CP y DQ4CP1 .7, aunque se aprecia que en el primer entrepiso la 
demanda de ductilidad del modelo DQ4CP1.7 es menor (figura 25). En direccion y las 
demandas de ductilidad y degradacion de rigidez se abaten considerablemente, permaneciendo 
la estructura practicamente elastica con demandas maximas de ductilidad apenas mayores a la 
unidad en algunos entrepisos y valores de rigidez minima y promedio practicamente iguales a la 
rigidez elastica (figura 26), sin apreciar diferencias significativas entre los modelos DQ4CP y 
DQ4CP1.7. 

En direccion x, la mayoria de las fluencias se presento en las trabes, fluyendo 
practicamente todas, con exception de las que se encuentran en los ultimos dos niveles y 
algunas del doceavo nivel de los marcos exteriores. Todas las columnas ubicadas en planta 
baja fluyeron en su base, lo cual se considera aceptable para los niveles de desplazamiento 
que el sismo maximo contemplado por el reglamento le demanda a la estructura, dado que la 
magnitud de la rotacion plastica maxima observada es reducida (figura 27). 

Tambien se presentan fluencias en la base de algunas columnas ubicadas en el 
segundo entrepiso, presentando el modelo DQ4CP1.7 menos fluencias en las columnas de este 
nivel que el modelo DQ4CP. Las columnas B1, Cl, B5 y C5 del segundo entrepiso del modelo 
DQ4CP1.7, a diferencia de estas mismas columnas en el modelo DQ4CP, no tuvieron trabajo 
inelastico, seguramente como consecuencia de haber resultado con refuerzos distintos (figura 
27 y tablas 1 y 3). 

Adicionalmente, se observaron en ambos modelos fluencias en la base de las columnas 
A1 , D1 , A5 y D5 del tercer entrepiso y en la parte superior de las columnas A1 , D1 , A2, D2, A3, 
D3, A4, D4, A5 y D5 de los entrepisos seis y siete asi como algunas del octavo entrepiso, 
apareciendo tambien en el entrepiso once fluencias en la base de algunas de las columnas 
ubicadas en los marcos exteriores, nuevamente, la magnitud de la rotacion plastica maxima es 
reducida (figura 27). 
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Figura 27 Rotaciones plasticas acumuladas en los miembros de los modelos DQ4CP y 
DQ4CP1.7 en direccion x, considerando resistencias nominales 
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Figura 28 Rotaciones plasticas acumuladas en los miembros de los modelos DQ4CP y 
DQ4CP1.7 en direccion y, considerando resistencias nominales 
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De hecho, todas las columnas que presentaron fluencias en los modelos DQ4CP y 
DQ4CP1.7, y que arriba se acaban de enumerar, presentaron demandas de rotacion plastica 
acumulada inferiores a 0.01 rad; mientras que las trabes presentaron valores cercanos a 0.09 
rad, observandose claramente una concentracion de trabajo inelastico en los niveles tres, 
cuatro y cinco y en particular en los marcos interiores (figura 27). Estas diferencias en las 
magnitudes de las rotaciones plasticas maximas observadas en vigas y columnas, ademas de 
su correspondiente mapeo, indican que el comportamiento inelastico observado corresponds en 
efecto al asociado a un diseno con mecanismo de falla ultimo de viga debil - columna fuerte. 

Las distintas magnitudes de las rotaciones plasticas en las vigas en elevacion indican 
que estas no son uniformes y que salvo para los niveles inferiores (2 a 8), no son muy altas, lo 
que ilustra claramente la crudeza de la recomendacion de las NTCC y del ACI de considerar 
para fines de analisis y diseno un momento de inercia agrietado en vigas igual al 50% de su 
inercia bruta, lo cual es muy conservador, sobre todo cuando se consideran fuentes de 
sobrerresistencia relacionadas con su detallado como marco ductil, donde las demandas de 
rotacion en vigas disminuyen notablemente, como se aprecia en la figura 29. 

En la direccion y fluyeron todas las columnas del modelo DQ4CP del primer entrepiso 
en su base, mientras que en el modelo DQ4CP1.7 no fluyeron las columnas B1, Cl, B5 y C5 
del primer entrepiso, siendo estas columnas mas reforzadas que las correspond ientes 
columnas del modelo DQ4CP como consecuencia de la aplicacion de la disposicion relacionada 
con el incremento en 70 por ciento de la fraccion de carga axial debida a sismo (figura 28 y 
tablas 1 y 3). En ambos modelos fluyeron en su base las columnas A2, D2, A4 y D4 del 
segundo entrepiso, y en el sexto entrepiso las columnas Al y A5 en ambos extremos, en el 
extremo superior las columnas D1 y D5, y en el extremo inferior las columnas A2, A3 y A4. 

Al igual que en direccion x, las columnas no presentaron rotaciones plasticas 
acumuladas mayores a 0.01 rad, mientras que las vigas no rebasaron el valor de 0.02 rad, 
concentrandose tambien el trabajo inelastico en los niveles dos y tres, predominando este en 
los marcos interiores. 

Al considerar las fuentes adicionales de sobrerresistencia en direccion x, dejan de fluir 
en su base dos columnas del modelo DQ4CP (Al y D1) y cuatro columnas del modelo 
DQ4CP1.7 (Al, B1, Cl y D1) ubicadas todas en el primer entrepiso, estableciendose esta 
diferencia por el diseno distinto de las columnas B1 y Cl de estos modelos (figura 29 y tablas 1 
y 3). En el segundo entrepiso disminuye mas el numero de fluencias en columnas para el 
modelo DQ4CP que para el DQ4CP1.7, ya que en este ultimo permanecen doce columnas 
articuladas en su base, mientras que en el modelo DQ4CP solo ocho. En ambos modelos, se 
mantienen las fluencias en la parte superior de las columnas Al, D1, A2, D2, A3, D3, A4, D4, 
A5 y D5 (entrepiso seis) y A2, D2, A3, D3, A4 y D4 (entrepiso siete). 

Las demandas de rotacion plastica acumulada en columnas considerando la 
sobrerresistencia se mantuvieron tambien por debajo de 0.01 rad. El numero de fluencias y 
demandas de rotacion en trabes disminuyo notablemente, en particular para aquellas trabes 
ubicadas arriba del octavo nivel, disminuyendo todas las rotaciones plasticas acumuladas de tal 
suerte que se aprecian valores maximos de 0.04 rad en algunas trabes de los entrepisos tres y 
cuatro (figura 29). 

Cuando se considerara la sobrerresistencia en direccion y (figura 30) la estructura 
responde casi de manera elastica, presentandose unicamente 12 fluencias en vigas y ninguna 
en columnas, concentrandose estas en los niveles dos, tres, cuatro y ocho del marco interior, 
con valores de rotacion plastica acumulada inferiores a 0.01 rad. 
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Figura 29 Rotaciones plasticas acumuladas en los miembros de los modelos DQ4CP y 
DQ4CP1.7 en direccion x, considerando sobrerresistencia 
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Figura 30 Rotaciones plasticas acumuladas en los miembros de los modelos DQ4CP y 
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Cabe senalar que la distribucion de las fluencias en las columnas del edificio en 
elevacion no es aleatoria, coincidiendo la ocurrencia de estas precisamente con los entrepisos 
en que se redujo en la etapa de diseno la seccion de las columnas (tablas 2 y 4), lo cual lleva a 
pensar que la fluencia de las columnas es entre otras cosas, una consecuencia de la 
discontinuidad en rigidez y resistencia en elevacion. 



6 CONCLUSIONES 

Como se demuestra en este estudio, existen muy pocas diferencias entre los disenos 
que contemplan e ignoran las disposiciones de las secciones 7.3 y 7.3.5. 1 de la propuesta de 
NTCC-01, relacionadas con el incremento en 70 por ciento de la fraccion de la carga axial 
debida a sismo para la obtencion de la carga axial de diseno empleada en el dimensionamiento 
de las columnas de marcos ductiles disenados con Q=4. 

Estas pequenas diferencias se presentan fundamentalmente en las columnas ubicadas 
en los primeros entrepisos y en la periferia del edificio, es decir, en las columnas de borde, pero 
no en las columnas de esquina. En las columnas centrales no se observaron diferencias en las 
demandas de refuerzo longitudinal, independientemente de su ubicacion en elevacion, ni en las 
columnas de los niveles superiores sin importar su ubicacion en planta. 

El hecho que la aplicacion de las disposiciones en cuestion solo tengan repercusion en 
las columnas ubicadas en el primer entrepiso y periferia hace que las diferencias en el 
comportamiento inelastico sean minimas y de caracter local. 

Los resultados de este estudio y los de investigaciones previas (Luna 2000) motivaron 
la cancelacion de este requisito, que aparecio por primera vez en 1987 y se mantuvo vigente 
hasta la propuesta de 2001, de las NTCC-04 vigentes. En trabajos previos (Luna y Tena, 2002) 
se demostro que tal disposicion no esta basada en estudios que metodologicamente sean 
apropiados para justificarla. De hecho, como se ha comentado en otros estudios, esta 
disposicion no existe en ningun otro reglamento de vanguardia a nivel mundial para el diseno 
de marcos ductiles, y dado que muchos de los programas de analisis y diseno mas utilizados en 
Mexico, al ser de procedencia extranjera, no contemplan esta disposicion exclusiva de la 
propuesta de las NTCC hasta antes de 2004, lo que complicaba enormemente el proceso de 
diseno de marcos ductiles y desincentiva el que se disene con marcos ductiles de concreto 
reforzado en Mexico, que es lo que precisamente se desea fomentar. 

Finalmente, el estudio tambien ilustra que aun considerando resistencias nominales, 
una estructura de concreto disenada con criterios para marcos ductiles sujeta a un 
acelerograma intenso asociado al espectro de diseno, como el registrado en SCT durante el 
sismo del 19 de septiembre de 1985, desarrolla magnitudes notablemente distintas de 
rotaciones plasticas en vigas en su elevacion, demostrando que estas no son uniformes y salvo 
para niveles inferiores e intermedios, estas no son altas, lo que ilustra claramente la crudeza de 
la recomendacion de reglamentos de concreto reforzado como las NTCC y el ACI de considerar 
para fines de analisis y diseno una inercia efectiva en vigas igual al 50% de su inercia bruta, lo 
cual es conservador, sobre todo cuando se consideran fuentes de sobrerresistencia 
relacionadas con su detallado como marco ductil, donde las demandas de rotation en vigas 
disminuyen notablemente. 
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RESUMEN 

Se comparan las respuestas en el tiempo, que se obtienen en tres estructuras con 
disipadores de energfa visco elasticos de goma, de tres, seis y nueve pisos, empleando el 
Metodo de la Energfa Modal de Deformacion, el Metodo del Espectro de Capacidad y el 
Metodo N2. En el primer metodo, una vez que se obtiene los factores de amortiguamiento, 
equivalente de la estructura con los disipadores, para cada modo de vibracion se realiza un 
analisis elastico. En cambio en los otros dos metodos se obtiene primero la curva de capacidad 
sfsmica resistente mediante un analisis no lineal estatico, tambien conocida como pushover, en 
que se cambia la matriz de rigidez de acuerdo al dano que experimentan sus elementos y 
luego se aplican los respectivos algoritmos de calculo. 

Cada estructura se analiza ante la accion de tres eventos sfsmicos y son el sismo de El 
Centro de 1940, el de Chile de 1985 y el de Mexico de 1985. Los dos primeros son de alta 
frecuencia y el tercero de baja frecuencia. La comparacion de las respuestas en el tiempo se 
las efectua a nivel de desplazamientos laterales maximos en cada piso y de derivas de piso. 



Palabras claves: Metodo N2, metodo del espectro de capacidad, visco elastico, pushover. 



ABSTRACT 

The answers are compared in time, obtained by three structures with energy 
dissipaters, of three, six and nine floors, using the Deformed Model of Energy Method, the 
Spectral Capacity Method and Method N2. In the first method, once the shock valves are 
obtained, equal to the structure with the dissipaters, for each mode vibration of vibration an 
elastic analysis is carried out however, using the other two methods the resistant seismic 
capacity curve is obtained first by a non-static linear analysis, also know as pushover, in which 
the original rigidity changes according to the damage its elements experimented and after the 
calculus algorithms are applied. 

Each structure is analyzed before the action of the three seismic events and they are 
the seismic of The Centro of 1940, Chile in 1985, and Mexico of 1985. The first two are high 
frequency and the third low frequency. The comparison of the answers at the time they are 
realized at maximum lateral displacement levels in each floor and derived floor. 



Key words: Method N2, Spectrum capacity method, elastic dissipaters, pushover. 
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1. ANTECEDENTES 



Una forma eficiente de reforzar estructuras ante acciones sfsmicas, es mediante la 
incorporacion de disipadores de energfa visco elasticos como el mostrado en la figura 1, que 
esta compuesto por un perfil cerrado y en sus caras exteriores se tiene goma y sobre ellas 
perfiles abiertos. De tal manera que la goma se halla entre dos perfiles de acero, en cada cara 
y la disipacion de energfa se produce por las deformaciones de corte de la goma que esta 
vulcanizada a los perfiles. Aguiar (2007, 1). 

Una de las ventajas que se tienen al utilizar el disipador de energfa visco elastico, 
descrito en la figura 1, es que al incorporarse a la estructura su analisis se realiza con Teorfa 
Lineal, aplicando el Metodo de la Energfa Modal de Deformacion, MEMD. Inaudi et al ( 1993) 
luego que se determina el factor de amortiguamiento equivalente del sistema en cada modo de 
vibracion de la estructura con disipadores de energfa. Aguiar y Prado (2007). 

En el laboratorio se encontro que efectivamente los lazos de histeresis del disipador de 
energfa indicado en la figura 1 , pueden ser descritos por el modelo Kelvin Voight. Aguiar y 
Jimenez (2007). Este modelo establece que el disipador tiene dos componentes una de tipo 
viscoso definida por el producto del amortiguamiento por la velocidad y otra elastica definida 
por el producto de la rigidez por el desplazamiento. 

Al incorporar el disipador de energfa visco elastico a la estructura, se incrementa la 
rigidez del sistema. Aguiar y Beltran (2007), se incrementa el amortiguamiento del sistema lo 
que ocasiona que exista una mayor disipacion de energfa. 




Figura 1 Disipador de energfa utilizado en el estudio. 



En resumen al incorporar el disipador visco elastico, la estructura se mantiene lineal 
pero se ha incrementado la rigidez y amortiguamiento del sistema. Jara y Casas (2002). 

El objetivo de este artfculo es comparar la respuesta en el tiempo, mediante un analisis 
sfsmico piano en tres edificios de tres seis y nueve pisos que se halla mediante el MEMD, con 
los que se obtienen al aplicar el Metodo del Espectro de Capacidad MEC y el Metodo N2. 
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2. METODOS DE ANALISIS 



Al final del apartado anterior, se indico los metodos de analisis sfsmico que se utilizan 
en este artfculo, los aspectos mas importantes de cada uno de ellos se describen a 
continuacion. 



2.1 Metodo de la Energfa Modal de Deformacion, MEMD 



En Aguiar y Beltran (2007) se describe, paso a paso, la forma como se obtiene la 
matriz de rigidez lateral de un portico piano con el disipador de energfa mostrado en la figura 1 
y posteriormente se presenta el calculo de la matriz de rigidez en coordenadas de piso. 



! 

La rigidez axial equivalente k e de la diagonal que contiene al disipador visco elastico 
se halla con las siguientes ecuaciones: 




Donde r/ e es el factor de perdida de la diagonal que contiene al disipador; 77 es el 

1 

factor de perdida de la goma; k es la rigidez al corte de la goma; kj es la rigidez de la 
diagonal. 



K d 



E 

I± + L2 

Ai A 2 



t-ZiA 



A = 4bL 



Siendo E el modulo de elasticidad del acero; Ei,L 2 son las longitudes del perfil 
cerrado y del perfil canal; G a es el modulo de corte de la goma; A es el area de la goma que 

trabaja al corte; e es el espesor de la goma; b es el ancho de la goma y L es la longitud de la 
goma. 



La ecuacion fundamental del MEMD es la determinacion del factor de amortiguamiento 
£ j donde el subfndice j representa el modo de vibracion, con la siguiente expresion: 



~ 1 ' J r\ 



J n 
W j 



2 \ 



W / 



Wj J 
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Donde /?. es el factor de amortiguamiento de la estructura en el modo j; r] e es el 
factor de perdida equivalente descrito en la ecuacion (2); Wj la frecuencia natural de la 
estructura sin disipadores de energfa en el modo j y Wj la frecuencia natural de la estructura 
con disipadores en el modo j . 

El primer termino de la ecuacion ( 5 ) corresponde al amortiguamiento de la estructura, 
sin disipador y el segundo termino al amortiguamiento esperado en los disipadores. En Aguiar y 
Prado (2007) se detalla el calculo. 

Una vez que se tiene definida las matrices de rigidez K , de masas M y el factor de 
amortiguamiento equivalente gj se halla la matriz de amortiguamiento C empleando el 

algoritmo de Wilson y Penzien (1972). 



C = X C i 

i - 1 

Mi 



M* =</>’ M<f>, 



Siendo C,- la matriz de amortiguamiento en el modo i; ^ el factor de amortiguamiento 
hallado con la ecuacion (5); w ; - es la frecuencia natural de la estructura con disipadores de 
energfa en el modo i; <f>j es el modo de vibracion i de la estructura con disipadores. 

Finalmente para encontrar la respuesta en el tiempo se aplico el Metodo denominado 
Procedimiento de Espacio de Estado. Se pudo aplicar otro Metodo pero el P.E.E. es un metodo 
estable. La deduccion y aplicacion de este metodo al igual que el Algoritmo de Wilson y 
Penzien se encuentra en Aguiar (2007, 2). 



2.2 Metodo del Espectro de Capacidad, MEC 



El MEC fue propuesto por Freeman (1975), Freeman et al (1978) y luego fue olvidado, 
en la ultima decada del siglo XX nuevamente se lo ha vuelto a utilizar con una serie de 
modificaciones. El ATC-40 (1996) presento tres procedimientos de calculo con este metodo los 
mismos que han sido criticados por otros investigadores quienes han presentado nuevas 
propuestas a este metodo. Chopra y Goel (1999, 2000, 2002). 

La base del MEC, al igual que del metodo N2, que se vera en el proximo sub apartado, 
es la curva de capacidad sfsmica resistente, que relaciona el desplazamiento maximo en el 
tope, D t de una estructura con el cortante basa, V . Esta curva se halla mediante analisis 
estatico incremental no lineal. Aguiar (2007, 3) desarrollo el programa DISIPA que sirve para 
encontrar la curva de capacidad sfsmica en porticos pianos con disipadores de energfa visco 
elasticos de goma, de la forma indicada en la figura 1 . 
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En la figura 2 se ilustra la forma como se obtiene el Espectro de Capacidad E.C. de 
una estructura, a la izquierda se presenta un portico al cual se le aplican cargas estaticas 
incrementales, en la parte inferior esta el cortante basal V que es igual a la suma de estas 
fuerzas y en la parte superior el desplazamiento lateral D t que generan estas fuerzas. Al 
centra se tiene la curva de capacidad sfsmica resistente y se ha senalado un punto cualquiera 
que tiene coordenadas (D t j,Vj) y a la derecha el E.C. que relaciona el desplazamiento 

espectral Sj con la aceleracion espectral S a . Cada punto del E.C. se halla con las siguientes 
ecuaciones: 



Donde, D rj y, V. son el desplazamiento lateral maximo y el cortante basal asociado, 

para el punto j de la curva de capacidad resistente; y } , es el factor de participacion dinamica 

para el primer modo en el punto j . <j> a , es la amplitud en el tope asociado al primer modo; a l , 

es el factor de participacion del primer modo en el cortante basal y, M T , es la masa total de 
la estructura. En Aguiar (2002) se detalla la forma de calculo del E.C. 

Por otra parte, la demanda sfsmica representada por los espectros, sean estos de 
respuesta o de diseno, vienen en el formato que relaciona el perfodo T con la aceleracion 
espectral S a . Esta demanda se presenta en el formato Sj - S a utilizando la definicion de 
Seudo Espectro. 

Finalmente, se coloca en una sola figura la capacidad de la estructura, representada 
por su E.C. y la demanda sfsmica representada por el espectro de demanda E.D. y se aplica en 
este caso el M.E.C. y en el proximo sub apartado se aplicara el Metodo N2. 

Una vez que se ha colocado en el mismo formato el E.C. y el E.D. se construyen 
espectros inelasticos para diferentes demandas de ductilidad //, empleando las siguientes 



Dt 




Figura 2 Curva de Capacidad y Espectro de Capacidad 




M T a j 




ecuaciones: 
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Donde S de ,S ae son las ordenadas del espectro de demanda elastico; S ai , S di son las 
ordenadas del espectro de demanda inelastico para una ductilidad // y es el factor de 
reduccion por ductilidad. 

Existe una gran cantidad de trabajos sobre R M , un acopio de los mismos se halla en 

Aguiar (2007, 4), Leandro y Santana (2004) pero en este articulo se utiliza la ecuacion 
propuesta por Aguiar et al (2007) quienes encontraron este factor en base a 112 sismos 
artificiales compatibles con los espectros del Codigo Ecuatoriano de la Construccion, CEC- 
2000. Las ecuaciones a las que llegan, son: 

aT (1-0.165//) 
oT (1-0.165 //)+ 4900_ 



Donde el parametro a depende del tipo de suelo y esta definido en la tabla 1 , para los 
cuatro perfiles de suelo que contempla el CEC-2000. 



Tabla 1 Valores c 


el parametro a , Aguiar et al (2007, 1) 


Variable 


Suelo SI 


Suelo S2 


Suelo S3 


Suelo S4 


a 


100500 


91000 


73600 


38900 



R„=1 + (m-1) 



En la figura 3 se describe el MEC, no se presenta el espectro de demanda elastico 
pero si se presentan los espectros de demanda inelasticos para ductilidades de 2.38 y 4. Se 
presenta tambien la curva de capacidad sfsmica. Para el espectro de demanda con // = 2.38 se 

observa que este intercepta el espectro de capacidad en un desplazamiento d t = 0.083 m. Por 
otro lado, al dividir esta cantidad para el desplazamiento de fluencia de la estructura que es 
0.0351 m., se halla que es aproximadamente igual a la demanda de ductilidad /u . 




Figura 3 Description del Metodo del Espectro de Capacidad 



Comparacion de tres metodos de analisis sfsmico de estructuras con disipadores de energfa... 71 



Por lo tanto, el punto de desempeno se halla en forma iterativa hasta cuando la 
ductilidad del espectro de demanda es aproximadamente igual a la ductilidad del espectro de 
capacidad. 



La definicion de espectro esta asociada a un sistema de un grado de libertad. Por lo 
tanto, el desplazamiento que se encuentra en el MEC o en el Metodo N2 esta asociado a un 
grado de libertad. Para encontrar el desplazamiento maximo en el tope del edificio D t hay 
varias ecuaciones de paso, una de ellas es la propuesta por Algan (1982) mediante el 
parametro /?, 



D t = A d, 



2N + 1 



Donde D t , d t son los desplazamientos laterales maximos en el sistema de multiples 

grados de libertad y en el de un grado de libertad; N es el numero de pisos. La ecuacion de 
Algan fue obtenida a partir del analisis, de una viga de corte de seccion constante, modelada 
como un continuo. Existen otros trabajos para determinar p x como el desarrollado por Miranda 
y Reyes (2002) en que analizan una viga de corte acoplada a una de flexion pero de seccion 
variable. 

Otra forma de hallar p x , proviene de la Dinamica de Estructuras, mediante el Factor de 
Participacion Modal, FPM, para seguir con la misma nomenclatura se denomina p\ y se halla 
con la siguiente expresion: 

I * 

A= A, , 



Donde M t es la masa del nivel i; (j) i es la ordenada de la forma modal fundamental en 
el nivel i; N es el numero de pisos; <j) m es la amplitud del modo fundamental en la azotea. 



Metodo N2 



En el Metodo N2 , Fajfar (2000) se obtiene un modelo bilineal del E.C. cuya rama 
elastica cumpla con la condicion de que intercepte la curva del E.C. en un punto cuya ordenada 
espectral sea aproximadamente igual al 60% de la resistencia de plastificacion. A la izquierda 
de la figura 4 se presenta el E.C. y el modelo bilineal de una estructura. 
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0.06 
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Sd <m) 



Figura 4 Description del Metodo N2 



El modelo bilineal y el espectro de demanda elastico se colocan en formato S d - S a 
como se ilustra a la derecha de la figura 4. El punto de desempeno se halla prolongando la 
rama elastica del E.C., hasta que intercepta el espectro de demanda elastico. Para la 
estructura de la figura 4, se halla que el desplazamiento d t = 0.069 m. 



Finalmente se obtiene el desplazamiento en el ultimo piso de la estructura D t 
aplicando la ecuacion de Algan o el F.P.M. Con este desplazamiento D t se hallan los 
desplazamientos en los pisos inferiores ingresando al programa con el cual se obtuvo la curva 
de capacidad sfsmica resistente. 



Cuando se aplica la tecnica del pushover se aplican cargas laterales en cada piso y se 
hallan los desplazamientos laterales en cada piso, que estas producen. Esto para cada 
incremento de carga. Si bien es cierto se obtiene el desplazamiento en el tope del edificio D t y 

el cortante basal V , no es menos cierto que para cada D t se tienen los desplazamientos 
laterales en los pisos inferiores. Las cargas laterales que se aplicaron en este estudio fueron 
proporcionales al primer modo de vibracion. 



3. ESTRUCTURAS ANALIZADAS 



Se analizaron tres estructuras de tres, seis y nueve pisos, cuya configuracion en planta 
es la indicada en la figura 5. Los disipadores de energfa se colocaron en los porticos exteriores 
y se han presentado con linea continua en la figura 5. 
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Figura 5 Planta tipo de las estructuras de 3, 6 y 9 pisos. 




Figura 6 Dimensiones de un portico exterior y ubicacion de los disipadores de energfa. 



La distribucion de los disipadores de energfa visco elasticos, en elevacion se indica en 
la figura 6, para el caso de un edificio de tres pisos. Lo importante es notar la forma en que se 
colocan las diagonales; de igual manera se colocaron para los edificios de seis y nueve pisos. 

El disipador esta compuesto por un perfil tubular cuadrado de 10 cm., de lado y 3 mm., 
de espesor que tiene un area transversal de 1 1 .41 cm 2 . La longitud del perfil tubular es de 4.0 
m. La goma tiene un ancho de 8 cm., un espesor de 2 cm., y una longitud de 2.0 m. Los perfiles 
canales tipo “U” son de 10 cm./ 5 cm., y un espesor de 2 mm; cada perfil tiene un area de 3.87 
cm 2 . Por consiguiente el area de los cuatro perfiles es de 15.48 cm 2 . En la figura 7 se presenta 
la seccion transversal del disipador y a la derecha la longitudinal. La goma tiene un modulo de 

corte G a =1 kg / cm 2 y un modulo de perdida G p =\Akg I cm 2 . Este disipador se utilizo 

tambien en los edificios de seis y nueve pisos. 
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Figura 7 Seccion transversal y longitudinal del disipador de energia. 



Se realizo un analisis sismico piano de cada estructura, para lo cual se analizo el 
portico exterior, que contiene los disipadores de energia; las dimensiones de las vigas y 
columnas al igual que la carga muerta considerada se indica en la tabla 2. 



Tabla 2 Dimensiones de vigas y columnas, de portico exterior. Carga muerta considerada. 



Edificio 


Base y altura de columnas (superior) y vigas (inferior) (cm.) 


Carga 

Muerta 

(T/m2) 




1 nivel 


2 nivel 


3 nivel 


4 nivel 


5 nivel 


6 nivel 


7 nivel 


8 nivel 


9 nivel 


3 pisos 


30/30 

(25/30) 


30/30 

(25/30) 


25/25 

(25/30) 














0.75 


6 pisos 


40/40 

(30/40) 


40/40 

(30/40) 


35/35 

(30/35) 


35/35 

(30/35) 


25/25 

(25/30) 


25/25 

(25/30) 








0.75 


9 pisos 


50/50 

(30/40) 


50/50 

(30/40) 


50/50 

(30/40) 


45/45 

(30/40) 


45/45 

(30/35) 


45/45 

(30/35) 


40/40 

(30/35) 


40/40 

(30/35) 


30/30 

(25/30) 


0.75 





Columnas 1 er nivel y 2do nivel Columnas 3er nivel 

Figura 8 Armadura de columnas de portico exterior de edificio de tres pisos. 



Para aplicar los metodos MEC y N2 se necesita conocer la armadura longitudinal y 
transversal de las columnas y vigas. Para el caso de la estructura de tres pisos, en la figura 8 
se indica la armadura longitudinal de las columnas y en la figura 9 de las vigas. 
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Las estructuras sin disipadores de energia son bastante flexibles, razon por la cual han 
sido reforzadas con disipadores de energia visco elasticos de goma. Lo propio sucede con la 
seccion transversal de las columnas y vigas de los edificios de seis y nueve pisos. En Mieles 
(2007) se indica la armadura de estos edificios. 



3 Pisos g Pisos 






0 0.1 0.2 0.3 0.4 0.5 0.6 



Desplazamiento (m.) 



Figura 10 Capacidad Sismica resistente de Estructuras de 3, 6 y 9 pisos 



En la figura 10 se presentan las curvas de capacidad sismica resistente de las 
estructuras analizadas. La curva inferior corresponde a la estructura sin disipadores de energia 
y las curvas superiores a las estructuras con disipadores de energia, con diferentes espesores 
de goma. Esto unicamente con el objeto de ver como influye el espesor de la goma en la 
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capacidad sfsmica, las restantes dimensiones del disipador se mantuvieron constantes y son 
las indicadas en la figura 7. 

A mayor espesor de goma, la rigidez inicial del sistema disminuye pero la resistencia es 
mayor. Cuando el espesor de la goma es mayor a 3 cm., se debe colocar en la parte central de 
la goma una placa de acero. 

Como se indico anteriormente, en el estudio realizado se trabajo con gomas de 2 cm., 
de espesor. Con esta acotacion, en la figura 11, se indican los E.C. de las estructuras 
analizadas. 




Figura 1 1 Espectro de capacidad de estructuras de 3, 6 y 9 pisos. 



4. SISMOS DE ANALISIS 



Las estructuras van a ser analizadas ante tres eventos sismicos que son: el sismo de 
El Centro, California, registrado el 18 de mayo de 1940, se trabaja con la componente N-S que 
tuvo una aceleracion maxima del suelo igual a 0.348 g., siendo g., la aceleracion de la 
gravedad; se lo identificara con la letra C. 




Perfodo ( s ) 

Figura 12 Espectros de respuesta elasticos de sismos seleccionados para el estudio. 

Ante el sismo de Chile registrado el 3 de marzo de 1985, se trabaja con un registro de 
Vina del Mar, que tuvo una aceleracion maxima del suelo igual a 0.363 g., se lo identifica como 
Ch. y ante el sismo de Mexico del 19 de septiembre de 1985, se trabaja con el registro de la 
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componente E-W, de la Secretarfa de Comunicaciones y Transporte, del Distrito Federal y tuvo 
una aceleracion maxima de 0.184 g., se lo identifica con la letra M. Los espectros, de respuesta 
elastica, en el formato clasico de estos acelerogramas se los indica en la figura 12, para un 
factor de amortiguamiento £ = 0.05 

Los sismos de El Centro y Chile., son de alta frecuencia; en cambio, el sismo de 
Mexico, es de baja frecuencia. Los espectros de la figura 12 se pasan al formato S d -S a y 
son los indicados en la figura 13. 




Figura 13 Espectros de demanda de sismos seleccionados en el formato S d - S a 




Figura 14 Punto de Demanda de las estructuras analizadas por el MEC. 
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5. METODO DEL ESPECTRO DE CAPACIDAD 



En la figura 14 se presenta el punto de desempeno encontrado en cada uno de los 
edificios analizados, en la parte superior se muestra la respuesta para el edificio de 3 pisos, en 
la parte central para el de 6 pisos y en la inferior para el de 9 pisos. Se destaca que es 
respuesta maxima en un sistema de un grado de libertad. 

Por otra parte, en la primera columna de la figura 1 , se tiene la respuesta ante el sismo 
del Centro, en la segunda columna ante el sismo de Chile y en la tercera columna ante el sismo 
de Mexico. 

En cada una de las graficas de la figura 14, se indica la ductilidad con la cual se 
alcanzo el punto de desempeno y el correspondiente desplazamiento d t . No se dibujo todo el 
E.C. y el E.D., sino unicamente de la parte en la cual se encuentra el punto de desempeno. 

En la tabla 3 se presenta los valores de /u y d t , correspondientes al punto de 
demanda, se indica tambien el desplazamiento lateral maximo en el tope del edificio D t que se 
halla al utilizar la ecuacion de Algan y con el F.P.M. 



Tabla 3 Desplazamientos maximos hallados con el MEC. 



Edit. 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 




M 


d t 


D t 


D t 


M 


d t 


D t 


D t 




d t 


D t 


D, 






(m.) 


Algan 

(m.) 


FPM 

(m.) 




(m.) 


Algan 

(m.) 


FPM 

(m.) 




(m.) 


Algan 

(m.) 


FPM 

(m.) 


3 P. 


2.38 


0.083 


0.107 


0.104 


1.50 


0.051 


0.066 


0.064 


0.83 


0.026 


0.034 


0.033 


6 P. 


2.50 


0.116 


0.161 


0.157 


2.55 


0.114 


0.158 


0.154 


1.15 


0.055 


0.076 


0.075 


9 P. 


1.10 


0.095 


0.136 


0.130 


3.60 


0.33 


0.469 


0.445 


2.20 


0.284 


0.284 


0.270 



6. METODO N2 



En el mismo formato de la figura 14, se presenta en la figura 15 el punto de demanda 
encontrado aplicando el metodo N2 y en la tabla 4 se presentan las respuestas maximas tanto 
en el sistema de un grado de libertad d t como en el sistema de multiples grados de libertad 

D t - 



Tabla 4 Desplazamientos maximos hallados con el Metodo N2 



Edit. 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 




M 


d t 


D t 


D t 


M 


d t 


D t 


D t 


M 


df 


D t 


D, 






(m.) 


Algan 

(m.) 


FPM 

(m.) 




(m.) 


Algan 

(m.) 


FPM 

(m.) 




(m.) 


Algan 

(m.) 


FPM 

(m.) 


3 P. 


2.37 


0.069 


0.089 


0.087 


1.37 


0.040 


0.051 


0.050 


0.72 


0.021 


0.027 


0.026 


6 P. 


2.35 


0.113 


0.156 


0.153 


2.80 


0.135 


0.187 


0.182 


1.04 


0.050 


0.069 


0.068 


9 P. 


1.41 


0.102 


0.145 


0.138 


3.25 


0.235 


0.334 


0.317 


1.59 


0.115 


0.163 


0.155 
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Figura 15 Punto de Demanda de las estructuras analizadas por el metodo N2 



7. METODO DE LA ENERGIA MODAL DE DEFORMACION 



Con el metodo MEMD se encuentra directamente los desplazamientos en cada uno de 
los pisos de la estructura. Las respuestas maximas de desplazamiento, se indican en la figura 
16. A diferencia de los metodos anteriores, en el MEMD no se necesita conocer la armadura 
longitudinal y transversal de sus elementos, unicamente con las dimensiones de vigas y 
columnas se realiza un analisis elastico siguiendo el procedimiento descrito en el apartado 2.1, 
de ahf que es el mas utilizado. 

Por lo que se realiza un analisis elastico es el interes que se tiene en comparar las 
respuestas de este metodo con el MEC y el N2, que de alguna manera incorporan el analisis 
no lineal. 



80 



Roberto Aguiar y Yordi Mieles 









7“ 

■ ■ 






\ 

\ 




^ 

\ 


y 







M 





J 


"7 


/ 


\ 




y 







D.DD 0.02 



0 . 0 * 0.06 

C 



0.08 



0.00 







* 








/ 






/ 








/ 






/ 








/ 









0.02 0 . 0 * 

Ck 



0.00 







7 








/ 






/ 








/ 






/ 








/ 









0.01 0.02 

M 



0.03 







r 








/ 






/ 














f 








/ 









0.00 0.05 



3 
8 
T 
6 
5 
* 

3 
2 
1 
0 
0.00 



0.10 

C 



0.15 0.20 





7 




/ 




y 




/ 


/ 




y 




y 




y 









8 * 



0.00 0.05 0.10 0.15 0.20 0.00 0.02 0 . 0 * 0.06 0.08 

M 



0.10 



0.20 







/ 






y 






y 






y 




y 






y 










y 






y 









7 




t 


y 




y 




/ 




/ 




/ 




/ 




y 





*.*• *. 1 * *. 2 * *. 3 * 0.00 

DesDlazami&nto (m.) 



o.io 



o.ao 



Figura 16 Desplazamientos maximos, en cada piso, hallados con el MEMD. 

8. COMPARACION DE DESPLAZAMIENTOS MAXIMOS 



En la figura 17 se presentan los desplazamientos laterales maximos, en cada uno de 
los pisos, hallados con los tres metodos de analisis sfsmico: MEMD, MEC y N2. Para los dos 
ultimos metodos se indica la respuesta encontrada con la ecuacion de Algan y con el F.P.M. 

En la mayor parte de los casos se aprecia, que los menores desplazamientos se 
encuentran con el MEMD, razon por la cual en las tablas 5 a 7, se presenta en porcentaje la 
diferencia de desplazamientos que se halla con los otros metodos, con relacion al MEMD. 



Tabla 5 Diferencia de desplazamientos con los metodos MEC y N2 



PIS 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 


os 


MEC 


N2 


MEC 


N2 


MEC 


N2 




ALGA 


F.P.M 


ALGA 


F.P.M. 


ALGA 


F.P.M 


ALGA 


F.P. 


ALGA 


F.P. 


ALG 


F.P. 




N 




N 




N 




N 


M. 


N 


M. 


AN 


M. 


3 


25% 


25% 


8% 


8% 


31% 


31% 


4% 


4% 


35% 


31% 


16% 


12% 


2 


23% 


8% 


5% 


5% 


30% 


8% 


4% 


4% 


33% 


8% 


15% 


12% 


1 


28% 


4% 


11% 


11% 


34% 


4% 


9% 


9% 


34% 


4% 


17% 


14% 
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Figura 17 Comparacion de desplazamientos maximos, obtenidos con los tres Metodos. 



En las tablas 5 a 7 se aprecia, que las respuestas maximas halladas con el MEC y con 
el N2, son bastante parecidas y cuando hay diferencia el metodo N2 reporta menores valores 
que el MEC. 



Tabla 6 Diferencia de desplazamientos con los metodos MEC y N2 
con relacion a MEMD, para edificio de 6 pisos. 

PIS I SISMO CENTRO ^ SISMO CHILE ^ SISMO MEXICO 

OS MEC ~1 N2 MEC T N2 MEC I ~N2 

ALGA I F.P.M ALGA I F.P.M. ALGA f F.P.M ALGA f F.P. ALGA I F.P. ALG I F.P. 

N . N N . N M. N M. AN M. 

6 20% 20% 20% 20% 9% 9% 9% 9% 26% 25% 19% 17% 

5 22% 22% 22% 22% 12% 12% 12% 12% 27% 26% 20% 18% 

4 22% 22% 22% 22% 13% 13% 13% 13% 30% 29% 22% 21% 

3 21% 21% 21% 21% 12% 12% 12% 12% 31% 30% 24% 22% 

2 25% 25% 25% 25% 16% 16% 16% 16% 35% 33% 27% 26% 

1 29% 29% 29% 29% 20% 20% 20% 20% 37% 36% 29% 28% 
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Por otra parte al comparar los resultados que se hallan con la ecuacion de Algan o con 
el F.P.M., para hallar los desplazamientos en el sistema de multiples grados de libertad a partir 
de los resultados del sistema de un grado de libertad. Se aprecia que son bastante parecidos a 
excepcion de los resultados obtenidos en el sismo de 9 pisos en que hay una notable diferencia 
en el Metodo N2 para el caso del sismo de Chile. 



Tabla 7 Diferencia de desplazamientos con los metodos MEC y N2 



con relacion a 



MEMD, para edificio de 9 pisos. 



PIS 

OS 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 


MEC 


N2 


MEC 


N2 


MEC 


N2 


ALGA 

N 


F.P.M 


ALGA 

N 


F.P.M. 


ALGA 

N 


F.P.M 


ALGA 

N 


F.P. 

M. 


ALGA 

N 


F.P. 

M. 


ALG 

AN 


F.P. 

M. 


9 


-17% 


-17% 


-8% 


-8% 


37% 


37% 


22% 


5% 


64% 


64% 


36% 


36% 


8 


-15% 


-15% 


-6% 


-6% 


38% 


38% 


23% 


6% 


64% 


64% 


35% 


35% 


7 


-12% 


-12% 


-4% 


-4% 


38% 


38% 


24% 


7% 


64% 


64% 


35% 


35% 


6 


-7% 


-7% 


1% 


1% 


39% 


39% 


25% 


8% 


64% 


64% 


34% 


34% 


5 


-5% 


-5% 


3% 


3% 


40% 


40% 


26% 


10% 


64% 


64% 


33% 


33% 


4 


-6% 


-6% 


1% 


1% 


41% 


41% 


28% 


11% 


64% 


64% 


32% 


32% 


3 


-8% 


-8% 


0% 


0% 


43% 


43% 


29% 


13% 


64% 


64% 


30% 


30% 


2 


-10% 


-10% 


-2% 


-2% 


45% 


45% 


32% 


16% 


65% 


65% 


28% 


28% 


1 


-3% 


-3% 


4% 


4% 


55% 


55% 


43% 


28% 


69% 


69% 


31% 


31% 










/' N ' 
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*. 
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Figura18 Comparacion de derivas de piso, halladas con los tres Metodos. 
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9. COMPARACION DE DERIVAS DE PISO MAXIMA 



A nivel de derivas de piso, el comportamiento es similar a la de los desplazamientos, 
las menores derivas de piso se hallan con el MEMD y las mayores derivas con el MEC. Los 
valores que reporta el metodo N2, se encuentran en la mitad de estos valores, como se aprecia 
en la figura 18. Las diferencias de valores de los metodos MEC y N2 con respecto al MEMD se 
indican en las tablas 8 a 10 para los edificios de tres, seis y nueve pisos, respectivamente. 

Tabla 8 Diferencia de derivas de piso con los metodos MEC y N2 



con relacion al MEMD, para edificio de 3 pisos. 



PISOS 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 


MEC 


N2 


MEC 


N2 


MEC 


N2 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


3 


31% 


31% 


19% 


19% 


36% 


36% 


1% 


1% 


40% 


32% 


16% 


13% 


2 


17% 


17% 


-3% 


-3% 


25% 


25% 


-1% 


-1% 


32% 


30% 


14% 


10% 


1 


28% 


28% 


11% 


11% 


34% 


34% 


9% 


9% 


34% 


33% 


17% 


14% 



Tabla 9 Diferencia de derivas de piso con los metodos MEC y N2 
con relacion al MEMD, para edificio de 6 pisos. 



PISOS 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 


MEC 


N2 


MEC 


N2 


MEC 


N2 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


6 


-8% 


-8% 


-8% 


-8% 


-25% 


-25% 


-25% 


-25% 


17% 


16% 


10% 


8% 


5 


22% 


22% 


22% 


22% 


5% 


5% 


5% 


5% 


15% 


13% 


7% 


5% 


4 


26% 


26% 


26% 


26% 


17% 


17% 


17% 


17% 


27% 


25% 


19% 


17% 


3 


15% 


15% 


15% 


15% 


5% 


5% 


5% 


5% 


25% 


23% 


17% 


16% 


2 


22% 


22% 


22% 


22% 


13% 


13% 


13% 


13% 


33% 


32% 


26% 


24% 


1 


29% 


29% 


29% 


29% 


20% 


20% 


20% 


20% 


37% 


36% 


29% 


28% 



Tabla 10 Diferencia de derivas de piso con los metodos MEC y N2 
con relacion al MEMD, para edificio de 3 pisos. 



PISOS 


SISMO CENTRO 


SISMO CHILE 


SISMO MEXICO 


MEC 


N2 


MEC 


N2 


MEC 


N2 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


ALGAN 


F.P.M. 


9 


-54% 


-54% 


-43% 


-43% 


21% 


21% 


6% 


-12% 


61% 


61% 


41% 


41% 


8 


-53% 


-53% 


-39% 


-39% 


28% 


28% 


13% 


-4% 


63% 


63% 


43% 


43% 


7 


-50% 


-50% 


-38% 


-38% 


31% 


31% 


12% 


-6% 


62% 


62% 


38% 


38% 


6 


-20% 


-20% 


-11% 


-11% 


34% 


34% 


17% 


-1% 


63% 


63% 


39% 


39% 


5 


0% 


0% 


8% 


8% 


35% 


35% 


21% 


4% 


64% 


64% 


38% 


38% 


4 


-2% 


-2% 


6% 


6% 


37% 


37% 


23% 


7% 


64% 


64% 


37% 


37% 


3 


-5% 


-5% 


3% 


3% 


38% 


38% 


24% 


8% 


63% 


63% 


34% 


34% 


2 


-14% 


-14% 


-6% 


-6% 


37% 


37% 


23% 


5% 


61% 


61% 


25% 


25% 


1 


-3% 


-3% 


4% 


4% 


55% 


55% 


43% 


28% 


69% 


69% 


31% 


31% 



Una de las causas por las cuales los desplazamientos y derivas de piso son mayores 
con el MEC y con el N2 radica en que estos metodos parten de la curva de capacidad sfsmica 
resistente, que relaciona el cortante basal con el desplazamiento lateral maximo, La misma que 
se obtiene aplicando analisis no lineal estatico. 
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En el analisis no lineal estatico se va cambiando la matriz de rigidez de cada elemento 
de acuerdo al nivel de solicitaciones que gravitan sobre ella, de acuerdo al grado de dano, de 
ahf que es no lineal. Pero al ser estatico se aplican cargas laterales incrementales en cada piso 
y se encuentran los desplazamientos en funcion de la matriz de rigidez que tiene el sistema. 

En ningun momenta interviene la matriz de amortiguamiento del sistema en el analisis 
no lineal estatico y una de las caracteristicas fundamentales de los disipadores de energia es 
que se incrementa el factor de amortiguamiento del sistema y por ende la matriz de 
amortiguamiento. En el MEMD si se considera la matriz de amortiguamiento en el analisis 
sismico por lo que se tienen menores desplazamientos y derivas de piso con relacion a los 
otros metodos. 



10. CONCLUSIONES 



Se ha realizado el analisis sismico piano de tres estructuras, de tres, seis y nueve 
pisos; compuestas por vigas, columnas y disipadores de energia visco elastico, ante la accion 
de tres eventos sismicos, empleando los siguientes procedimientos: Metodo de la Energia 
Modal de Deformacion; Metodo del Espectro de Capacidad y Metodo N2; con el proposito de 
comparar las respuestas maximas en desplazamientos y en derivas de piso. Del estudio 
realizado se desprenden las siguientes conclusiones: 

■ El Metodo del Espectro de Capacidad es el que reporta los mayores desplazamientos 
laterales y derivas de piso. Le sigue el Metodo N2. Esto se debe a que los dos metodos 
parten de la curva de capacidad sismica resistente la misma que se obtiene sin tomar en 
cuenta a la matriz de amortiguamiento ya que se la halla mediante un analisis no lineal 
estatico. 

■ En estructuras con disipadores de energia la matriz de amortiguamiento tiene un papel 
fundamental en la respuesta sismica ya que los disipadores incrementan esta matriz y por 
consiguiente se obtienen menores respuestas como se ha visto al aplicar el Metodo de la 
Energia Modal de Deformacion. 

■ Tanto en el Metodo del Espectro de Capacidad como en el Metodo N2, se obtiene el 
desplazamiento maximo en un sistema de un grado de libertad y para encontrar el 
desplazamiento maximo en el sistema de multiples grados de libertad se puede utilizar la 
ecuacion propuesta por Algan o mediante el calculo del factor de participacion modal, que 
dan resultados parecidos. Con el Factor de Participacion Modal se obtienen ligeramente 
menores resultados. 
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RESUMEN 

En este trabajo se estudia el comportamiento sismico de los edificios con forjados 
reticulares, tan habituales en Espana, y se compara con el de otra tipologia estructural 
desaconsejada en zonas sismicas: la de los edificios con vigas planas. Como referencia se 
utilizan los edificios con porticos resistentes a momento tipologia con un adecuado 
comportamiento sismico esperado. Los edificios mencionados se han proyectado aplicando los 
requisitos de la instruccion espanola (EHE) y de la norma NCSE-02. El estudio se Neva a cabo 
analizando la respuesta no lineal de dichos edificios que permite obtener los valores de 
ductilidad asi como verificar los valores de los factores de reduction aplicados en la 
determination de las fuerzas sismicas y los valores de reserva de resistencia. En el proyecto 
sismorresistente de edificios es necesario conocer previamente el valor de la ductilidad que 
estos podran llegar a alcanzar al ser sometidos a movimientos sismicos fuertes. Los valores de 
ductilidad de referencia estan incluidos en las normas de diseno sismorresistente. Su 
estimation tambien puede hacerse bajo criterio de expertos o de la observation de la 
respuesta que los edificios han tenido ante determinados terremotos. Sin embargo, 
generalmente el problema estudiado se ha centrado en la respuesta de edificios con respuesta 
ductil, como son los edificios porticados resistentes a momentos, sin que existan muchos datos 
que avalen la respuesta de los edificios de ductilidad limitada, entre los que se encuentran los 
edificios con forjados reticulares. Finalmente, en el articulo se incluyen comentarios sobre 
defectos de proyecto conceptual observados en algunos edificios existentes que, si bien no 
han sufrido danos por terremotos, se consideran particularmente vulnerables debido a dichos 
defectos de proyecto. 

Abstract 

In this article, the seismic behaviour of waffled slab buildings, very frequent in Spain, is 
studied and compared with that of the buildings with wide beams. Both of them are not 
adequate in seismic areas. As a reference for an adequate seismic design is used the building 
typology with moment resisting frames. The mentioned buildings have been designed by using 
the Spanish codes for reinforced concrete structures (EHE) and for earthquake resistant design 
(NCSE02). The study has been performed by analyzing the non linear response of the 
buildings, what allowed obtaining their ductility, as well as verifying the reduction coefficients 
applied in the calculation of the seismic forces and the overstrength. In the earthquake resistant 
design, it is necessary to know a priori the ductility which the buildings would reach subjected to 
strong ground motions. Reference ductility values are included in the seismic codes. They can 
be estimated using expert opinions or the observed behaviour of buildings during certain 
earthquakes. Nevertheless, most of the studies focussed on the behaviour of buildings with a 
ductile seismic response, like the buildings with moment resisting frames. There are not 
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available data on the response of limited ductility buildings, such as the buildings with waffled 
slabs and the study of these aspects is the main objective of the article. Finally, the article 
includes comments on conceptual design errors observed in existing waffled slabs buildings 
which are found to be extremely vulnerable to seismic actions. 



1. INTRODUCCION 



El desarrollo de procedimientos claros con base en prestaciones (ATC-40, FEMA-273) 
aplicables a la adecuacion y proyecto sismorresistente de edificios, ha incrementado el interes 
de los investigadores sobre la respuesta estatica no lineal de las estructuras (Bertero y Bertero 
2002). Entre las tipologias mas estudiadas se encuentra la de los edificios porticados con vigas 
de canto (Elnashai y Mwafi 2002; Fragiacomo et al. 2006); sin embargo, los edificios que se 
proyectan para que tengan una respuesta poco ductil, denominados edificios de ductilidad 
limitada, entre los que se encuentran los edificios con forjados reticulares y los edificios 
porticados con vigas planas, han sido menos estudiados (Erberik y Elnashai, 2006). 

Existen dos caracteristicas de la respuesta no lineal que permiten establecer si la 
respuesta estructural es adecuada para una determinada amenaza sismica: la ductilidad 
estructural y la reserva de resistencia. Las normas de proyecto sismorresistente, en general, y 
la norma espanola NCSE-02, en particular, contemplan valores de ductilidad muy bajos para 
los edificios de ductilidad limitada, los cuales han sido fijados considerando que dichas 
estructuras disponen de una baja capacidad de disipacion de energia. Por otro lado, a 
excepcion del UBC-97 y del IBC-2003, no existen en las normas de proyecto sismorresistente 
de edificios referencias directas a los valores de la reserva de resistencia que son tan 
importantes en la determination de los factores de reduction de respuesta, (Vielma et al. 2006 
y 2007). El principal objetivo de este trabajo es el de calcular los valores de la ductilidad y de la 
reserva de resistencia de los edificios de ductilidad limitada, aplicando el analisis estatico no 
lineal con control de fuerzas para obtener las curvas de capacidad para las estructuras. Los 
desplomes correspondientes al punto de plastificacion se obtienen utilizando la forma bilineal 
idealizada de dicha curva (Park, 1988). Tambien se utiliza el concepto de indice de dafio 
estructural con el fin de obtener los valores de los desplomes ultimos de los edificios. Conocida 
la respuesta no lineal, se estudia el efecto beneficioso que sobre esta tendria la mejora de las 
caracteristicas de ductilidad del acero de refuerzo, del confinamiento longitudinal y transversal 
y la aplicacion de tipologias estructurales que logren mejorar las caracteristicas de los edificios 
de ductilidad limitada con las de los edificios porticados resistentes a momentos. Finalmente, 
se realiza una comparacion entre la respuesta no lineal de los edificios con ductilidad limitada y 
la respuesta de un edificio porticado resistente a momentos encontrando que, en el caso de los 
primeros, no se satisfacen algunas de las hipotesis consideradas en la etapa de proyecto 
sismorresistente. 



2. DESCRIPCION DE LOS EDIFICIOS ESTUDIADOS 



A continuation se consideran tres edificios de hormigon armado proyectados para 
distintos niveles de ductilidad, conforme a la norma sismorresistente espanola NCSE-02. El 
primer edificio es de forjados reticulares de 30 cm de canto, cuyos nervios se encuentran 
orientados segun las h'neas que unen los extremos de los pilares. El edificio tiene tres niveles, 
el primero con una altura de 4,5 m, mayor que la del resto de los niveles, que tienen 3,0 m. Los 
pilares de este edificio no se encuentran alineados, como puede observarse en la Figura la. 
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Figura 1. Planta tipica del edificio con forjados reticulares (igual en todos los niveles). 



El segundo edificio, proyectado para ductilidad de 2, es porticado con vigas planas y 
forjados unidireccionales orientados segun el eje y se muestra en la Figura 1b. El tercer edificio 
es porticado con vigas de canto, proyectado para ductilidad de 4, con forjados unidireccionales 
orientados segun el eje y, tal como puede verse en la Figura 1c. 



Los edificios descritos anteriormente se han calculado mediante modelos 2D, 
definiendose porticos representatives para cada uno de ellos. Dada la presencia de pilares no 
alineados en el edificio con forjados reticulares, se ha modelado uno de los porticos exteriores 
de este edificio, ya que el resto de los elementos no constituyen un sistema estructural que se 
pueda modelar como piano en sentido estricto, pues aparece el efecto de torsion. 
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3. ANALISIS NO LINEAL 

Efectuando el analisis modal previsto en la norma NCSE-02, se han determinado las 
fuerzas sismicas de calculo a aplicar en los nudos de los porticos del edificio. Para ello se 
utiliza el espectro de inelastico que puede verse en la Figura 2 que se calcula de manera 
simplificada reduciendo el espectro de proyecto elastico por un factor correspond iente a la 
ductilidad y al amortiguamiento (Barbat et al. 2005). Dicho espectro inelastico es el mismo que 
se utiliza en el calculo de las fuerzas sismicas del edificio con forjados reticulares. Las fuerzas 
sismicas maximas probables se calculan mediante la aplicacion del metodo de la raiz cuadrada 
de la suma de los cuadrados (SRSS). 




0,0 1,0 2,0 
Periodo (s) 

Figura 2. Espectro inelastico obtenido de acuerdo con la norma NCSE-02. 



El analisis no lineal tiene como principal objetivo el de hacer una evaluacion mas 
ajustada a la realidad de la respuesta de los edificios proyectados de acuerdo con el metodo de 
calculo lineal elastico previsto en la norma NCSE-02. De esta manera se podra ilustrar, por una 
parte, la manera en que el proyecto sismorresistente mejora la capacidad ductil de las 
estructuras y, por otra, de que manera la respuesta no lineal cuestiona ciertas simplificaciones 
que se efectuan en el analisis elastico. Se debe mencionar tambien que dichas simplificaciones 
no siempre tienen un significado del todo claro para el proyectista de estructuras, (Barbat et al. 
2005). 



3.1. Procedimiento de calculo: analisis al empuje incremental 

Para evaluar la respuesta inelastica de las tres estructuras consideradas se aplica un 
procedimiento de empuje incremental, habitualmente conocido como pushover analysis. Dicho 
procedimiento consiste en someter la estructura a un patron de fuerzas laterales, 
representativas de la accion sismica, distribuidas en altura de acuerdo con diferentes criterios, 
e incrementarlas hasta que se produce el colapso estructural. Entre los criterios mencionados, 
el mas habitual es el de distribucion de acuerdo con la forma del primer modo de vibracion, en 
el que las fuerzas de cada nivel son directamente proporcionales a la altura del edificio (vease 
la Figura 3), debido a que este patron de fuerzas suele ser el mas desfavorable y representa de 
forma bastante correcta la respuesta de las estructuras regulares en planta y elevacion. Las 
fuerzas laterales se incrementan progresivamente, desde el valor nulo, pasando por el que 
produce el cambio de comportamiento estructural elastico a plastico, hasta alcanzar el valor de 
desplome ultimo, a partir del cual la estructura ya no es capaz de soportar ningun incremento 
de carga y su colapso es inminente. Antes de aplicar las cargas laterales representativas de la 
accion sismica, la estructura es sometida a la accion de las cargas de gravedad, de acuerdo 
con las combinaciones aplicadas en el analisis elastico. El metodo descrito de calculo no lineal 
estatico no permite considerar los efectos torsionales. 
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Figura 3. Distribucion de fuerzas aplicada en el analisis no lineal. 



A partir de la respuesta estatica no lineal obtenida mediante la aplicacion de tecnicas 
de elementos finitos, se obtiene la expresion idealizada bilineal de la Figura 4, que tiene un 
segmento representative del comportamiento elastico que va desde el origen hasta el valor que 
corresponde al 75% del cortante maximo en la base (Park 1988). Para obtener la segunda 
recta, que representa la rama de comportamiento plastico, ha sido necesario encontrar la 
intersection de este segmento con uno horizontal, correspondiente al valor maximo del 
cortante en la base. Mediante este procedimiento de compensation se logra garantizar que la 
energia disipada por el sistema modelizado y el ideal sean iguales (Figura 4). 





Figura 4. Esquema para la determination del indice de dano de un elemento estructural. 



Del analisis no lineal interesan dos coeficientes que caracterizan la calidad de la 
respuesta sismorresistente de los edificios. La primera de estas es la ductilidad estructural, 
definida como (Barbat et al. 2007) 




( 1 ) 



y que se calcula a partir de los valores del desplome de plastificacion, Ay, y del desplome 
ultimo, Au, de la curva de capacidad idealizada que aparecen en la Figura 4. 



Tambien resulta interesante comparar los valores del de cortante V p de proyecto en la 



92 



Alex H. Barbat, Juan Carlos Vielma y Sergio Oiler 



base, con el coeficiente de cortante V y de plastificacion en la base, ambos tambien 
representados en la Figura 4. La relacion entre ambas fuerzas cortantes define la reserva de 
resistencia del edificio, R r 



Rr 




( 2 ) 



Para efectos de proyecto, la reserva de resistencia R r puede interpretarse como un factor de 
seguridad. 



3.2. Respuesta no lineal del edificio con forjados reticulares 

En la literatura especializada no abundan los calculos de la respuesta no lineal de este tipo 
de edificios, dada la naturaleza especial del sistema estructural. En este caso se ha adoptado 
un modelo mecanico similar al de portico equivalente descrito en la Norma ACI-318 (2005). 

Este tipo de estructuras se proyecta para un nivel de ductilidad muy bajo en relacion con el 
permitido para otros sistemas estructurales como, por ejemplo, el clasico porticado. En la 
Figura 5 se muestra la curva de capacidad calculada para este edificio, utilizando elementos 
finitos con modelos constitutivos no lineales de dano y plasticidad y la teoria de mezclas de 
materiales. Para controlar la disipacion de energi'a y su correcto comportamiento se han 
introducido valores medios aproximados de la resistencia y la energi'a de fractura de cada 
material constituyente (acero-hormigon) (Car et al. 2000 y 2001). 




Figura 5. Idealizacion bilineal de la curva de capacidad del portico exterior del edificio con 
forjados reticulares. 



La ductilidad estructural calculada para este portico exterior mediante la Ecuacion (1) 
de la curva de capacidad idealizada que puede verse en la Figura 5 es de: 



= A^ = 2,912381 
A y 1,853333 



(3) 



La ductilidad obtenida para este tipo de estructuras, |a = 1,57 , se considera muy baja, 
pues es inferior a \x = 2 que da la norma espanola NCSE-02 para este tipo de estructuras. 
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Para este edificio la reserva de resistencia es: R r = 1,94 . 



Si bien el valor de la ductilidad calculado para este edificio parece ser una sugerencia 
a revisar los valores del factor de ductilidad, tambien es cierto que los resultados muestran que 
este edificio exhibe una alta reserva de resistencia. 

En la Figura 6 se muestra la evolucion del indice de dano del edificio con forjados 
reticulares que cuantifica la perdida de rigidez del conjunto de elementos de la estructura que 
resisten la carga o conjunto de cargas que conducen al fallo de la misma. Dicho indice esta 
calculado mediante un programa de elementos finitos (PLCd, 1991-2007) con un modelo 
constitutive de dano y plasticidad y permite correlacionar el dano con los desplomes. El punto A 
indica la aparicion de las primeras micro fisuras en la estructura que se incrementan hasta 
alcanzar la rotula en las vigas (punto B) y, a partir de aqui, comienzan a aparecer las rotulas en 
pilares hasta que se alcanza el estado de umbral de colapso (punto C de la Figura 6). 



1 




0 0,3 0,6 0,9 1,2 1,5 1,8 2,1 2,4 2,7 3 

Desplome (%) 



Figura 6. Evolucion del indice de dano global en el portico exterior del edificio con forjados 
reticulares. 



La respuesta poco ductil del edificio con forjados reticulares puede atribuirse a la 
aparicion de rotulas plasticas en los puntos de transicion entre los abacos y los nervios del 
forjado del primer nivel. Debe recordarse que, al estar los elementos de los forjados sometidos 
a flexion inducida por cargas de gravedad, ademas de la debida a las fuerzas sismicas de 
calculo, las zonas que requieren un especial armado son las proximas a los nudos y al centra 
del vano, en donde frecuentemente se producen los mayores momentos. Se observa que es 
complicado poder llevar a cabo un eficiente confinamiento en la zona central de los forjados, lo 
que en cierto sentido explica su posible mecanismo fallo en el caso sismico y, por ende, el bajo 
nivel de ductilidad de la estructura. 



3.3. Respuesta no lineal del edificio con vigas planas 

Los edificios cuyos porticos tienen vigas planas reciben en la norma espahola NCSE-02 
una consideracion diferente a la del resto de los edificios porticados con vigas de canto, en lo 
que a ductilidad se refiere. Es notable la dificultad tecnica de llevar a cabo las disposiciones 
normativas de armado para garantizar una respuesta ductil de los elementos del portico. En la 
Figura 7 se muestra la respuesta global de la estructura hasta alcanzar el desplome ultimo 
(desplome previo al colapso total) que, junto con el valor del desplome de plastificacion, 
permite calcular el valor de la ductilidad estructural. 
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Figura 7. Curva de capacidad idealizada del portico exterior del edificio con vigas planas. 



La Figura 7 indica que el comportamiento se mantiene elastico hasta un valor 
relativamente bajo del coeficiente de cortante en la base (inferior a un valor de 0,10). Sin 
embargo, no es este valor el que se considera en el calculo de la ductilidad, ya que se utiliza el 
valor de plastificacion obtenido de la forma bilineal idealizada que se muestra en la misma 
Figura 7. En esta, tambien puede verse el calculo de la ductilidad del edificio a partir de la 
forma bilineal idealizada. Tal como puede verse, la ductilidad que se obtiene es de 1,55 lo que, 
sin duda alguna, llama la atencion debido a que el valor utilizado como factor de reduccion de 
respuesta, de acuerdo con las recomendaciones normativas, es de 2. Finalmente, en la Figura 
8 se muestra la evolucion del indice de dano global del portico en funcion del incremento del 
desplome que se produce al aplicar las cargas horizontales. En esta figura, el punto A, es el de 
aparicion de microfisuras en vigas, el B es el de formacion de rotulas plasticas en pilares y C el 
umbral de colapso. La evolucion del dano muestra que la rigidez del portico se degrada 
rapidamente, incluso para desplomes relativamente pequefios. Este hecho condiciona 
fuertemente la respuesta no ductil de estos edificios. 




Figura 8. Evolucion del indice de dano global en el portico. 



La respuesta del edificio con vigas planas muestra que la estabilidad de la estructura 
depende del fallo de sus vigas. Esta observacion es importante a la hora de decidir entre 
utilizar un sistema con vigas de canto o uno con vigas planas, puesto que este ultimo dispone 
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de una ductilidad inferior a la de proyecto y, por ende, de un factor de reduccion de respuesta R 
menor. 



3.4. Respuesta no lineal del edificio con porticos resistentes a momentos 

A continuation se muestran los resultados de la respuesta no lineal del edificio cuya 
estructura principal esta constituida por porticos resistentes a momento. Al igual que en los dos 
edificios anteriores, se ha analizado el portico externo del edificio de manera que sea posible 
comparar los resultados obtenidos en cada caso. En la Figura 9 se muestra la curva de 
capacidad en la que se puede apreciar que este tipo de edificio es capaz de sostener una 
respuesta ductil estable, evidenciada por alto valor del desplome ultimo. 

De la curva bilineal idealizada de la Figura 9 se obtiene una ductilidad de 5,15, que es 
un valor mayor que el considerado en el proyecto sismorresistente de la estructura, que es de 
4. Esto significa que los edificios con vigas de canto disponen de suficiente capacidad de 
respuesta ductil ante la action de fuerzas sismicas y ademas, de una adecuada reserva de 
resistencia. Finalmente, en la Figura 10 se muestra la evolucion del indice global de dafio para 
este tipo de edificio. 
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Figura 9. Curva de capacidad idealizada del portico externo del edificio con porticos 
resistentes a momentos. 
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La respuesta no lineal del edificio con porticos resistentes a momentos corresponde a 
la tfpica respuesta que exhiben los edificios de baja altura de hormigon armado, en los que el 
comportamiento general incluye la formacion de rotulas plasticas en las bases de los pilares del 
primer nivel. Este comportamiento general se debe a que no es facil proyectar edificios con 
pilar fuerte-viga debil, fundamentalmente por el predominio de las acciones de gravedad sobre 
las vigas que requieren secciones que, finalmente, resultan mayores que las de los pilares. 



4. COMPARACION DE LA RESPUESTA NO LINEAL DE LOS TRES EDIFICIOS 

Una vez obtenidos los resultados del analisis no lineal de los tres edificios, es 
interesante destacar los aspectos generales que caracterizan la respuesta de cada tipologia 
estructural. Para ello se han representado graficamente las respuestas de los mencionados 
edificios en la Figura 11, en la que puede observarse que los tres alcanzan un cortante en la 
base superior al de proyecto. Esto indica que los tres edificios satisfacen este objetivo inicial del 
proyecto, aunque es notoria la diferencia entre las reservas de resistencia de cada uno de 
ellos. Tambien se puede observar que los edificios porticados con vigas de canto son los que 
disponen de mayor capacidad ductil y que esta es varias veces superior a la de los edificios 
con vigas planas o con forjados reticulares. En resumen, los edificios porticados con vigas de 
canto son los unicos capaces de garantizar un comportamiento ductil y tener, al mismo tiempo, 
una reserva de resistencia satisfactoria. 




Figura 11. Comparacion de la respuesta no lineal de los tres edificios analizados. 



5. DEFECTOS DE PROYECTO CONCEPTUAL 

Son muchos, y de diferente naturaleza, los errores conceptuales que se pueden 
cometer en la fase de proyecto. A continuacion se describen y comentan algunos de los mas 
usuales, encontrados en edificios que, si bien no han sufrido danos por acciones sismicas, no 
estan exentos de sufrirlos durante terremotos futuros. 

La respuesta dinamica de los edificios puede alterarse cuando se incrementa la masa 
con la altura, caso que puede verse en la Foto 1 . Este hecho hace que en los niveles de mayor 
masa se incrementen las fuerzas de piso, ademas de producir el incremento en los efectos 
torsionales sobre los pilares. Otro caso de variation de la masa de la planta con la altura se 
aprecia en la Foto 2. Es de hacer notar que la excentricidad en planta no solo se incrementa 
con el crecimiento asimetrico de las dimensiones de esta, hasta cierto punto controlable en las 
etapas de proyecto de la edification, sino tambien por la concentration de cargas debida al 
almacenamiento de materiales. 

Otros errores de proyecto consisten en el uso de pilares cortos. El mecanismo de fallo 
por cortante que se produce en este caso se origina en el hecho de que los momentos en los 
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extremos del pilar se deben equilibrar con los cortantes en la longitud efectiva del pilar, que es 
una distancia reducida. Esto hace que los cortantes en los extremos del pilar superen a los de 
proyecto, correspondientes al resto de los pilares del mismo nivel, que tienen una longitud 
normal. En la Foto 3 se observa un pilar construido entre una pantalla, que es un elemento muy 
rfgido, y un forjado. Es evidente que los pilares con esta caracterfstica quedan expuestos al 
fallo por cortante. 



Foto 1. Irregularidad en planta: crecimiento de las plantas con la altura. 



Foto 2. Variacion de la masa de las plantas y discontinuidad en lineas resistentes verticales. 



Foto 3. Pilares cortos en un edificio. 



En el caso particular de los edificios con sistema estructural de forjados reticulares, el efecto de 
pilar corto puede ser inducido por los cambios de nivel de los forjados contiguos, agravandose 
el problema debido al bajo grado de redundancia de estos sistemas, que no permite transmitir 
las cargas bajo el hipotetico fallo de alguno de estos pilares. En la Foto 4 puede verse un caso 
extremo, en el que el pilar corto se reduce a una seccion sometida a corte directo. 
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Foto 4. Juntas pilar-forjado que inducen al comportamiento de corte directo. 



En la foto 4 tambien se puede apreciar otro aspecto que suele ser poco cuidado en el 
proyecto, como el de los nucleos de escaleras o de ascensores, que son de vital importancia a 
la hora de realizar cualquier labor de evacuacion de un edificio despues de una catastrofe. La 
mayor parte de los problemas que se producen se deben al hecho de que las subestructuras 
de circulacion suelen tener periodos de vibracion muy diferentes a las del resto del edificio, lo 
que puede conducir a desplomes seguidos de fallos localizados en los elementos verticales 
(pilares o pantallas a cortante). Otro de los efectos indeseados que puede ser originado por los 
nucleos de ascensor es que dichos nucleos, que se resuelven habitualmente mediante 
pantallas a cortante que tienen mayor rigidez que los pilares de la planta, pueden alterar la 
ubicacion de los centros de rigidez, provocando el incremento de la excentricidad de la planta y 
de la torsion global del edificio. 

El comportamiento inadecuado de los pilares puede provocar desplomes excesivos de 
estos a nivel de forjados. Este comportamiento se puede atribuir a la falta de arriostramiento 
lateral del pilar en una direccion, como se muestra en la Foto 5, en la que uno de los pilares se 
encuentra arriostrado mediante viga en una sola direccion, pero en la otra direccion la viga fue 
suprimida por requisitos arquitectonicos. 




Foto 5. Falta de arriostramiento del pilar al nivel de 
forjado. 



Un defecto grave consiste en no prever un sistema sencillo de transmision de las 
cargas, inclusive las sismicas, hasta la cimentacion. Por ejemplo, en la Foto 6 se muestra un 
pilar que no tiene continuidad en el nivel inferior, lo que puede llegar a crear problemas si el 
forjado que lo sostiene llega a danarse, originando un mecanismo que puede desencadenar el 
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colapso de la estructura entera. 




Foto 6. Discontinuidad de un pilar en las plantas inferiores de edificio. 



En la Foto 7 se muestra un edificio con planta baja debil. Esta situation se debe a la 
necesidad de contar con espacios abiertos de mayor altura en la planta baja, destinados 
principalmente al uso comercial. En la misma foto puede observarse como los paneles de 
mamposteria (elementos no estructu rales) no estan confinados por los elementos resistentes 
verticales que no son ni pilare ni pantallas, sino mas bien un intermedio entre ambos lo que, de 
por si, es un error. 



Foto 7. Planta baja libre en edificio. 

Frecuentemente se proyectan estructuras con pilares dimensionados con areas 
suficientes para controlar eficientemente los desplomes en una de las dos direcciones 
ortogonales de un edificio. Sin embargo, es necesario tener en cuenta la direccionalidad del 
sismo, por lo que se deben proyectar soportes verticales con suficiente rigidez en dos 
direcciones ortogonales. En la Foto 8 se muestra un edificio cuyos pilares muestran una gran 
diferencia de rigidez en una y otra direction. Tambien puede observarse otro defecto: no se 
respeta la continuidad de los ejes resistentes del edificio, como consecuencia de la ubicacion 
de los pilares en la planta sin que se permita la formation de porticos ortogonales. El hecho de 
desplazar los pilares por razones de espacio hace que algunas vigas se apoyen sobre otras 
vigas, lo que no permite la redistribution de momentos cuando plastifica alguno de los 
extremos de las vigas. Puede observarse como uno de los pilares esta colocado debajo de una 
vigueta y no debajo del “nudo” formado por la intersection de dos vigas planas. 



100 



Alex H. Barbat, Juan Carlos Vielma y Sergio Oiler 




Foto 8. Deficiente distribution de rigidez en los pilares. 



En la Foto 9 se muestra el inadecuado apoyo de un forjado sobre una viga metalica 
-siendo las viguetas simplemente apoyadas sobre la misma- lo que no asegura suficiente 
rigidez en piano del forjado. 




Foto 9. Apoyo inadecuado de un forjado. 
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Foto 10. Estructura singular, forjado postensado sobre juntas de neopreno. 

La construccion de estructuras singulares, llamadas asi por su particular configuracion, 
por solucionar casos poco usuales o por utilizar tecnicas constructivas novedosas, no deben 
realizarse al margen de los principios de proyecto sismorresistente. Por ejemplo, en la Foto 10 
se observa un edificio con una losa postensada que se sustenta sobre apoyos de neopreno 
colocados en la parte superior de los pilares de la planta baja. 

6. CONCLUSIONES 

En este trabajo se han podido verificar algunos aspectos del proyecto sismorresistente de 
edificios de hormigon armado, analizando su respuesta obtenida mediante un analisis estatico 
no lineal: 

• Se aplica el procedimiento de analisis estatico no lineal con control de fuerzas para el 
cual se fija como criterio de comportamiento para el desplome ultimo, o umbral de 
colapso de la estructura, en un punto en el que se alcanza un nivel preestablecido del 
fndice de dano. 

• Los desplomes de plastificacion de los porticos se obtienen mediante formas bilineales 
idealizadas, aplicando el procedimiento propuesto por Park (1988). Esta forma bilineal 
demuestra ser adecuada para igualar la energfa de la respuesta no lineal obtenida del 
analisis estatico no lineal y la respuesta no lineal idealizada, sin embargo no permite 
determinar valores adecuados del indice de dano D . 

• En cuanto a los resultados de la respuesta no lineal de los casos estudiados, los 
edificios porticados con vigas de canto disponen de suficiente ductilidad y reserva de 
resistencia para garantizar que un comportamiento estable, incluso para ductilidades 
superiores a las de proyecto. 

• Los edificios de ductilidad limitada tienen altas reservas de resistencia siempre y 
cuando, al ser sometidos a un terremoto, sean capaces de mantener la respuesta en el 
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rango elastico. 

• Los edificios con forjados reticulares pueden mejorar su respuesta ductil si se preven 
vigas de canto que conecten entre si, en las dos direcciones ortogonales, los extremos 
de los pilares de la misma planta. Esta mejora tambien se refleja en un incremento de 
la reserva de resistencia. 

• La respuesta ductil de los edificios con vigas planas no mejora sustancialmente si se 
mejora el confinamiento de sus elementos. En consecuencia, para garantizar su 
respuesta ductil, un adecuado confinamiento es requisito necesario pero no suficiente, 
siendo tambien necesaria una adecuada configuracion y un optimo dimensionamiento 
de las secciones transversales. El buen confinamiento solo es aprovechable en el caso 
de los edificios con configuracion y tipologia que permitan un comportamiento ductil, tal 
como son los edificios porticados con vigas de canto. 

• Se concluye que la respuesta ductil de los edificios porticados con vigas de canto 
resulta adecuada para edificios esenciales que requieren mantenerse en servicio luego 
de la ocurrencia de un terremoto, como es el caso de hospitales, estaciones de 
bomberos, etc. Sin embargo, es necesario que esta tipologia sea proyectada con un 
adecuado armado, prestando especial atencion a las zonas de los nudos y otras zonas 
sensibles a cortante y con aceros con caracteristicas especiales de ductilidad. 

• Es fundamental la cuidadosa aplicacion de las recomendaciones de disefio conceptual, 
de manera que los edificios proyectados no incorporen irregularidades que hagan que 
su respuesta sea compleja fragil e incluso inestable. 
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